AD-A267  984 


UMENTATION  PAGE 


form  Approved 
Of/B  .Vo.  0704-0188 


>or»  estimated  tc  a.eraq**  t  ni.ur  oer  "?sporve  including  me  lime  »o r  .. «;» .uK.n«.  vearcning  ev^tmg  aata 

dieting  ana  reviewing  me  -:olle<Tiono<  in«'>rmat.an  Sena  cowmen’s  regard,  r  g  th.v  buraen  estimate  or  anv  other  asoeci  o*  this 
ducmq  rn<>  ou^cen  to  w\ash»r  gton  »-reaaauarteri  Services.  O  mrrorare  ?or  information  Of^eration>  and  *teoort*.  I2i5  je^ler^on 
ana  tc  th*»0#*«f*0<  Management  ana  budget.  Paperwork  fledurtiOn  Proiect  (0704-0188).  Washington,  DC  20503 


4.  TITLE  AND  SUBTITLE 


3.  REPORT  TYPE  AND  DATES  COVERED 

THESIS/DISSERTATION 


5.  FUNDING  NUMBERS 


Modeling  Teh  Minimum  Energy  State  Of  The  Earth's 
Magnetotail 


6.  AUTHOR(S) 

Capt  Markus  Stephen  Sorrells 


7.  PERFORMING  ORGANIZATION  NAMt(S)  AND  ADDRESS(ES) 

AFIT  Student  Attending.  Utah  State  University 


8.  PERFORMING  ORGANIZATION  i 

REPORT  NUMBER  I 


AFIT/CI/CIA-  93-123  V 


9.  SPONSORING/MONITORING  AGtNCY  NAME(S)  AND  ADDRESS(ES) 

DEPARTMENT  OF  THE  AIR  FORCE 

AFIT/CI 

2950  P  STREET 

WRIGHT-PATTERSON  AFB  OH  45433-7765 


10.  SPONSORING  /  MONITORING 
AGENCY  REPORT  NUMBER 


12a.  DISTRIBUTION/  AVAILABILITY  STATEMENT 

12b.  DISTRIBUTION  CODE 

l 

i 

Approved  for  Public  Release  IAW  190-1 

Distribution  Unlimited 

MICHAEL  M.  BRICKER,  SMSgt,  USAF 

Chief  Administration 

1 

I  1 

I 

! 

! 

* 

» 

1 

1 

i 

1 

13.  ABSTRACT  (Maximum  200  words) 


$ 


DTIC 

jELECTE 
kAUG  17,1993 

1  B  . 


D 


93-18994 

IIIUKII 


14.  SUBJECT  TERMS 


15.  NUMBER  OF  PAGES 

130 


17.  SECURITY  CLASSIFICATION  118.  SECURITY  CLASSIFICATION  119.  SECURITY  CLASSIFICATION  I  20.  LIMITATION  OF  ABSTRACT 

OF  REPORT  |  OF  THIS  PAGE  I  OF  ABSTRACT  I 


NSN  7540  01-280-5500 


Standard  Form  298  (Rev  2  89) 

e.evnbro  by  ansi  mo  /(v-t 
JV810.’ 


MODELING  THE  MINIMUM  ENERGY  STATE  OF  THE 
EARTH’S  MAGNETOTAIL 


by 


Markus  Stephen  Sorrells 


A  diesis  proposal  submitted  in  pardal  fulfillment 
of  the  requirements  for  die  degree 

of 

MASTER  OF  SCIENCE 

in 

Physics 


Approved: 


Major  Professor 

Dr.  Kent  L  Miller  T  ~~ 
Committee  Member 


UTAH  STATE  UNIVERSITY 
Logan,  Utah 


ii 


ACKNOWLEDGMENTS 

I  would  like  to  thank  Dr.  W.  Farrell  Edwards  for  all  the  guidance  and  support  he  gave  me 
during  my  work  on  this  thesis.  I  would  also  like  to  thank  my  other  committee  members.  Dr.  Jan 
J.  Sojka  and  Dr.  Kent  L.  Miller,  for  all  their  assistance  throughout  the  entire  process  of  picking  a 
thesis  topic  to  actually  writing  and  defending  it  Also,  I  would  like  to  thank  Dr.  Anthony  T.  Y. 
Lui  for  providing  me  with  data  depicting  typical  magnetotail  profiles  during  solar  quiet  times  for 
the  research  in  this  thesis. 

I  also  wish  to  thank  Imad  Barghouthi  for  providing  assistance  with  the  computer  model 
programs  in  this  diesis.  Without  the  use  of  his  considerable  programming  experience  on  die 
VAX  computer  system.  I’m  convinced  I  would  still  be  programming  the  computer  today. 

I  give  a  special  thanks  to  my  family,  friends,  and  colleagues  for  their  encouragement,  moral 
support,  and  patience  as  I  worked  my  way  from  die  initial  proposal  writing  to  this  final 
document.  I  could  not  have  done  it  without  all  of  you. 


Markus  S.  Sorrells 


DTIC  QTTALTTY  tnspotbd  3 


Aooesslon  For 

HTIS  GRAAI 

DTIC  TAB 

□ 

Unannounced 
Justification _ 

□ 

By - 

Distribution/ 


Availability  Codes 


Avail  and/or 
Special 


CONTENTS 


Page 

ACKNOWLEDGMENTS  .  ii 

LIST  OF  TABLES  .  iii 

LIST  OF  FIGURES  .  iv 

CHAPTER 

I.  INTRODUCTION  . 1 

Earth’s  Magnetotail  .  1 

Magnetospheric/tail  Models  .  2 

Magnetospheric/tail  Physics  .  6 

Magnetotail  Energy  and  Storage  and  Dissipation  Mechanisms  .  16 

Magnetotail  Equilibrium/Steady  State  Condition  .  18 

II.  MINIMUM  TOTAL  ENERGY .  20 

Systems  in  Thermodynamic  Equilibrium  .  20 

Systems  Not  in  Thermodynamic  Equilibrium  .  21 

Calculus  of  Variations .  22 

m.  MAGNETOTAIL  MODEL  .  28 

Model  Description  ... . 28 

Magnetotail  Model’s  Total  Energy  Equation  .  30 

Equilibrium  Solutions  .  37 

Maxwell’s  Equilibrium  Equations  .  37 

Kinetic/Force  Equilibrium  Equations  .  41 

Boundary  Conditions  .  46 

Number  Density  Boundary  Conditions  .  46 

Electric  Field  Boundary  Conditions  .  47 

Magnetic  Field  Boundary  Conditions  . 47 

Particle  Conservation  Equations  . 47 

Minimizing  Model’s  Total  Energy  Equation  .  50 

Calculus  of  Variations  Calculations  .  50 

Integration  by  Parts  .  55 


IV 


IV.  MODEL  PROGRAM  .  59 

Program  Description  .... . . . . .  59 

DVCPR  Routine  .  59 

Dimensional  Transformation  . 60 

Model  Output  ... . 66 

V.  DATA  ANALYSIS  _ 82 

Observed  Magnetotail  Parameters  _ ...... _ ... _ ....... _ .... _ ... _ ....  82 

Comparison  to  Model  Output  . 82 

VI.  CONCLUSIONS  _  89 

REFERENCES  _  91 

APPENDIXES  .  93 

A.  IMSL  DVCPR  Program  Information  94 

B.  Minimum  Total  Energy  Program  . 99 

C.  Model  Data  _ _ _ ..................... - - — - - — ......... - - - - -  116 

D.  Bibliography  . . .. .  126 


V 


LIST  OF  TABLES 

Table  Page 

I  Potential  Energy  Functions  . . . . .  20 

II  Scaling  Data  Corresponding  to  Fig.  (24)  .  69 

III  Scaling  Data  Corresponding  to  Fig.  (25)  . 70 

IV  Scaling  Data  Corresponding  to  Fig.  (26)  .  71 

V  Scaling  Data  Corresponding  to  Fig.  (27)  .  72 

VI  Scaling  Data  Corresponding  to  Fig.  (28)  .... . . . 73 

VII  Scaling  Data  Corresponding  to  Fig.  (29)  .  75 

VIII  Scaling  Data  Corresponding  to  Fig.  (30)  .  76 

IX  Scaling  Data  Corresponding  to  Fig.  (31)  . .  77 

X  Scaling  Data  Corresponding  to  Fig.  (32)  _ ... _ _ _ ..... _ ............... . .  78 

XI  Scaling  Data  Corresponding  to  Fig.  (33)  . 79 

XII  Scaling  Data  Corresponding  to  Fig.  (34)  .  80 

XIU  Profile  Comparison  Summary  Between  Model  and  Magnetotail  Systems  .  86 


vi 


LIST  OF  FIGURES 

Figure  Page 

1  The  original  undisturbed  model  with  magnetic  field  lines 

extending  into  a  vacuum  _ _ _ ......... . .  2 

2  The  closed  magnetospheric  model  in  which  the  solar  wind 

distorts  die  earth’s  magnetic  field  .  3 

3  Wavy  structure  of  the  heliomagnedc  current  sheet  due 
to  die  current  sheet’s  inclination  to  the  sun’s 

rotational  equator  . .  4 

4  An  open  magnetospheric  model  in  which  the  terrestrial 

and  interplanetary  field  lines  connect  .  5 

5  The  interaction  of  the  earth’s  dipole  magnetic  field  with 
(a)  southern  component  of  the  IMF  and  (b)  the  northern 

component  of  the  IMF  . . .. .  6 

6  Diagram  showing  die  structure  of  die  earth’s 

magnetosphere  system  . . ....... -  7 

7  Schematic  showing  how  the  Lorentz  force  creates 

the  magnetopause  current  .  8 

8  Flow  patterns  of  the  two  principal  current  systems 
which  determine  the  configuration  of  the 

magnetosphere/tail  .................................... _ ...... — - ...... - ...... -  8 

9  Cross  section  of  the  earth’s  magnetotail  showing 
the  primary  directions  of  die  magnetopause 

and  crosstail  currents  .  9 

10  A  view  of  the  dayside  plasmasphere  showing  die 
creation  of  a  westward  ring  current  by 

particle  drift  .......................... _ ..... ..... .................... — . — ......................... - - -  10 

11  Field  aligned  currents  - - -  10 

12  Schematic  showing  magnetic  merging  process  .  11 

13  Diagram  showing  access  of  solar  wind  particles 
through  die  northern  and  southern  portions 

of  the  magnetosphere  . . . . .  13 


vii 

14  Access  through  the  LLBL  .  13 

15  Sketch  of  the  equatorial  section  of  the  earth’s 
magnetosphere  looking  from  above  die  north 

pole  showing  idealized  plasma  flow  .  13 

16  Schematic  showing  the  electric  field  set  up  as  a 

result  of  die  Lorentz  force  and  £“-9x6  .  15 

17  Schematic  showing  the  earthward  9x6  drift  .  15 

18  Schematic  showing  the  equatorial  9x6  drift  .  16 

19  (a)  Schematic  of  how  solar  wind  energy  was  thought 
to  be  transported  to  die  earth,  (b)  Current  view 

of  how  solar  wind  energy  is  thought  to  get  to  the  earth .  17 

20  Diagram  showing  how  the  magnetosphere's  energy 
storing  and  release  processes  resemble  a 

faucet  drip  . 17 

21  Diagram  showing  that  a  system  can  have  many 

equilibrium  states  . .  25 

22  Since  die  value  for  U(  a  >  0  )  will  always  be 
greater  than  U(  a  =  0  ),  we  can  be  assured 
that  we  will  find  the  equilibrium  condition 
with  the  least  possible  energy  when  a  goes 

to  zero  .  25 

23  Drawing  showing  the  layout  of  die  magnetotail  model  .  28 

24  N-S  profiles  of  number  densities  (n^  and  n2),  x-component 
of  the  magnetic  field  (Bx),  and  the  z  component  of  the 
electric  field  (Eg)  generated  from  die  minimum  total 

energy  program  . 69 

25  N-S  profiles  of  number  densities  (nj  and  n2),  x-component 
of  the  magnetic  field  (Bx),  and  die  z  component  of  the 
electric  field  (Eg)  generated  from  the  minimum  total 

energy  program  _  70 

26  N-S  profiles  of  number  densities  (nj  and  n2),  x-compoi 

of  the  magnetic  field  (BJ,  and  the  z  component  of  die 
electric  field  (Eg'  ‘ '  * 

energy  program 


71 


viii 


27  N-S  profiles  of  number  densities  (n^  and  n 2),  x -component 
of  the  magnetic  field  (Bx),  and  the  z  component  of  the 
electric  field  (Eg)  generated  from  th'  minimum  total 

energy  program  .  72 

28  N-S  profiles  of  number  densities  (nj  and  n2),  x-component 
of  the  magnetic  field  (Bx),  and  die  z  component  of  the 
electric  field  (Ex)  generated  from  the  minimum  total 

energy  program  . 73 

29  N-S  profiles  of  number  densities  (nj  and  n2),  x-component 
of  the  magnetic  field  (Bx),  and  the  z  component  of  tire 
electric  field  (Eg)  generated  from  the  minimum  total 

energy  program  .  75 

30  N-S  profiles  of  number  densities  (nj  and  n2),  x-component 
of  the  magnetic  field  (Bx),  and  the  z  component  of  the 
electric  field  (Eg)  generated  from  the  minimum  total 

energy  program  . .. .  76 

31  N-S  profiles  of  number  densities  (n^  and  n2),  x-component 
of  the  magnetic  field  (Bx)f  and  the  z  component  of  die 
electric  field  (Eg)  generated  from  the  minimum  total 

energy  program  . . .  77 

32  N-S  profiles  of  number  densities  (n^  and  t^),  x-component 
of  die  magnetic  field  (Bx),  and  the  z  component  of  die 
electric  field  (Eg)  generated  from  the  minimum  total 

energy  program  .  78 

33  N-S  profiles  of  number  densities  (nj  and  n2),  x-component 
of  the  magnetic  field  (Bx),  and  the  z  component  of  die 
electric  field  (Eg)  generated  from  the  minimum  total 

energy  program  ..... .  79 

34  N-S  profiles  of  number  densities  (nj  and  n2),  x-component 
of  the  magnetic  field  (Bx),  and  the  z  component  of  the 
electric  field  (Eg)  generated  from  die  minimum  total 

energy  program  ............ _ .............................................................. _ ......... _ _ _ _ _ _ _  80 

35  North-south  profiles  of  magnetic  field  (B),  number 
density  (n),  plasma  bulk  flow  speed  (v),  and 

ion  temperature  (Tp  in  the  midnight  region 

of  the  magnetotail .  83 

36  Number  dinsity  (n)  profile  from  Fig.  (35) 

using  a  linear  number  density  scale  ................. . . .  83 

37  Plasma  bulk  flow  speed  (v)  profile  from  Fig.  (35) 
using  a  linear  velocity  scale 


84 


LX 


38  Ion  temperature  (Tj)  profile  from  Fig.  (35) 
using  a  linear  temperature  scale  _ 


84 


X 


ABSTRACT 

Modeling  die  Minimum  Energy  State  of  the 
Earth’s  Magnetotail 

by 

Markus  S.  Sorrells,  Master  of  Science 
Utah  State  University,  1993 

Major  Professor  Dr.  W.  Farrell  Edwards 
Department:  Physics 

For  a  system  that  remains  in  thermodynamic  equilibrium,  stable  equilibria  are  determined  by 
minimizing  an  appropriate  potential  energy  function  such  as  the  Gibb’s  free  energy.  However, 
when  a  system  does  not  remain  in  thermodynamic  equilibrium  (Le.  a  radiating  system),  one 
cannot  use  a  potential  function  to  derive  a  state  of  stable  equilibrium.  If  we  assume  this  is  die 
case  for  the  earth’s  magnetotail,  then  we  must  come  up  with  another  method  to  determine  stable 
equilibrium.  We  conjecture  that  this  will  be  a  minimum  total  energy. 

By  combining  Poynting*s  theorem  with  die  Loren tz  force  equation  and  Maxwell’s  equations, 
one  can  account  for  the  energy  advections  into  and  out  of  a  system  containing  a  fully  ionized 
plasma,  as  well  as  for  the  energy  fluctuations  caused  by  changing  electric  and  magnetic  Adds. 

By  minimizing  the  total  energy  equation  representing  an  idealized  model  of  the  earth’s 
magnetotail,  one  should  be  able  to  calculate  its  stable  equilibrium  state. 

Although  the  idealized  minimum  energy  model  presented  in  this  paper  depicts  many  trends 
observed  in  the  earth’s  magnetotail,  this  thesis  concludes  that  more  modifications  are  needed 
before  it  becomes  a  useful  tool  for  analyzing  the  earth’s  magnetotail.  (140  pages) 


CHAPTER! 


INTRODUCTION 
Earth’s  Magnetotail 

This  thesis  attempts  to  develop  a  model  to  determine  if  a  minimum  energy  equilibrium  state 
establishes  itself  within  the  earth’s  magnetotail.  During  quiet  conditions,  the  earth’s  magnetotail 
receives  and  stores  energy  from  the  solar  wind.  At  periodic  intervals,  the  magnetotail  releases 
all  or  a  portion  of  this  stored  energy  into  the  earth’s  ionosphere  and  interplanetary  space  in  the 
form  of  geomagnetic  storms,  substorms,  and  plasmoid  ejections.  At  this  point  in  time,  no  one  is 
quite  certain  how  or  why  the  magnetotail  periodically  releases  its  energy  although  there  are  a 
couple  of  theories  available.  One  theory  suggests  that  the  solar  wind  may  deposit  energy  into 
the  magnetotail  at  a  rate  faster  than  the  magnetotail  can  absorb  it,  causing  the  excess  energy  to 
be  discarded  during  the  energy  storage  cycle.  Another  theory  states  that  the  magnetotail  has  a 
finite  energy  storage  capacity  and  when  this  critical  level  is  reached,  some  unknown  triggering 
mechanism  causes  the  magnetotail  to  release  its  energy.  Whatever  the  case,  this  energy  release 
is  probably  die  result  of  the  magnetotail  trying  to  establish  an  equilibrium  state.  It  seems 
reasonable  that  die  magnetotail  would  ultimately  be  trying  to  reach  an  equilibrium  state  having 
die  lowest  possible  total  energy.  This  diesis  attempts  to  find  out  if  the  magnetotail  releases 
enough  of  its  energy  at  one  time  to  ever  reach  this  minimum  total  energy  state. 

The  minimum  energy  model  developed  in  this  diesis  calculates  die  minimum  energy  profiles  of 
key  parameters  of  the  magnetotail  such  as  magnetic  held  strength,  ion  and  electron  number 
densities,  temperatures,  and  bulk  flow  velocities.  By  comparing  the  model’s  profiles  with  those 
observed  in  the  magnetotail  at  the  end  of  an  energy  release  cycle,  one  should  be  able  to 
ascertain  whether  die  model  works  or  not  Certain  insights  to  die  earth’s  magnetotail  energy 
cycle  may  be  gained  if  the  model  adequately  describes  die  above  situation.  Before  getting  into 
the  minimum  energy  model,  a  brief  background  into  die  magnetospheric/tafl  physics  and 
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structure  will  be  discussed  in  this  chapter. 

Magnetospheric/taQ  Models 

Prior  to  the  advent  of  space  exploration,  scientists  had  to  rely  solely  on  ground  measurements 
to  observe  and  come  up  with  theories  to  explain  the  structure  of  die  earth’s  magnetic  field.  As  a 
result,  early  scientists  pictured  die  earth’s  magnetic  field  as  a  perfect  magnetic  dipole  (Fig- 1). 

In  1931,  Chapman  and  Ferraro  brought  the  acceptance  of  a  solar  wind  into  die  scientific 
community  (Akasofu,  1981).  Solar  wind  plasma  originates  from  die  sun’s  corona  and  flows 
radially  outward  at  ’supersonic"  speeds  that  vary  between  300  to  800  kilometers  per  second 
(Gosling,  1984).  As  the  solar  wind  reaches  die  earth,  it  has  a  typical  density  of  about  10 
particles  per  cubic  centimeter  (Gosling,  1984)  as  compared  to  die  3  x  10^  particles  per  cubic 
centimeter  found  in  die  atmosphere  at  the  earth’s  surface.  During  solar  quiet  times,  the  energy 
of  the  solar  wind  particles  is  normally  a  few  eVs. 


FIG.  1.  The  original  undisturbed  model  with  magnetic  field  lines  extending  into  a  vacuum  [Gosling, 
1984]. 
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Scientists  had  to  take  into  account  the  effects  of  the  solar  wind  how  against  the  earth's 
magnetic  held.  They  came  up  with  what  is  now  called  the  closed  magnetospheric  model  in 
which  magnetic  held  lines  left  one  hemisphere  of  the  earth  and  returned  to  the  other 
hemisphere  and  were  modified  by  the  solar  wind  how  into  hie  shape  depicted  in  Fig.  2. 

This  model  worked  relatively  well  for  describing  some  of  the  magnetospheric  processes  that 
were  observed  before  the  advent  of  space  exploration.  In  1955,  a  British  geophysicist  by  the 
name  of  James  W.  Dungey  attempted  to  apply  hydromagnetic  theory  to  calculate  and  determine 
the  shape  of  the  magnetotail  (Carovillano  et  a].,  1967).  As  it  turned  out,  his  calculations  worked 
fairly  well  on  the  dayside  portion  of  the  magnetotah  but  it  predicted  that  the  magnetotail  would 
close  off  relatively  dose  behind  the  earth  with  something  like  a  Mach  angle.  Of  course  current 
observations  do  not  support  this  theory.  As  more  magnetospheric  processes  were  observed  with 
the  use  of  satellites,  it  became  apparent  that  the  dosed  magnetospheric  model  was  not  suffident. 
For  example,  it  could  not  explain  how  solar  wind  particles  gained  relatively  easy  access  to  the 
earth’s  polar  regions  nor  could  it  explain  the  sharp  outer  boundary  layer  of  the  plasma  sheet 
which  separated  the  anti-paralld  magnetic  fidd  lines  of  the  magnetotail. 


Closed  Magnetosphere 


FIG.  2.  A  closed  magnetospheric  model  in  which  the  solar  wind  distorts  the  earth’s  magnetic  field. 
All  the  earth’s  magnetic  lines  start  and  end  on  the  earth.  (For  simplicity,  the  tilt  of  the  earth’s  axis 
is  not  shown)  (Gosling,  1984]. 


4 


In  the  early  1960’s  scientists  discovered  that  the  solar  wind  contained  a  magnetic  field  Since 
the  sun  has  a  dipolar  magnetic  structure,  the  solar  wind  carries  "frozen*  remnants  of  the  sun’s 
magnetic  field  within  its  flow.  This  "frozen  in*  magnetic  field  is  known  as  die  Interplanetary 
Magnetic  Field  or  IMF.  Due  to  the  dipole  nature  of  the  sun,  die  IMF  is  either  a  sunward  or  anti- 
sunward  orientation.  A  heliomagnetic  current  sheet  separates  these  oppositely  aligned  magnetic 
fields,  which  prevents  them  from  making  contact  and  mutually  annihilating  one  other.  Due  to 
die  rotation  of  the  sun  and  the  current  sheet’s  7°  inclination  (Tascione,  1988)  to  the  solar 
rotational  equator,  the  current  sheet  incurs  a  wavy  structure  depicted  in  Fig.  3.  Because  of  the 
wavy  nature  of  this  current,  the  earth’s  orbit  passes  through  it  at  least  twice  during  each  solar 
rotation.  As  the  earth  passes  through  this  structure,  the  IMF  will  either  have  a  northern  or 
southern  component  to  it  depending  on  whether  the  earth  is  above  or  below  the  heliomagnetic 
current  sheet.  The  acceptance  of  a  magnetic  field  (or  IMF)  in  the  solar  wind  led  scientists  to 
develop  an  open  magnetospheric  model.  In  this  model,  some  of  the  magnetic  field  lines  flow 
from  the  earth  out  into  interplanetary  space  instead  of  flowing  back  towards  the  earth  at  the 


FIG.  3.  Wavy  structure  of  the  heliomagnetic  current  sheet  due  to  the  current  sheet’s  inclination  to 
the  sun’s  rotational  equator  [after  National  Research  Council,  1981]. 
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opposite  hemisphere  (Pig.  4).  This  model  is  able  to  explain  many  of  die  phenomena  observed 
within  die  magnetosphere/tail  that  die  dosed  model  could  not,  such  as  how  solar  wind  plasma 
and  energy  are  able  to  move  into  and  through  die  magnetosphere  (Gosling,  1984).  In  fact,  this 
model  can  also  attempt  to  explain  the  dependence  of  die  earth’s  magnetic  activity  on  die  IMF. 

When  the  IMF  points  southward,  as  seen  in  Fig.  5(a),  its  magnetic  field  can  easily  merge  with 
the  earth’s  own  fidd.  According  to  the  model,  magnetic  merging  is  one  of  the  primary  processes 
by  which  the  solar  wind  transfers  a  portion  of  its  mass/energy  to  the  earth’s  magnetosphere/tail. 
As  a  result,  the  earth’s  magnetosphere/ tail  is  rdatively  active  during  southward  IMF. 

However,  magnetic  merging  occurs  less  often  when  the  IMF  points  northward  (Fig.  5(b)). 

This  means  the  amount  of  solar  wind  mass/energy  transferred  into  the  earth’s 
magnetosphere/ tail  is  reduced.  As  a  result,  die  earth’s  magnetosphere/ tail  experiences  less 
geomagnetic  activity  when  the  IMF  is  pointing  northward  rather  than  southward. 


Open  MtQneioiphtre 


FIG.  4.  An  open  magnetospheric  model  in  which  die  terrestrial  and  interplanetary  fidd  lines 
connect.  (Tilt  of  the  earth’s  axis  is  not  shown)  (Gosling,  1984]. 
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FIG.  5.  The  interaction  of  the  earth’s  dipole  magnetic  field  with  (a)  southern  component  of  the  IMF 
and  (b)  the  northern  component  of  the  IMF. 


Even  with  all  its  good  points,  the  open  magnetospheric  model  does  not  explain  every  situation 
within  the  magnetosphere/tail  so  most  scientists  today  use  a  combination  of  the  open  and  dosed 
models  to  describe  magnetospheric  processes. 

Magnetospheric/tail  Physics 

The  "supersonic*  plasma  of  the  solar  wind  encounters  the  presence  of  the  earth’s 
magnetosphere  at  a  boundary  referred  to  as  the  "bow  shock"  (Fig.  6).  This  encounter  is 
analogous  to  a  shock  wave  created  by  an  object  moving  at  supersonic  speeds  through  a  fluid 
medium.  The  region  directly  behind  the  bow  shock  is  called  the  "magnetosheath"  and  forms  a 
physical  boundary  between  the  solar  plasma  at  the  bow  shock  and  the  magnetospheric  plasma  at 
the  magnetopause.  As  the  solar  wind  particles  continue  to  flow  past  the  earth’s  northward 
pointing  magnetic  field,  they  are  acted  upon  by  the  Lorentz  force  P  ■}(  V  xfi)  which 
deflects  the  positive  ions  towards  the  dusk  side  of  the  earth  and  the  negative  electrons  towards 
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MEDIUM  MAGNETOSHEATH 


FIG.  6.  Diagram  showing  the  structure  of  the  earth's  magnetosphere  system  [National  Research 
Council,  1981]. 


the  dawn  side  (Fig.  7).  This  force  sets  up  a  magnetopause  current  that  follows  die  path  taken 
by  the  positive  ions  and  flows  around  the  outside  of  the  magnetosphere  and  eventually  merges 
with  die  magnetotail  current,  which  flows  around  die  exterior  of  the  magnetotail  from  dusk  to 
dawn  as  depicted  in  Fig.  8.  The  magnetotail  current  then  merges  with  the  current  or  plasma 
sheet  current  (Fig.  9),  which  flows  from  dawn  to  dusk  in  the  interior  of  die  magnetotail.  The 
magnetotail  and  plasma  sheet  currents  are  enhanced  by  the  earth’s  magnetic  field,  which  points 
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FIG.  7.  Schematic  showing  how  the  Loren tz  force  creates  the  magnetopause  current  [from 
McPherron,  1991]. 


FIG.  8.  Flow  patterns  of  the  two  principal  current  systems  which  determine  the  configuration  of  the 
magnetosphere/tail  [adapted  after  Axford,  1965]. 


sunward  in  the  northern  half  of  the  magnetotail  and  anti-sunward  in  the  southern  half  (Fig.  6). 
The  ring  current,  which  constitutes  the  earthward  most  regions  of  the  plasma  sheet,  is  due  to 
drift  of  charged  particles  across  the  earth's  magnetic  field  gradient  The  drifts  are  charge 
dependent  so  that  positive  particles  drift  westward  and  negative  particles  drift  eastward 
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FIG.  9.  Cross  section  of  the  earth's  magneto  tail  showing  the  primary  directions  of  the  magnetopause 
and  crosstail  currents  [following  Axford  et  al.,  1965]. 


(Fig.  10)  producing  a  net  westward  current  (Williams,  1987). 

The  field  align  currents  (FACs),  also  known  as  Birkeland  currents,  connect  the  earth’s 
ionosphere  to  the  magnetosphere  and  interplanetary  space.  The  net  current  usually  flows  into 
die  ionosphere  at  die  morning  sector,  across  the  polar  cap  (Auroral  electrojet),  and  out  of  the 
:  osphere  at  the  evening  sector  (Fig.  11).  The  magnetotail  forms  when  die  earth's  magnetic 
field  lines  are  stretched,  as  far  as  1000  earth  radii  or  downstream  by  the  solar  wind.  The 
process  starts  when  die  IMF  and  the  earth's  magnetic  field  merge  and  connect  on  the  sunward 
side  of  die  earth  near  die  magnetopause.  The  solar  wind  is  then  able  to  drag  the  interconnected 
magnetic  field  lines  from  die  dayside  portion  of  the  magnetosphere  down  die  magnetotail  to  the 
point  where  they  reconnect  (Fig.  12)  (Sibeck,  1990). 

As  stated  above,  hydromagnetic  theory  was  not  able  to  explain  die  long  nature  of  die 
magnetotail.  The  problem  was  that  it  assumed  the  solar  wind  slid  smoothly  along  the  boundary 
of  the  magnetotail  without  agitating  it  Scientists  soon  concluded  that  there  had  to  be  some 
force  exerted  on  the  tail  to  allow  it  to  be  dragged  out  for  such  a  distance.  One  explanation  is 
that  there  must  be  a  pressure  pj  within  the  tail  region  which  tends  to  push  out  with  just  enough 


FIG.  11.  Field  aligned  currents  [from  Clauer  and  McPherron,  1974] 
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FIG.  12.  Schematic  showing  magnetic  merging  process  [Aeronomy  Lecture,  1992]. 

pressure  to  equal  die  inward  magnetic  tension  (or  pressure)  of  die  tail  (Carovillano  et  al.,  1967). 
This  effect  is  expressed  in  Eq.  (1). 
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However,  there  is  also  a  pressure  outside  the  tail  P2>  which  combines  with  the  inward  magnetic 
tension  of  the  tail  to  push  against  die  internal  pressure.  Therefore,  in  order  for  a  transverse 
equilibrium  condition  to  exist,  the  total  pressure  inside  the  tail  p  j  must  balance  both  the  inward 
magnetic  tension  of  die  tail  plus  the  pressure  outside  die  tail  p2-  This  pressure  balance  is 
expressed  in  Eq.  (2). 
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Enhanced  mixing  of  the  solar  wind  plasma  with  the  magnetopause  boundary  is  also  used  to 
explain  die  magnetosphere's  long  tail.  It  is  thought  that  this  mixing  could  be  due  to  surface 
instabilities  on  die  magnetotail  boundary  or  other  unknown  mechanisms.  Whatever  die  cause, 
this  mixing  would  allow  particles  from  the  solar  wind  to  penetrate  into  the  tail  at  large 
distances.  The  inertia  of  this  penetration  could  then  be  the  mechanism  that  stretches  die 
magnetotail  out  to  the  long  distances  that  we  observe.  The  solar  wind  penetration  can  also  be 
the  source  of  die  tension  along  die  magnetic  field  lines.  Calculations  have  shown  that  only 
1/10*  of  the  total  solar  wind  momentum  is  required  to  produce  die  necessary  shape  observed 
and  tension  required  within  the  earth's  magnetotail  (Carovdlano  et  al.,  1967).  It  is  these 
calculations  which  have  caused  the  latter  explanation  to  be  widely  accepted. 

In  order  for  solar  wind  plasma  to  interact  within  the  earth’s  magnetic  field,  it  must  first  get 
into  the  magnetospheric/tail  system.  One  way  this  can  happen  is  through  the  magnetic  merging 
of  the  earth’s  magnetic  field  with  an  IMF  as  described  above  CTascione,  1988).  As  the 
interconnected  magnetic  fields  are  pushed  along  by  die  solar  wind,  some  of  die  solar  wind 
particles  are  dragged  along.  When  die  magnetic  field  lines  reconnect,  a  portion  of  the  field  line 
gets  propelled  earthward  by  the  inward  tension  of  the  magnetic  field,  thereby  injecting  some  of 
the  trapped  solar  wind  particles  into  die  earth’s  magnetosphere  (Fig.  12). 

Solar  wind  particles  can  also  get  into  the  earth’s  magnetotail  through  diffusion  across  the 
magnetosheath  boundary  layer.  A  portion  of  the  solar  wind  particles  get  through  the  bow  shock 
and  are  able  to  diffuse  through  the  magnetopause  as  they  flow  along  it.  Particles  flowing  across 
die  top  and  bottom  of  die  magnetosphere  must  travel  through  the  magnetotail’s  lobes  before 
they  reach  die  plasma  sheet  (Fig.  13).  Particles  flowing  along  the  equatorial  region  of  the 
magnetosphere,  also  known  as  die  Low-Latitude  Boundary  Layer  (LLBL),  are  able  to  diffuse 
through  die  magnetopause  and  enter  die  plasma  sheet  from  die  sides  (Fig.  14)  (Lundin  et  al., 
1991).  Once  die  solar  wind  plasma  is  inside  the  magnetosphere/tail  system,  it  initially  flows 
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Lobe/Plasma  Mantle  Source 


FIG.  13.  Diagram  showing  access  of  solar  wind  particles  through  the  northern  and  southern 
portions  of  the  magnetosphere  [Pilipp  and  Morfil,  1978,  Cowley  and  Southwood,  1980]. 


Low-Latitude  Boundary  Layer  Source 


FIG.  14.  Access  through  die  LLBL  [Heikkila,  1982,  Eastman  et  al.,  1985]. 


FIG.  15.  Sketch  of  die  equatorial  section  of  the  earth’s  magnetosphere  looking  from  above  the  north 
pole  showing  idealized  plasma  flow  [after  Axford,  1964]. 
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anti-sunward  in  die  magnetosheath.  Most  of  the  plasma  within  the  magnetosphere  moves  into 
the  plasma  sheet  located  in  the  equatorial  region  of  the  magnetotail.  Once  in  the  plasma  sheet, 
the  plasma  flow  is  generally  sunward  (Fig.  15)  from  die  earthward  side  of  the  magnetic  field 
reconnection  point  (about  100  Rg).  The  plasma  flow  within  the  plasma  sheet  is  due  to  either 
magnetic  field  line  reconnections  or  die  inward  solar  wind  pressure  on  the  earth’s  magnetic 
field. 

As  die  plasma  begins  to  flow  earthward  from  the  magnetic  field  reconnection  site,  several 
things  happen  to  it  First,  the  plasma  again  experiences  the  Loren tz  force  P  «  j  (  ?  x  8  ), 
which  forces  die  positive  ions  to  flow  towards  the  dawnside  of  die  magnetotail  and  the  negative 
electrons  to  flow  towards  the  duskside.  The  net  effect  of  this  action  is  to  set  up  an  electric  field 
across  die  magnetotail  which  points  from  dawn  to  dusk  (Fig.  16).  Another  way  of  looking  at 
this  is,  since  the  plasma  can  be  considered  a  collisionless  medium,  its  motion  sets  up  an  electric 
field  which  satisfies  the  equation  fi  «  -  9  x  S. 

Since  this  plasma  can  be  considered  collisionless,  it  also  experiences  an  fi  x  fi  (electric  field 
vector  crossed  with  the  magnetic  field  vector)  drift  (Hones,  1986).  Near  the  earthward  side  of 
die  reconnection  site,  the  earth’s  magnetic  field  has  a  strong  northerly  component  to  it  As  a 
result,  the  £  x  3  drift  directs  the  plasma  towards  the  earth  (Fig.  17).  Midway  between  die 
reconnection  site  and  the  earth,  the  fi  x  8  drift  forces  the  plasma  towards  the  equatorial 
regions  of  the  magnetotail  due  to  the  anti-parallel  orientation  of  earth’s  magnetic  fields  at  this 
point  (Fig.  18).  As  a  result  of  this  action,  the  bulk  of  die  magnetotail’s  plasma  and  energy  is 
tied  up  within  a  relatively  thin  equatorial  region  called  die  central  plasma  sheet  This  is  where 
key  processes  such  as  energy  storage  and  dissipation  usually  take  place. 
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FIG.  16.  Schematic  showing  die  electric  field  set  up  as  a  result  of  die  Lorentz  force  and 
-  ?xS- 
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FIG.  17.  Schematic  showing  die  earthward  fi  x  6  drift 
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FIG.  18.  Schematic  showing  the  equatorial  fi  x  6  drift. 

Magnetotail  Energy  Storage  and  Dissipation  Mechanisms 

The  solar  wind/magnetospheric  interaction  basically  converts  die  solar  wind’s  kinetic  energy 
into  die  electrical  and  magnetic  energies  observed  within  the  magnetosphere/tail.  The  solar 
wind  puts  energy  into  the  earth’s  magnetosphere/tail  system  at  a  approximate  rate  of  1019 
ergs/sec  (Lanzerotti  and  Krimigis,  1985).  It  appears  that  the  earth’s  magnetospheric  system  has 
two  models  (Fig.  19)  to  describe  the  way  it  processes  the  energy  it  receives  horn  the  solar  wind. 
One  way  is  by  direct  dissipation  of  the  solar  wind  energy  into  die  earth’s  upper  atmosphere 
where  it  manifests  itself  in  the  form  of  auroral  storms  (i.e.  northern  lights).  This  process  is 
called  the  direct  driven  model  and  is  depicted  in  Fig.  19(a)  (Akasofu,  1987). 

The  other  way  die  magnetospheric  system  handles  the  solar  wind’s  energy  is  by  storing  it  and 
releasing  it  periodically.  This  is  called  the  driven  reconnection  model  and  is  depicted  in  Fig. 
19(b)  (Akasofu,  1987).  This  model  is  analogous  to  a  dripping  faucet  (Fig.  20)  where  'solar  wind 
plasma  seeps  into  the  magnetosphere  all  along  its  boundary,  accumulating  in  die  tail  until  a 
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FIG.  19.  (a)  Schematic  of  how  solar  wind  energy  was  thought  to  be  transported  to  the  earth,  (b) 
Current  view  of  how  solar  wind  energy  is  thought  to  get  to  die  earth  [Akasofu,  1987]. 


FIG.  20.  Diagram  showing  how  the  magnetosphere’s  energy  storing  and  release  processes  resemble 
a  faucet  drip  [Gosling,  1984]. 
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portion  breaks  off  like  a  swollen  drop  of  water  from  a  faucet"  (Gosling,  1984,  p.  33).  In  this 
process,  a  portion  of  the  stored  energy  is  shed  while  another  portion  is  directed  towards  the 
earth’s  poles.  The  shedding  process  has  been  directly  observed  in  comets  but  not  in  the  earth’s 
magnetotail. 

Most  of  the  solar  wind’s  energy,  and  particles,  appears  to  be  stored  within  die  central  plasma 
sheet  of  the  magnetotail.  Observations  indicate,  however,  that  plasma  sheet  energies  are  in  the 
keV  range  while  solar  wind  particles  are  only  in  the  eV  range.  That  means  there  must  be  some 
mechanism  within  the  magnetotail  which  somehow  accelerates  the  low  energy  solar  wind 
particles  to  die  higher  plasma  sheet  energies  observed.  Magnetic  field  reconnection  may  be  one 
possible  mechanism.  Scientists  have  predicted  that  when  two  anti-parallel  field  lines  of  the 
magnetotail  lobes  come  into  contact,  they  "explosively*  cancel  each  other  out  with  a  resulting 
release  of  energy  in  the  5  kev  range.  This  energy  might  be  transformed  into  the  high  plasma 
kinetic  energy  observed  within  the  central  plasma  sheet  (Tascione,  1988). 

Magnetotail  Equilibrium/Steady  State  Condition 

As  it  interacts  with  the  solar  wind,  the  magnetospheric/tail  system  appears  to  be  trying  to 
establish  an  equilibrium  when  it  periodically  releases  some  of  its  excess  energy  during  magnetic 
storms.  It  seems  possible  that  die  magnetospheric/tail  system  may  be  trying  to  reach  an 
equilibrium  state  having  the  lowest  total  energy.  This  is  the  assumption  of  this  thesis  and  will 
be  discussed  in  detail  later. 

In  order  to  test  this  assumption,  one  needs  to  have  measurements  of  the  magnetotail  while  it 
is  in  a  steady-state  equilibrium  condition.  Unfortunately,  due  to  the  constant  influx  of  energy 
from  the  solar  wind,  this  condition  rarely  exists,  but  it  may  exist  for  some  special  cases.  One 
case,  as  mentioned  above,  may  exist  for  a  brief  period  of  time  immediately  after  a  magnetic 
storm  when  the  magnetotail  has  released  its  excess  energy.  Another  case  may  occur  when  the 


19 


earth’s  magnetic  field  is  exposed  to  a  large  northerly  component  of  the  IMF  for  an  extended 
period  of  time.  When  die  IMF  is  pointing  northward,  die  solar  wind  has  minimal  interactions 
with  the  earth's  magnetic  field.  By  looking  at  the  specific  parameters  within  the  magnetotail 
system,  one  may  be  able  to  determine  whether  the  system  is  in  a  minimum  energy  equilibrium 
state  or  not  This  diesis  explores  this  possibility. 

Chapter  II  of  this  thesis  explains  a  technique  for  finding  die  total  minimum  energy  of  a  system 
not  in  thermodynamic  equilibrium.  This  chapter  discusses  how  the  "Calculus  of  Variations"  is 
used  to  determine  which  systems  of  equilibrium  states  are  also  minimum  energy  states. 

Chapter  III  describes  an  idealized  magnetotail  model  used  in  this  research  and  discusses  how 
die  energy  minimizing  technique  described  in  Chapter  11  is  applied  to  it.  Also,  this  chapter  will 
explain  what  assumptions  were  used  to  create  die  idealized  model  and  why. 

Chapter  IV  addresses  the  computer  program  used  to  solve  a  series  of  differential  equations 
which  evolved  from  the  calculus  of  variations  technique.  It  also  describes  the  model  outputs  and 
any  associated  graphs  and  plots. 

Chapter  V  compares  the  model’s  results  with  real  data  from  the  actual  earth’s  magnetotail.  If 
the  results  approximate  typical  conditions  observed  in  the  magnetotail,  then  we  will  be  able  to 
say  that  we  will  have  a  working  model  for  describing  the  earth’s  magnetotail  energy  state  while 
it  is  in  equilibrium.  If,  on  the  other  hand,  die  model  does  not  adequately  simulate  die 
magnetotaQ’s  energy  state,  then  we  will  have  eliminated  one  possible  method  for  describing  the 
earth’s  magnetotail  energy  structure. 

Chapter  VI  draws  conclusions  based  upon  die  results  of  this  research  and  points  out  possible 


areas  for  further  research. 
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Chapter  II 

MINIMUM  TOTAL  ENERGY 
Systems  in  Thermodynamic  Equilibrium 

If  a  system  is  completely  isolated  from  its  surroundings,  it  is  in  thermodynamic  equilibrium. 
Determining  stable  equilibria  of  systems  in  these  states  is  done  by  minimizing  the  appropriate 
potential  energy  function  for  that  thermodynamic  equilibrium.  Several  potential  energy 
functions  and  their  specific  uses  are  listed  below  (Adkins,  1968). 

TABLE  I.  Potential  Energy  Functions. 

Function  Description 

Enthalpy  •  can  be  used  for  determining  die  maximum 

amount  of  heat  content  available  within  a 
given  system. 

Helmholtz  function  -  can  be  used  for  determining  the  maximum 
amount  of  mechanical  work  which  can  be 
extracted  from  a  given  system. 

Gibbs  function  -  can  be  used  for  determining  the  maximum 
amount  of  free  energy  available  within  a 
given  system. 

Therefore,  in  thermodynamic  equilibria,  minimizing  a  particular  potential  energy  function 
yields  a  stable  equilibrium  state. 


21 


Systems  Not  in  Thermodynamic  Equilibrium 

But  what  if  die  system  does  not  remain  in  thermodynamic  equilibrium.  Such  is  the  case 
where  radiation  is  present  (Pippard,  1957).  We  will  assume  that  these  equilibria  have  minimum 
total  energy  rather  than  minimum  thermodynamic  free  energy  (or  potential  energy).  For  such 
systems,  we  will  attempt  to  find  their  stable  equilibria  by  using  calculus  of  variations  to 
determine  their  minimum  total  energy  states. 

The  first  step  prior  to  using  this  technique  is  to  find  an  equation  which  adequately  describes 
the  energy  of  the  system  in  question.  The  total  energy,  U,  of  our  system,  can  be  represented  by 
adding  up  all  of  the  applicable  energy  contributions,  assuming  that  there  are  no  dissipation  and 
advection  terms: 
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(3) 


The  next  step  is  to  identify  the  constraining  equations  which  apply  to  the  problem.  These 
include  Maxwell’s  electrodynamic  equations  and  equations  of  state  for  the  plasmas  in  question. 

The  final  step  of  this  technique  is  to  determine  which  functions  appearing  in  the  energy 
equations,  Eq.  (3),  within  die  limits  imposed  by  the  constraint  equations,  yield  a  minimum  total 
energy.  These  functions  include  variables  such  as  magnetic  fields,  electric  fields,  particle 
densities,  etc.,  which  are  dependent  upon  the  independent  space  variables.  One  could  randomly 
insert  different  functions  for  each  of  die  dependent  variables  in  Eq.  (3),  subject  to  die 
appropriate  constraints,  until  the  lowest  possible  value  for  U  was  reached.  Unfortunately,  this 
approach  is  very  tedious  and  time  consuming  and  would  not  guarantee,  with  any  certainty,  that 
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another  combination  of  values  may  not  generate  an  even  lower  value  for  U.  It  is  due  to  this 
uncertainty  that  another  technique  must  be  used  to  find  the  lowest  possible  value  for  U,  or 
minimum  total  energy.  The  technique  uses  the  calculus  of  variations;  it  works  as  follows. 

Calculus  of  Variations 

Suppose,  for  example,  we  want  to  find  fire  minimum  total  energy  of  a  given  system.  Since  the 
minimum  total  energy  of  a  system  is  a  special  case  of  die  system’s  total  energy,  we  will  look  at 
the  total  energy  first  without  restricting  it  to  being  a  minimum.  To  illustrate  this  method,  let  us 
assume  the  total  energy  of  die  system  can  be  represented  as  a  function  of  two  dependent 
variables,  A  and  B  *=  dh/dz,  which  are  functions  of  one  independent  variable  z.  With  this  in 
mind,  we  can  express  die  total  energy  of  a  given  system  as. 


TOTAL 

ENERGY 


(U) 


[y  f  (  A  (  z  ).  B  (  z  ),  z  )  dz  . 


(4) 


In  order  to  find  the  minimum  total  energy  of  die  system,  we  must  determine  what  functions  A 
and  B  in  Eq.  (4)  must  be  in  order  for  U  to  achieve  the  lowest  value  possible.  The  trick  is  to 
discover,  with  some  degree  of  certainty,  what  these  functions  are.  Remember,  die  values  of  A 
and  B  may  be  constrained  by  one  or  more  external  equations  as  well  as  boundary  conditions  that 
may  be  imposed  on  the  system.  Therefore,  in  order  to  find  A(z)  and  B(z),  we  must  use  the 
calculus  of  variations. 

To  use  this  technique  we  first  assume  that  we  know  what  die  ’correct’  values  of  A  and  B  are, 
and  express  diem  as  Arniq  and  Bmjn.  Substituting  these  new  variables  into  Eq.  (4)  yields 
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MINIMUM 

TOTAL  (U)  -  f  ((Amki(z),BmlJz),z)<lz.  (5) 

ENERGY  w 


But  how  do  we  determine  what  die  values  of  A^,,  and  Btnin  are?  To  find  out,  we  select  some 
arbitrary  functions,  a(z)  and  b(z),  multiply  them  by  an  infinitesimal  dimensionless  parameter  a, 
and  add  them  to  A^;„  and  Bn,;n  to  generate  two  entirely  new  functions.  A’  and  B\ 


A'  (  z,  «  )  -  A z,  0  )  ♦  am  (  z  )  ; 
B;  (  z,  a  )  *  B "Ji  z,0)  +  ab(z). 


One  of  the  condidons  of  die  calculus  of  variations  is  that  all  values  of  the  ’correct1  functions, 
Amtw  and  Bmln,  must  be  the  same  as  those  of  the  new  functions,  A  and  S’,  at  the  systems 
boundaries  or  endpoints.  In  this  case,  the  endpoints  are  at  the  top,  z,,  and  bottom,  Zj,  of  the 
system. 

Rewriting  Eq.  (4)  using  the  new  functions  A’(z/x)  and  ff  (z/0  gives  us  the  following 
expression. 


-  /yf(A/(z,«X  B'(z,«),  z)dz. 


(7) 


By  making  aa  and  ab  equal  zero  in  Eq.  (6),  we  see  that  A  and  S’  in  Eq.  GO  will  have  the  same 
values  as  A^  and  B,,^  in  Eq.  (5)  and  the  system  will  be  in  a  minimum  total  energy  state. 

We  begin  the  process  of  minimizing  by  looking  for  cases  where  die  derivative  of  U  with 


respect  to  a  equals  zero.  Assuming  dU/da  is  a  continuous  function,  the  system  will  be  in  an 
equilibrium  energy  state  when  dU/d«  equals  zero. 


da 


it  /  \  f  /  *  3A'  3f  3B'  . 

U(a)  "  fr(u'*r  +  W'*r)dz 


a  o 


(8) 


Unfortunately,  it  is  impossible  to  know  for  certain  whether  a  particular  equilibrium  coincides 
with  die  system’s  minimum  total  energy  state  without  doing  any  further  tests.  The  reason  for 
this  is  that  a  system  may  have  many  different  equilibrium  states  (Fig.  21). 

Another  way  to  insure  that  an  equilibrium  state  coincides  with  a  system’s  minimum  total 
energy  is  to  take  the  limit  of  Eq.  (7)  as  «  goes  to  zero. 


u«i. 


lim 

o-O 


U(«  ) 


0%fy  f<A',  B'.  z)dz; 

/y  f  (  A*.*  « a,  «b,  z  )  dz  . 


(9) 


We  can  see  as  a  goes  to  zero,  Eq.  (9)  begins  to  look  like  die  minimum  total  energy  equation 
(Eq.  (5)).  Therefore,  in  order  to  find  die  equilibrium  state  of  a  system  which  also  coincides  with 
its  minimum  total  energy  (U),  we  need  to  apply  both  the  partial  differentiation  and  limit  to  the 
energy  equation  at  die  same  time  as  shown  in  Eq.  (10)  and  Fig.  22. 


STABLE  BQUHJBRIUM 
POTENTIAL  ENERGY  -  f 


UNSTABLE  EQUILIBRIUM 
POTENTIAL  ENERGY  -  ? 


FIG.  21.  Diagram  showing  that  a  system  can  have  many  equilibrium  states  (stable  and 
unstable). 


ENERGY 


FIG.  22.  Since  die  value  for  U(  a  >  0 )  will  always  be  greater  than  U(  a  «  0  ),  we  can  be 
assured  that  we  will  find  the  equilibrium  condition  with  the  least  possible  energy  when  a  goes 
to  zero. 
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Equation  (10)  must  hold  independent  of  the  choice  of  functions  a  and  b  in  Eq.  (9)  because  as 
a  goes  to  zero,  A*  and  S’  take  on  the  values  of  and  Bj^.  By  taking  the  limit  as  «  goes  to 
zero,  however,  we  eliminate  the  values  of  the  arbitrary  functions,  a  and  b,  so  how  can  we  vary 
them  so  that  we  are  left  only  with  A*  *=  A^^  and  —  BTni-n?  The  answer  is  to  use  integration 
by  parts  (Goldstein,  1980). 

Integration  by  parts  allows  one  to  factor  the  arbitrary  functions,  a  and  b,  out  of  Eq.  (6)  and 
leave  only  a  differential  equation  within  the  integral  based  on  the  values  A’  and  B\  Because 
ff  =  dA/dz,  the  second  integral  of  Eq.  (8)  is 


f  ( JL  ®:)dz  =  r>(JLE±L)dz; 

Jv  aB'  da  3B'  *a« 


ar  a2  a  7 


now  we  integrate  by  parts 


p  ( JL  )  dz  = 

■»*,  ae'  drift  3B'  da  *>  Jh  dz  dB'  da 


(11) 


Since  the  correct  and  arbitrary  functions  must  be  equal  at  the  system’s  endpoints/boundaries,  z  j 
and  Z2«  tire  partial  derivative  of  A  with  respect  to  a  in  the  first  part  of  Eq.  (11)  will  disappear, 
leaving  us  with 


df  3B 
3B'  da 


)  dz 


(  df  x  dAf 
J*,  dz  1  m'  *  da 


dz  . 


(12) 
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Therefore, 


£”(*> 


</v 


aA' 


/  «JL -AJL )(*L'))dz 

'v  aA'  dz  ae'  Sa 


or. 


^-U<«)  -  fh((P(*)X^))dz, 
da  J*t  da 


(13) 


where. 


F<*) 


JL  -  A  JL 

dA'  a*  a}'  " 


(14) 


We  now  have  an  expression  F(z)  that  is  written  in  terms  of  A*  and  IF  only.  Since  die  partial 
derivative  of  A  with  respect  to  «  can  be  any  arbitrary  value,  F(z)  must  be  equal  to  zero  in  order 
to  satisfy  the  condition  we  made  in  Eq.  (8).  Therefore,  by  solving  the  differential  equation  F(z) 

-  0  for  A’  and  IF,  we  are  able  to  calculate  die  minimum  total  energy  of  die  system.  F(z)  -  0  in 
Eq.  (14)  is  called  the  Euler-Lagrange  equation.  This  thesis  will  apply  this  method  to  an  idealized 
model  of  the  earth’s  magnetotail  and  compare  die  results  with  real  data  of  the  earth's  actual 
magnetotail. 
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Chapter  III 

MAGNETOTAIL  MODEL 

Model  Description 

The  model  we  use  in  this  thesis  to  calculate  die  minimum  total  energy  of  the  plasma  sheet 
within  the  earth’s  magnetotail  is  depicted  as  a  rectangular  slab  (Fig.  23)  of  infinite  length  and 
width  and  a  finite  height,  h.  The  axes  of  this  model  magnetotail  are  defined  by  a  right-handed 
Cartesian  system  where  die  positive  x-axis,  which  represents  the  length  of  the  slab,  points 
towards  the  sun.  The  y-axis,  which  represents  the  width  of  the  slab,  is  directed  perpendicular  to 
the  x-axis  in  die  duskward  direction.  Finally,  die  z-axis,  which  represents  the  height  of  die  slab, 
points  northward,  perpendicular  to  both  the  x  and  y  axis. 

By  using  some  key  assumptions  and  applying  an  appropriate  set  of  equilibrium  equations  to 
this  magnetotail  model,  we  should  be  able  to  calculate  its  minimum  total  energy.  Comparing 
die  results  of  the  model’s  calculations  to  certain  parameters  observed  in  die  plasma  sheet  of 


FIG.  23.  Drawing  showing  die  layout  of  the  magnetotail  model. 
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the  earth’s  magnetotail  during  quiet  (steady  state)  conditions  will  help  determine  die  validity  of 
this  model. 

As  mentioned  above,  several  assumptions  are  used  in  the  model  to  keep  it  within  the  scope  of 
this  thesis.  The  assumptions  used  in  die  model  are  as  follows: 

a)  System  is  in  a  steady-state  condition  (d/dt  -  0). 

b)  System  is  in  equilibrium. 

c)  Particle  number  densities,  ions  (n^),  and  electrons  (no)  vary  in  the  z  direction  only. 
In  other  words,  die  particle  number  densities  remain  constant  within  die  x-y  planes. 

d)  Ion  velocities  are  restricted  to  die  x-y  plane  to  keep  previous  assumption  true. 

e)  Only  one  species  of  particles  (ions  or  electrons)  moves,  the  other  remains  stationary, 
only  ions  possess  velocity.  This  assumption  is  used  to  simplify  calculations  within  this  diesis. 

0  The  x  and  y  components  of  the  electric  field  equal  zero  (E^  =  Ey  =  0).  Only  the  z 
component  of  the  electric  field  has  a  value  which  varies  only  in  the  z  direction  (E^z)  «*  value). 
Without  this  assumption,  the  ions,  within  the  model,  would  experience  a  continuous  acceleration 
within  the  infinite  x-y  planes  due  to  a  changing  electric  field.  This  condition  would  eventually 
cause  the  ions  to  reach  an  infinite  velocity,  an  unrealistic  situation. 

g)  The  y  and  z  components  of  die  magnetic  field  equal  zero  (By  =  Bz  *=  0).  Only  the  x 
component  of  the  magnetic  field  has  any  value  and  this  value  varies  only  in  die  z  direction 
(Bx(z)  *=  value).  Tht  y  and  z  components  of  the  magnetic  field  are  assumed  to  be  zero  because 
the  model  is  approximating  the  center  of  the  earth’s  magnetotail  near  the  midnight  sector.  In 
this  region  of  the  earth’s  magnetotail,  the  magnetic  field  lines  are  very  nearly  parallel  to  the  x 
axis  with  negligible  y  and  z  components. 

h)  Polarization  and  magnetization  of  the  plasma  are  assumed  to  be  negligible  (M  =  0, 

P  =  0). 

i)  The  plasma  within  the  magnetotaQ  model  is  fully  ionized. 
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By  using  these  simplifying  assumptions,  we  are  left  with  a  very  simple  or  zero  order  model 
which  may  not  include  all  the  necessary  details  to  adequately  describe  the  process  observed 
within  the  central  plasma  sheet  of  the  earth’s  magnetotaO. 

Magnetotail  Model's  Total  Energy  Equation 

As  stated  in  Chapter  II,  we  must  develop  a  total  energy  equation  to  describe  the  model  system 
above  which  is  not  in  thermodynamic  equilibrium.  To  do  this,  we  must  determine  what  energy 
terms  we  want  the  total  energy  equation  to  have.  Since  we  are  attempting  to  model  a  region 
where  electric  and  magnetic  fields  are  prevalent,  it  seems  logical  that  the  total  energy  equation 
should  contain  terms  describing  the  energy  changes  associated  with  these  fields.  We  also  need  a 
kinetic  energy  term  and  an  energy  term  that  account  for  moving  plasma  particles  (ions  or 
electrons)  and  kinetic  pressures  found  within  the  earth’s  magnetotail.  Therefore,  we  want  to 
develop  a  total  energy  equation  for  the  model  system  that  has  the  following  form: 

E  Field  B  Held  Kinetic  Energies 
Energy  Energy  Energy  Associated 

Per  ♦  Per  ♦  Per  with  )  dx  . 

Unit  Unit  Unit  Kinetic 

Vol  Vol  Vol  Pressure 

Now  we  need  an  expression  for  each  term  in  die  total  energy  equation  shown  above.  We  can 
get  expressions  for  die  electric  and  magnetic  field  energies  by  combining  Poyn ting’s  theorem 
with  Maxwell’s  equations  which  results  in  the  expression  shown  in  Eq.  (15). 
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Energy 

Change 

Per 

Unit  Time 

dt 


Electric/ 
Magnetic 
Field  Strength 
Change  Pa 
Unit  Time 


Power 

+  Per  ♦ 

Unit  Vol 


Energy 

Advecdoa 

Per 

Unit  Vol 


»  f  -Si)*./  <J.8>*./,<^t>a. 

at  Jv  v  2  2|10  J*  1 *  Po 


(IS) 


where, 

Total  Energy  of  a  System; 
Electric  Field; 

Magnetic  Field; 

Electric  Current, 
Permittivity  of  Free  Space; 
Permeability  of  Free  Space; 
Unit  Volume; 

Unit  Area. 


It  is  standard  practice  to  define  the  «qE2/2  and  B 2/2iiq  terms  in  Eq.  (15)  as  electric  and 
magnetic  field  energy  densities,  respectively.  These  can  be  used  in  the  total  energy  equation. 
We  can  also  rewrite  the  energy  advection  term  in  Eq.  (15)  in  terms  of  kinetic  energy. 
Substituting  these  expressions  into  the  total  energy  equation  yields: 


Toad 

Energy 


U 


B  Field 

B  Field 

Kinetic 

Energies 

Energy 

Energy 

Energy 

Associated 

)dx 

■  /,< 

Per 

♦  Per 

♦  Per 

♦  with 

Unit 

Unit 

Unit 

Kinetic 

Vol 

Vol 

Vol 

Pressure 

■  /,< 

*oE* 

♦  *1 

5  nmqv* 

♦  OnT 

)* 

2 

2p0 

2 

(16) 
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where, 


n  —  Particle  Number  Density; 
nx  -  Ions; 
na  *  Electrons; 
m  =  Particle  Mass; 
q  =  Particle  Qon/Electron)  Velocity; 
C  =  Adiabatic  Constant; 

Y  -  Ratio  of  Specific  Heats. 


The  last  term  of  the  total  energy  equation,  CnT,  which  accounts  for  die  kinetic  pressure  found 
within  the  earth’s  magnetotail  system,  is  derived  from  the  Perfect  Gas  Law. 


PV  -  (-5-XRHT)  -  (n)(k)(T); 
na 

or 

cf  )(^)(T>; 

(»)  (k)  (T); 

where. 


P  - 

P  = 


p  =  Pressure; 

V  =  Volume; 

r  =  Universal  Gas  Constant, 

Nx  “  Avogadro’s  Number; 

X  *=  Temperature  (in  degrees  Kelvin); 
n  =  Number  Density; 

It  -  Boltzmann’s  Constant 


Assuming  we  have  adiabatic  conditions,  we  can  say  that  the  internal  energy  of  a  system  is  equal 
to  the  energy  extracted  from  the  system  by  die  performance  of  work  (First  Law  of 
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Thermodynamics) . 


Au  -  -  w  . 


However,  the  internal  energy  change  for  any  process  for  an  ideal  gas  can  be  rewritten  as. 


Au  "  n  Cy  AT  , 


where  cv  is  the  heat  capacity  of  a  system  at  a  constant  volume.  In  addition,  the  variable 
representing  the  performance  of  work  can  be  rewritten  as. 


w  -  PAV, 


where  we  assume  that  die  pressure  P  remains  approximately  constant,  for  a  very  small  change  in 
volume  8V.  Substituting  these  new  expressions  into  the  thermodynamic  equation  yields: 


n  Cy  AT  -  -  P  AV  . 


Now  we  want  to  write  die  above  equation  in  terms  of  temperature  and  volume  only.  To  do  this 
we  begin  with  the  following  form  of  die  perfect  gas  law. 


nRT 

V 
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and  substitute  it  into  the  rewritten  thermodynamic  equation  yielding: 


n  Cy  AT  -  -(“JLI)AV; 


AT 

T 


(A  )  AV  =  _(5L^V)  AV  .  .(T.1}  AY 
Cy  V  Cy  V  V 


Integrating  the  above  expression  yields: 


(  T  )  (  V  )  T'1 

n  R 

(PV)(V)H 

(P)(V)» 


Constant ; 

Constant  ; 

Constant  (  n  R  )  ; 
Constant ; 


but. 


v  -  (  R  ,  (rn)(T)  . 

'  *  /  n  \  * 


N, 


(P) 


(k) 

(a) 


dH-Ar  ) 

_ n  k 

(P) 


-i  . 

t 
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therefore, 


P<  V)* 
P((n)‘)T 
P(n)-t 

P 


Constant  : 


Constant  ; 


^ - *  —  j  . 

LOQSltZU  ) 


Constant  (  n  )  T 


Cn  »  . 


For  computing  purposes,  it  is  easier  to  deal  with  dimensionless  equations  than  those  with 
dimensions.  In  order  to  make  Eq.  (16)  dimensionless,  we  redefine  the  dependent  variables  in 
such  a  way  that  each  energy  term  becomes  dimensionless.  This  has  die  same  effect  as  setting 


m  -  1  ; 

q  -  i ; 
«0  ■  1  » 
1*0  ■  1  • 


which  gives  us  die  following  equation. 


U 


♦  ♦  Cn*  )  dt  . 

2 


(17) 


It  is  important  to  note  here  that  the  dependent  variables  in  Eq.  (17)  are  dimensionless  and  are 
not  die  same  as  die  dimensioned  variables  shown  in  Eq.  (16).  It  will  be  shown  later  in  this 
diesis  that  by  multiplying  the  dimensionless  variables  in  Eq.  (17)  by  an  appropriate  scaling 
quantity,  we  can  end  up  with  the  dimensioned  variables  shown  in  Eq.  (16).  Therefore,  unless 
specified,  die  dependent  variables  presented  in  this  paper  are  dimensionless. 


In  order  to  account  for  the  two  different  particle  species,  ions  (n^),  and  electrons  G^),  we 
must  apply  some  of  the  assumptions  discussed  earlier  in  this  chapter  to  Eq.  (17).  They  are  as 
follows: 
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Ex 

By 

v, 

n2vx 


Ey  «  o,  (E  Fid  only  has  a  z  component); 
Bs  -  o ,  (B  Fid  only  has  an  z  component); 
0 ,  (No  vertical  velocities); 

0 .  (Electron  velocities  are  zero). 


Applying  these  assumptions  to  Eq.  (17)  yields: 


U 


♦  Cfr1  )  <h  . 


(18) 


Rewriting  Eq.  (18)  in  an  expanded  form  gives  us  Eq.  (19).  This  new  equation  represents  the 
total  energy,  U,  of  the  model  magnetotail  system  described  earlier  in  this  chapter.  For  die 
remainder  of  this  thesis,  Eq.  (19)  will  be  referred  to  as  the  total  energy  equation. 


U 


nr  /:<¥♦¥ 


BjV, 


®tvy 


+  C,n,T  ♦  CjtijT  )  dxdydi  .  (19) 


In  order  to  come  up  with  the  lowest  possible  value  for  the  total  energy,  U,  of  die  model 
system,  we  must  first  find  the  minimum  values  of  each  of  die  terms  in  Eq.  (19).  Of  course  one 
way  to  do  this  would  be  to  arbitrarily  set  each  term  in  Eq.  (19)  to  zero  and  make  die  total 
energy,  U,  of  die  system  equal  to  zero,  but  this  would  depict  a  very  unrealistic  situation  in  the 
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real  world. 

Since  we  are  looking  for  a  minimum  total  energy  of  a  system  in  equilibrium,  the  variables 
within  Eq.  (19)  must  be  governed  by  one  or  more  equilibrium  equations.  These  equilibrium 
equations,  in  turn,  will  be  constrained  to  some  degree  by  conditions  outside  the  system. 
Therefore,  a  proper  set  of  boundary  conditions  must  be  defined  along  with  the  appropriate  set  of 
equilibrium  equations.  By  having  die  appropriate  set  of  equilibrium  equations  and  boundary 
conditions,  we  will  be  able  to  define  a  range  of  values  (Appendix  Q  for  each  term  in  Eq.  (19). 

Equilibrium  Solutions 


Maxwell's  Equilibrium  Equations 

The  equilibrium  equations  derived  from  Maxwell’s  equations  are  used  to  govem/restrict  the 
electric  and  magnetic  fields  defined  within  the  model  system.  The  first  of  these  equilibrium 
equations  is  Gauss’  Law  of  Electric  Charge  and  is  expressed  as 


V  •  £  -  ^  ^  +  3i  -  .  (20) 

dx  dy  dz  c0 


We  know  the  total  charge  density  in  Eq.  (20)  can  be  written  as 


P**i 


m 

vol 


(+q)(Pi)*(-q)(na)  _  n,  - 
vol  vol 


(21) 


Therefore,  by  substituting  Eq.  (21)  into  Eq.  (20)  and  applying  the  following  assumptions. 
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we  end  up  with  the  following  expression: 


3E, 

dz 


1*1  +n2- 


(22) 


Eq.  (22)  provides  us  with  the  first  equilibrium  equation  derived  from  Maxwell’s  equations  which 
governs  the  model  magnetotail’s  electric  field  and  particle  number  densities. 

The  second  equilibrium  equation  derived  from  Maxwell’s  equations  is  Faraday’s  Law  of 
Induction  and  is  expressed  as 


V  x  g 


ae  _  , 

<-=?  “-j?)* 
dy  dz 

cE  3E 
♦  (  — - i  )  * 

dz  dz  3 


(23) 


Using  the  assumption  that  the  z  component  of  the  electric  field  only  varies  in  the  z  direction 
(Ej(z))  and  that  the  model  system  is  in  a  steady  state,  we  can  see  that  Eq.  (23)  reduces  to 
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Vxfi  -  0. 


(24) 


The  third  equilibrium  equation  derived  from  Maxwell’s  equations  is  Gauss’  Law  of  Magnetism 
and  is  expressed  as 


V  •  $ 


-  0, 


(25) 


which  is  covered  by  the  assumptions  previously  mentioned.  This  equilibrium  equation  states 
that  magnetic  field  lines  do  not  start  or  stop  in  space. 

The  final  equilibrium  equation  derived  from  Maxwell’s  equations  is  Ampere’s  Law  and  is 
expressed  as 


Vxfi 


3B  50  3B  B  8B 

(_* - Z  )  *  +  (  — 5 - ~  )9  +  (  — l 

v  dy  &  az  dx  3  ax 

p0  ( i  *  v  x  a  *  ^  *  £ ) . 


(26) 


Using  fire  assumptions  that  the  y  and  z  components  of  file  magnetic  field  equal  zero  (By  =  Bz  = 
0),  the  x  component  of  the  magnetic  field  varies  in  the  z  direction  only  (Bx(z)),  file  polarization 
and  magnetization  aspects  of  file  electric  and  magnetic  currents  equal  zero  (M  =  P  =  0),  the 
permittivity  of  free  space  is  dimensionless  (|i0  =  1),  and  file  system  is  in  steady  state 
(d/dt  =  0),  we  can  rewrite  Eq.  (26)  as 
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3BX  ~ 

><ir)  * 


(  P  )  (  9  )  -  (  ®i  “  °a  )  (  Vji  ♦  v^  ♦  ) 


(27) 


At  this  point,  we  use  one  of  die  model's  assumptions  that  electrons  (n2)  have  no  velocity. 
Normally  when  considering  velocities  due  to  fi  x  3  drifts,  both  the  ions  and  electrons  within  a 
plasma  drift  However,  results  from  other  works  attempting  to  model  the  Venus  flux  ropes 
indicate  that  models,  such  as  the  one  presented  in  this  paper,  may  woric  better  by  keeping  one  of 
die  species  of  particles  stationary  (private  conversation  with  Dr.  Edwards).  With  this  in  mind, 
combining  like  unit  vector  terms  in  Eq.  (27)  yields 


(28) 


This  expression  provides  us  with  the  second  useful  equilibrium  equation  derived  from  Maxwell’s 
equations.  This  equation  governs  the  model’s  magnetic  field  and  ion  velocities  (or  kinetic 
energies). 

In  summary,  the  equilibrium  equations  derived  from  Maxwell’s  equations  are  shown  below. 


V-  g 

V  x  g 
V-  g 


"  “i  - 


as, 

dt 

0; 

0; 

»x 

-3T  = 


(29) 


Vx  g 
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Kmetic/Force  Equilibrium  Equations 

The  equilibrium  equations  derived  from  the  force  equation  are  used  to  govern/restrict  the 
kinetic  energies  and  number  densities  of  the  particles  contained  within  the  model  magnetotail. 
The  force  equation  is  expressed  as  follows: 


9  >  mi. 


(30) 


However,  since  we  are  dealing  with  a  specific  volume,  we  need  to  express  the  force  in  Eq.  (30) 
as  a  force  per  volume  (Eq.  (31)). 


?  _  m 

vol  vol 


(31) 


Rewriting  Eq.  (31)  so  that  ?  represents  the  force  per  volume  and  pm  represents  the  mass  per 
volume,  we  get 


?  -  pm  i  . 


(32) 


However,  we  also  know 


?  -  Lorentz  Force  +  Pleasure  Gradient  Force 
i  =  p.(fit?x§)  -  V?; 


(33) 
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and. 


P. 


m 

vcd 


(+q)(»i)*(-0(ni> 

vol 


Total  number 
N  -  of  particles 
per  unit  vol  , 


(34) 


and  the  acceleration  is  given  by  the  convective  derivative, 


i 


♦  ( 9 • V) V  )  . 


(35) 


Applying  Eqs.  (33,  34,  and  35)  to  Eq.  (32)  yields  die  following  equation(s). 


?  -  (P»)d); 

(pa)(fi  +  9xfi)-(tf)-(N  )(?  *(9-V)9). 

at 


(36) 


But  according  to  Eq.  (34),  em-N,  therefore  Eq.  (36)  can  be  written  as 


(N)(fi  +  9xfi)-(9P)  -  <  N  )  (  —  ♦  (  9  •  V  )  9  )  .  (37) 

dt 


Applying  the  Perfect  Gas  Law  (P  =  nkT  =  CnT)  to  the  pressure  gradient  force  in  Eq.  (37)  yields 
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(N)(fi  +  Vxfi)-(CY(N)T'1V(N)) 

-(N)(-f?*(¥*V)¥).  (38) 

ot 


Dividing  both  sides  by  the  total  number  density  (N)  yields 


(  g  +  V  x  g  )  -  (  Cy  (  N  )  T*2  V  (  N  )  )  -  (££♦(?•  V  )  S  )  .  (39) 

at 


Using  the  assumption  that  the  electric  field  has  only  one  component  z,  which  varies  only  in 
the  z  direction,  we  can  say 


g  -  (  E,  )  i  .  (4°) 

Since  we  have  also  assumed  the  y  and  z  components  of  the  magnetic  field  equal  zero  (By  =  Bz 
*  0)  and  the  x  component  varies  only  in  the  z  direction  (Bx(z)),  we  get  the  following 
equation(s): 

9  x  g  =  (  VjB,  -  v,By  )  i 

+  (  v,By  -  VjB,  )  t  ; 

-  (  )  i 


(41) 
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Since  die  pressure  gradient  force  only  varies  in  the  z  direction  in  this  model,  we  get 


(  Cy  (  N  )  T*2  V  (  N  )  )  =  (  Cy  (  N  )  T"2  fLOLl  )  ±  .  (42) 

CL 


Looking  at  die  right-hand  side  of  Eq.  (42)  only,  we  can  derive  the  following  expression 


V)  9  )  -  (0 ♦ ( * • V) V)  -  0. 


(43) 


Combining  Eqs.  (39,  40,  41,  and  42)  into  Eq.  (38)  yields 


(  E,  )  £  -  (  v  B  )  £  -  (  Cy  (  N  )  T  '2  — (  N  >  )  £  -  0  .  (44) 

'  dz 


Rewriting  Eq.  (44)  so  that  the  change  of  number  density  with  height  is  on  the  left  side  of  the 
equation,  we  are  left  with  the  following  scalar  expression: 


3(N) 

dz 


(  J_  (N)s-t)(Ei  _  VyB  ) 

Cy 


(45) 


Since  the  values  for  the  ion  or  electron  number  density  are  (transforming  to  dimensionless 


variables) 
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Nta.  -<♦*>(■,)  -(*l>(«i); 

Na«*«  -  <-q)(»j)  =  (  _1  )  (  Hj  )  J 


(46) 


they  can  be  substituted  into  Eq.  (45),  giving  vis  an  expression  for  the  change  of  the  ion  number 
density  across  the  magnetotail  from  the  southern  to  the  northern  boundaries. 


3(fli) 

dz 


*  <  ib  < ",  >: > < %  -  >  • 


(47) 


Note  that  we  have  assumed  the  electron  velocities  are  zero.  Using  the  same  steps  described 
above,  we  can  derive  an  expression  for  the  change  of  the  electron  number  density  across  the  tail. 


dz 


(48) 


Therefore,  we  have  a  set  of  equilibrium  equations  (Eqs.  (47)  and  (48)),  derived  from  die  force 
equation,  which  governs  die  change  of  ion  and  election  number  densities  across  die  tail  for  the 
model  magnetotail.  Summarizing,  the  force  equilibrium  equations  are 

£J_=l2  .  .(lm*i)(v.,»,); 

1 it2  ■ 


(49) 
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At  this  point,  we  have  a  set  of  equilibrium  equations  which  constrain  the  values  of  die  terms 
within  the  total  energy  equation  (Eq.  (19)).  We  must  now  specify  die  boundary  conditions. 
Additionally,  we  must  add  another  condition  which  will  prevent  the  particles  (ions  and 
electrons)  contained  within  the  model  system  from  arbitrarily  changing  species.  These 
conditions  are  discussed  below. 


Boundary  Conditions 

Equilibrium  equations  alone  will  not  restrict  the  value  range  for  each  of  the  terms  in  the  total 
energy  equation  (Eq.  (19)).  Their  values  face  restrictions  on  the  boundaries  of  the  model 
system.  Since  the  model  system  only  has  boundaries  on  the  top  and  bottom,  boundary 
conditions  will  only  be  applicable  at  these  locations. 

Number  Density  Boundary  Conditions 

For  this  model,  it  is  assumed  that  charged  particle  density  (ion  or  electron)  outside  the  system 
is  negligible  compared  to  the  number  densities  inside.  Therefore,  we  can  see  that  the  particle 
number  density  inside  the  model  must  approach  zero  as  one  gets  closer  to  the  top  or  bottom 
boundaries  of  the  model  system.  We  must  have  a  set  of  boundary  conditions  which  depict  this 
situation.  Because  boundary  conditions  on  nj  and  are  not  specified,  we  must  establish 
boundary  conditions  on  the  corresponding  Lagrangian  multipliers,  which  will  be  introduced  and 
discussed  later.  These  particular  conditions  are: 


MO)  -  o; 
MO)  =  0; 


x,  (h)  -  0; 

Mh)  «  0. 


(50) 
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Electric  Field  Boundary  Conditions 

Unlike  the  charged  particle  number  densities,  the  model  assumes  there  is  no  electric  field 
outside  the  model  system.  This  means  that  the  electric  field  must  equal  zero  at  the  boundaries 
of  the  system,  which  happens  to  be  the  top  and  bottom  planes  of  the  model  system.  The 
boundary  conditions  which  restrict  fire  equilibrium  equations  governing  the  model’s  electric  field 
are  depicted  below. 

E,(0)  -  0;  Et(h)  -  0.  (51) 

Magnetic  Field  Boundary  Conditions 

Unlike  the  charged  particle  number  densities  and  the  electric  field,  the  model  does  assume  a 
substantial  magnetic  field  does  exist  outside  the  model  system.  Therefore,  the  magnetic  field 
strength  just  inside  the  model's  boundaries  must  equal  whatever  magnetic  field  strength  lies  just 
outside  the  boundary.  The  boundary  conditions  which  depict  this  situation  within  the  model  are 
as  follows: 

8,(0)  -  0;  B,(h)  -  +8^.  (52) 

Particle  Conservation  Equations 

These  equations  are  used  to  ensure  the  number  density  of  the  individual  species  (ions  and 
electrons)  contained  within  the  model  magnetotail  remains  conserved.  So  far,  we  have  had 
Maxwell’s  equations,  which  only  guarantee  that  die  sum  total  of  the  charged  particles  within  the 
model  system  will  remain  conserved  (Eq.  (22)).  In  other  words,  a  charged  particle  (proton  or 
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electron)  will  not  suddenly  appear/disappear  within  the  model  system  without  advecdon.  These 
equations  ensure  that  individual  particles  within  die  model  do  not  suddenly  change  species,  that 
is,  a  proton  becomes  an  electron  or  vice  versa.  Before  we  derive  these  particular  conservation 
equations,  we  need  to  define  die  following  variables. 


q  =  Total  Flux  of  Charge  Particles  through  die 
Surface  of  a  Volume. 

=  Average  Number  of  Charged  Particles  per  Unit 
Volume. 

dx  =  Unit  Volume, 
da  =  Unit  Surface  Area. 

1  =  Length  of  Model  Volume. 

w  =  Width  of  Model  Volume, 
h  =  Height  of  Model  Volume. 


We  want  to  get  an  expression  that  shows  the  total  sum  of  a  particular  particle  (ion  (n-j)  or 
electron  (r^))  over  a  given  volume  is  equal  to  the  average  number  density  (n^)  multiplied  by 
the  dimensions  of  the  volume.  These  expressions  are  shown  in  Eqs.  (53  and  54). 


Jy(n1)dT-jy(n#)dr«(n0)hlw,  (53) 

or, 

/v  (  Dj  )  dt  -  ^  (  t^  )  dv  «  (  n„  )  h  1  w  .  (54) 


Taking  Eq.  (53),  let  us  arbitrarily  state 


(55) 
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where  Q  is  a  variable  defined  by  this  equation  and  does  not  necessarily  have  any  other  physical 
meaning.  Rewriting  Eq.  (53)  in  terms  of  Q  gives  us  the  following: 


/v  (  %  )  dt 


/,<f 

C(^)lwdz: 

1  w  /ofc  (  dQ  )  ; 

1  w  [  Q  \l ; 


••  fv  C  ni  )  dt 


lw[Q(h)  -Q(0)]  . 


(56) 


Substituting  Eq.  (53)  into  Eq.  (56)  gives  us  fixe  following: 


(Ho)hlw  =  lw[Q(h)  -  Q  (  0  )  ]  ; 
(n0)h  =  Q(h)  -  Q  (  0  )  . 


(57) 


Assuming  Q(h)  and  Q(0)  have  the  same  magnitudes,  they  must  have  the  following  values 


Q,(M  -  +  ; 

Q*  ( 0 )  =  0. 


(58) 


The  equations  above,  when  used  in  conjunction  with  the  other  equilibrium  equations,  ensure 
that  the  particles  contained  within  the  model  system  do  not  suddenly  change  species. 
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Minimizing  Model’s  Total  Energy  Equation 

We  now  have  a  set  of  differential  equations,  derived  from  the  appropriate  equilibrium 
equations,  which  govern  the  values  of  each  term  within  the  total  energy  equation  (Eq.  (19)).  In 
order  to  find  which  values  establish  a  minimum  total  energy  of  the  model  system,  we  must  solve 
these  differential  equations  using  the  calculus  of  variations. 

Calculus  of  Variations  Calculations 

The  first  thing  we  must  do  is  adjoin  the  constraint  equations  using  a  set  of  Lagrangian 
multipliers.  This  is  done  by  multiplying  die  constraint  equations  (written  in  the  form  0  =  ...)  by 
a  Lagrangian  multiplier  variable  (Eqs.  (59  -  64)),  and  adding  the  results  to  the  total  energy 
equation  (Eq.  (65  -  65e)).  The  constraint  equations  with  their  corresponding  Lagrangian 
multipliers  are 


dE 

V  -  IT  - 

dB. 

V  -  ~  W' 


(59) 


(60) 


dn, 

dz 


*  ^ < E- -  1 


(61) 


X 


dn^ 

dz 


(62) 


(63) 
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The  modified  energy  equation  is  as  follows: 


( — )  U  -  £ { I  \  *  %■ +  ♦  CJn2T  ] 

xy  Jo  2  2  2 

(64) 

dE, 

(64a) 

dB 

*  V  1  *  ”<,r  1 

(64b) 

* ‘ '  ,>B- > 1 

(64c) 

.  i  r  1  _  2-t  p  -i 

‘MTr  q;"  E-J 

(64d) 

*  \  i  'x  *  *■ ) 1  * 

(64e) 

We  now  have  the  total  energy  equation  (Eq.  (19))  written  as  a  function  of  a  set  of  energy  terms 
(Eg,  By  nj,  n2,  Qz,  Vy,  and  Lagrangian  multipliers  and  constants).  In  order  to  find  the 
minimum  total  energy  of  the  model  system,  we  must  find  die  particular  set  of  values  of  the 
energy  terms  that,  when  plugged  into  the  total  energy  equation  (Eq.  (19)),  yield  the  lowest 
possible  value  for  U. 

We  start  the  process  by  first  assuming  we  know  the  value  of  U  when  die  system  is  in  a 
minimum  total  energy  state.  Then  we  look  for  the  particular  set  of  values  that  will  give  us 
assumed  value  of  U.  We  do  this  by  first  expressing  the  dependent  variables  (Eg,  Bx,  nj,  n2,  Qz, 
VJ  as  a  sum  of  a  ’true'  function  and  a  Variation'  function.  The  Variation'  function  itself  is 


composed  of  an  arbitrary  component  multiplied  by  an  infinitesimal  dimensionless  parameter  a . 
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The  ’true’  function  is  defined  as  having  the  value  required  to  minimize  Eq.  (19)  while  the 
Variation’  function  is  defined  as  having  any  value. 


[  Variation  Function  ] 
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Expanding  Eq.  (64  -64e)  in  this  form  gives  us  Eq.  (65  -  65e). 

We  can  see  that  if  we  eliminate  all  die  Variation’  functions  from  Eq.  (65  -  65e)  we  would  have 
the  necessary  values  Ctrue’  functions)  of  the  energy  terms  to  give  us  the  minimum  total  energy 
or  assumed  value  of  U.  Unfortunately,  the  Variation’  function  can  have  any  value  due  to  its 
arbitrary  component,  so  the  only  way  we  can  eliminate  it  is  to  reduce  the  parameter  a  to  zero. 
We  can  do  this  by  taking  the  limit  of  Eq.  (65  •  65e)  as  a  goes  to  zero. 

However,  since  we  are  looking  for  a  minimum  total  energy  condition  of  the  model  system 
when  it  is  also  in  an  equilibrium  state,  we  must  ensure  die  resulting  expression  depicts  an 
equilibrium  condition  after  we  take  the  limit  as  a  goes  to  zero.  We  saw  from  Chapter  II  of  this 
diesis  that  an  equilibrium  condition  exists  when  d/da  equals  zero.  Therefore,  in  order  to 
maintain  all  the  conditions  of  the  model  system  that  we  want,  we  must  also  take  the  derivative 
of  Eq.  (65  -  65e),  with  respect  to  a  at  the  same  time  that  we  take  the  limit  as  a  goes  to  zero. 

The  results  of  performing  these  two  mathematical  operations  yield  Eq.  (66  -  66e).  Notice  that 
the  total  energy  of  the  model  system  is  now  expressed  as  a  derivative  with  respect  to  a 
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1 

( xy ) 


J  o  2 

BXJ  -  2BX«8BX  +  a  *8BX2 


(  5,  ♦  a 8a,  )  (  vy2  +  2vya8Vy  +  a  J8vy2  ) 
2 

+  C,  (  a,  +  «8n,  )  T 


♦  Cj  (  nj  «■  a8nj  )  T  ] 


(65) 


d(  EL  +  a8R,  )  _ 

►  X ^  [ - — —  ^  ■  +  (  n,  ♦  «8n,  )  -  (  Bj  *  aSn,  )  ] 


(65a) 


<K  B  ♦  a8B  ) 

►  \  g - —  ♦  (  n,  +  «8n,  )  (  Vy  +  «8Vy 


)] 


(65b) 


d(  E,  ♦  a8n,  )  (  E,  ♦  «8n,  )a_T 

►  A.  [ - - - —  +  — - - - - 

dz  C,y 

(  ^  ♦  a8^  )  -  (  vy  *  a8Vy  )  (  5,  ♦  a8Bx  ) 
C,Y 


(65c) 


X*[ 


d(  Ej  +  aSnj  )  (  Ej  +  a 8a,  )*** 


dz 


(  E  +  a8E  )  , 
x  — — - —  ] 


c* 


(65d) 


►  \  I  +  (  “i  +  a4nt  )  1  J  <k  • 


(65c) 
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(^)^<iu<a,)  ■  /.‘,[E-SE- 

♦  B.fiB, 

v  2fin,  _ 

+  — — -  +  n.v.ov 

2  1  y  y 

+  C,y  Hi7'1  6% 

♦  CjY  V'  «»2  1 

(66) 

d(  5E  ) 

+  k ^  [ - +  814  -  614  ] 

(66a) 

.  r  d(  8B  )  _  _ 

[-  <fax  ♦  %fivy  +  vy5nj  ] 

(66b) 
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+  fin,  1 )  dz  . 
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(dU/da)'  7}^  expression  depicts  how  much  any  function  varies  from  the  assumed  minimum 
total  energy  value,  U.  When  dU/d«  equals  zero,  this  variance  will  also  be  zero.  Therefore,  if 
we  can  solve  Eq.  (66  •  66e)  so  that  dU/da  equals  zero,  we  will  find  die  values  of  die  energy 
terms  necessary  to  depict  the  minimum  total  energy  state  of  the  model  system. 

If  we  now  reduce  the  value  of  the  Variation'  function  to  zero,  dU/da  will  also  become  zero, 
thereby  giving  us  a  minimum  total  energy  condition  for  the  model  system. 

Taking  a  closer  look  at  Eq.  (66  -  66e)  reveals  that  the  value  of  dU/da  depends  on  the  values 
of  the  arbitrary  components  (jex.  8Bx»  fiiij.  8nj»  and  8Qg)  of  the  Variation’  function.  This 
poses  a  problem  for  us  because  the  value  of  these  arbitrary  components  can  be  anything. 
Therefore,  we  must  somehow  separate  or  factor  out  these  arbitrary  components  from  die  rest  of 
the  equation  so  that  we  are  left  with  only  the  part  of  die  equation  we  need. 

Integration  by  Parts 

We  use  integration  by  parts  to  separate  the  arbitrary  components  from  the  rest  of  the 
expression  in  Eq.  (66  -  66e)  as  shown  here. 


Let  8X 
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Integral 

f  t  d  (  6X  ) 
J  dz 


Any  Arbitrary  Component ; 


Endpoints 
-  *  (  «X  ) 


+ 


Separated 

Integral 

d  (  X  )  . 
dz 


«X  )  . 


(67) 


Notice  the  term  outside  die  integral  on  the  right-hand  side  of  Eq.  (67)  depicts  die  endpoints  or 
boundaries  of  a  given  function,  x  in  this  case. 
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Applying  integration  by  parts  to  Eq.  (66  -  66e)  gives  us  an  equation  Eq.  (68  -  68e)  in  which 
the  energy  terms  are  written  in  the  form  as  shown  in  Eq.  (67).  By  doing  this,  we  can  solve  the 
remaining  part  of  the  equation  to  find  out  what  values  of  the  energy  terms  will  give  us  zero 
variance  (dU/da  =  0). 


.h  d(V) 

-VE« +  /0  [— +  ‘  V"*  1  *  • 


(68a) 


rt  d(]L) 

-  VB*  +  /o  ^ . 8Bx  +  Xal5i6v,  +  Vr8n«  1  *  ; 


(68b) 


-Mn,  +  ♦  ^^  ^E.-v^  -  B^  ) 

1  1  Jo  dz  ^  CjY 


C2-r)5^ 

ClY 


(68c) 


,  .  fh  ,<*(  *2  ),_  *A*‘T«* 

+  Jo  l~^~^ - ^ - 


Xa(  2-y  )  n21'T6n2Et 
CjY 


]  dz  : 


(68d) 


rfc  V- 


-V^  +  fo  +  8n«  J  ^ 


(68e) 


Consolidating  terms  yields: 
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■  ■ 1 V®*  *  VB-  *  V">  *  V”»  *  i! 

*/„'<[  E,«E,  .  B.8B,  .  ♦  %V,«V, 

+  CiY(n,)  T",6n,  +  CjY(«b)  t_lfin2 

d(V> 

♦  5E,  ♦  X^fin,  -  Xqfin, 
d(X.) 

♦  -57-“.  *  V>sv,  *  V»8nt 

,  £i!>s.  .  *,<  “i  eE.  -  v,tB.  -  B,cv.  > 

dz  ”*  CjY 

+  X,(  2  -  T  )  (  a,  )^(  ^  )  (  E,  -  ) 

ClY 

dCAj)  Xj(  iij  )  2‘T(  6E,  ) 

+  - ^ - 

2  -  Y  )  (  n,  )  "-*(  da,  )  (  E,  ) 

d(Xa) 

♦  — +  A^an,  ]  )  dz  . 

(69) 


It  is  important  to  note  here  that  one  of  the  constraints  required  by  the  calculus  of  variations  is 
that  the  value  of  any  function,  "true*  or  ’‘variation/  must  be  the  same  at  any  defined  endpoints 
of  boundaries  of  a  given  system.  Therefore,  the  first  line  of  Eq.  (69)  will  equal  zero.  In  order 
for  dU/da  to  equal  zero,  the  remaining  differential  equations  in  Eq.  (69)  must  have  die 


following  forms: 
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<K\) 

dz 


E  X,(ni)^  A2(na)a-t 

^  Cjr  qr  ’ 


(70) 


dz 


-  B. 


.  *t(  ■,)**(  V,  )  , 


C.Y 


(71) 


da,) 

dz 


\  "  ^Q,  • 


(72) 


d(A1) 

dz 


Mi. 


(73) 


dz 


=  0, 


(74) 


in  order  for  dU/da  to  equal  zero.  Additionally,  the  following  velocity  term  must  equal 


’  ^ *  Aa 


C»Y 


(75) 


We  now  have  a  set  of  differential  equations  represented  by  Lagrangian  multipliers  (Eq.  (70  - 
74))  in  addition  to  Eqs.  (59  -  63)  which,  when  solved,  will  provide  us  with  the  values  of  the 
energy  terms  that  will  minimize  the  minimum  total  energy  equation  (Eq.  (19)).  The  next  step  is 
to  solve  this  set  of  differential  equations  to  find  out  what  the  value  of  the  energy  terms  Ez,  Bx, 
nj,  1*2,  Qz,  and  Vy  should  be  in  order  to  make  dU/da  equal  zero  in  Eq.  (69).  This  is  where  a 
computer  program  becomes  useful. 
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CHAPTER  IV 
MODEL  PROGRAM 

Program  Description 

As  stated  in  Chapter  III,  we  now  have  a  set  of  ordinary  differential  equations  (see  Chapter  III) 
which  describe  an  idealized  model  of  the  earth's  magnetotail.  We  want  to  solve  these 
differential  equations  with  boundary  conditions  to  find  a  particular  set  of  values  of  die  energy 
terms  (Ez,  Bx,  nj,  n2,  Qz,  and  vy)  such  that  when  put  into  die  total  energy  equation  (Eq.  19), 
they  yield  the  lowest  possible  value  for  U.  This  then  becomes  the  minimum  total  energy  model 
for  the  magnetotail  system.  Attempting  to  solve  these  equations  analytically  would  prove 
difficult,  if  not  impossible.  Therefore,  we  must  use  a  computer  program  to  solve  these  equations. 

DVCPR  Routine 

In  order  to  quickly  and  efficiendy  solve  the  differential  equations  derived  in  Chapter  III,  we 
have  used  an  IMSL  flntemational  Mathematics  Science  Library)  computer  routine  -  DVCPR.  The 
DVCPR  program  solves  a  system  of  ordinary  differential  equations  with  boundary  conditions  at 
two  points,  by  using  a  variable  order,  variable  step  size,  finite  difference  method  with  deferred 
corrections. 

As  stated  above,  the  boundary  conditions  for  the  idealized  magnetotail  model  occur  at  die  top 
and  bottom  surfaces  of  the  rectangular  slab.  Basically,  the  DVCPR  tries  to  solve  die  differential 
equations  by  using  the  trapezoidal  rule  over  die  entire  height  of  the  magnetotail  model.  The 
accuracy  of  the  programs  calculations  increases  as  the  number  of  step  sizes  increases.  These 
step  sizes  represent  uniform  horizontal  slices  within  the  model.  The  details  of  how  this  program 
works  are  shown  in  Appendix  A. 
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Dimensional  Transformation 

We  now  have  a  program  which  will  solve  the  system  of  dimensionless  differentia]  equations 
we  derived  in  Chapter  III.  Before  we  enter  these  equations  into  the  program,  we  must  ensure 
the  dimensionless  forms  of  the  energy  terms  relate  to  each  other  in  the  same  manner  as  the 
dimensioned  terms  in  Eq.  (16)  do.  We  do  this  with  scaling  quantities  as  shown  in  Eq.  (76) 
below. 


Dimensioned  Value  -  Dimensionless  Value  x  Scaling  Quantity  . 


(76) 


In  order  to  determine  the  values  of  the  scaling  quantities,  we  rewrite  Eq.  (16)  and  die 
constraint  equations  in  the  following  form: 


m  (  n  )  (  vy  ) 
2 


Energy  Equation 

(B,2)  t  €,(£?) 

2p0  2 


*  (  C  )  (  n»  ) 


(77) 


(Note:  n  =  nj  or  n2  and  C  =  Cj  or  C^O 


Maxwell's  Equations 

d(Et)  q  (  n,  )  q  (  14  ) 
d(z)  %  "  ’ 

and  , 

3(  B„  ) 

■5777  ' 


(78) 


(79) 
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Force  Equations 


3QH)  _  (  n,  )*-*  (  E,  -  v,Bx  )  (8°J 

3<z)  (C,)y 

and  , 

a(n»)  ^  (n,)2-*  (E.)  (81) 

3(z)  "  (Cj)y 


Charge  Conservation  Equation 


3(Q,) 

3(0 


q  ( %  )  • 


(82) 


The  next  step  is  to  multiply  each  of  the  equations  above  (Eq.  (77  -  82))  by  a  corresponding 
scaling  quantity. 


Energy  Equation 


m  (  n  n,  )  (  vy  vy>  ?  +  (  B,  f  gp  (  ^  E^  )* 


2nfl 


(  C  C,  )  (  n  n,  V  . 


(83) 


Maxwell's  Equations 

3(^E^)  q  ( °i  ni, )  q  (°a°2.)  (84) 

3  (  z  z,  )  «o  «o 

and  , 

a  (  B,  Bs  ) 
a  (  z  r. ) 


-  h>  q  ( “i  \ )  ( v7  \  >  • 


(85) 
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3  (  n,  n^  ) 

3  (  *  *,  ) 


3  (  *>2  \  ) 
3  (  z  ^  ) 


Force  Equations 

(  n,  )**  (  (  E,  )  -  (  vy  vy  )  (  Bx  )  ) 

(C,  C^)  Y 

and  , 

(njii^-ME,  V 
(C2^)y 


(86) 


(87) 


Charge  Conservation  Equation 


3  ( Q,  ) 

3(h) 


q  (  n,  n^  )  . 


(88) 


In  Eqs.  (84  -  88),  we  can  collect  all  the  non-scaled  energy  terms  on  the  left  side  of  the 
equation(s)  and  set  them  equal  to  any  arbitrary  value  we  want.  In  this  case,  we  will  set  them 
equal  to  one.  Now  we  set  the  remaining  constants  and  scaled  quantities  on  the  right  side  of  die 
equations  equal  to  the  left  side  of  the  equations  to  get  the  following  equations: 


Maxwell's  Equations 


1 


1 


q  (“r.Hfi) 

+  —  - - - ; 

*0  (V 

q 

«o  <  V  ’ 

and  , 


<  “t.  >  <  vn  >  (  *.  > 

(  > 


(89) 


1 


(90) 
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Force  Equations 

Y  (  ) 

,  _  1  <v<v. 

Y  (Cla)  ' 
and  , 

J  _ 

Y  (<\) 


(91) 


(92) 


Charge  Conservation  Equation 


1 


+  q 


(\)  K) 

<V 


(93) 


We  basically  perform  die  same  type  of  operation  for  the  energy  equation  (Eq.  (83))  with  a  slight 
twist  The  first  thing  we  do  is  factor  the  constants  and  scaled  variables  from  the  first  term  on 
the  right  side  of  Eq.  (83).  We  then  divide  the  rest  of  the  terms  on  the  right-hand  side  of  the 
equation  by  the  same  factor  as  shown  in  Eq.  (94). 


u  = 


Energy  Equation 

(  X  j°v,‘  )  <  ^  X  jB-’  >  (  X  fe’  > 

(  mV,’  )  (  Wr!  )  <  nm»vr.2  ) 


(  Cat  X  Ca  T  ) 

♦  ^  , - -  .(94) 

(  mn*vy<  ) 
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But  since, 


<  5”',’  >  <  jB,1  )  (  iE,1  ) 

(  mn,vy*  )  (  unyr,2  )  (  rnn.v,2  ) 


(CnT) 

(  majr*  ) 


1 

(  mn,vyf  )  ' 


(95) 


we  can  say. 


m^vT. 


2£  . 

2p0 


C.n.T. 


(96) 


At  this  point  we  have  eight  scaling  equations  (Eqs.  (89  -  93,  96))  made  up  of  five  constants 
(<*.  Po.  t»  e,  and  m)  and  seven  unknown  scaling  variables  (Ej^,  B^,  ns,  zy  Q^,  Cs,  and  vy).  It 

is  important  to  note  that  we  are  talking  about  scaling  quantities  here,  not  scaled  quantities.  In 
the  thesis,  "scaling*  quantities  are  used  to  generate  dimensioned  parameters  by  multiplying  them 
with  a  given  dimensionless  number,  and  are  not  to  be  confused  with  a  "scale"  height  or  length 
that  represents  an  interval  in  which  some  parameter  changes  by  a  factor  of  e. 

Setting  one  of  the  unknown  scaling  variables  to  a  specific  value  with  units  (i.e.  vs  =  3.0  x  108 
m/s),  we  can  begin  rewriting  the  remaining  unknown  scaling  variables  in  terms  of  constants  and 
one  unknown  variable,  which  we  arbitrarily  pick  to  be  n.  Since  we  are  choosing  v$  to  be  equal 
to  die  speed  of  light,  we  will  denote  it  as  a  constant  c,  from  now  on. 

Therefore,  we  begin  to  solve  for  the  values  of  die  scaling  quantifies  by  starting  with 


^  «  x  . 


(97) 
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We  can  see  from  Eq.  (90)  that 


-  ^n0  n,  c  2  m  .  (98i 

At  this  point,  we  can  solve  for  any  of  the  scaling  quantities  from  Eqs.  (89  -  93).  We  arbitrarily 
solve  for  the  scaling  height  z$>  from  Eq.  (89).  Also  using  Eqs.  (97)  and  (98),  we  see  that 

^  «o(  V  _  «o  [  c  (  ^  Cp  [  c  (  n,  c  2  m  )  ]  (99) 

'  9«^  9^  9H, 


.2 

If  we  now  square  Eq.  (99)  and  factor  out  c  «  (  p0  ^  )  *»  we  get 


2 


m 

l»o  9  *  “» 


(100) 


Taking  the  square  root  of  Eq.  (100)  yields 


— ® —  ,  (101) 

»»«92n, 

which  is  written  in  terms  of  constants  and  one  scaling  quantity,  n$.  We  are  now  at  a  point 
where  we  can  rewrite  die  remaining  scaling  variables  in  terms  of  ns,  z$,  c,  and  other  known 
constants.  They  are  as  follows: 
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e*  -  f-  (  n.  O  ; 
B*  * 

Q,,  -  e  (  n,  ^  )  : 


C.  - 


«o  Y 


(102) 


We  now  have  all  the  scaling  quantities  written  in  terms  of  constants  and  two  scaling 
quantities,  n$  and  Zg.  Since  the  scaling  height,  Zg,  is  also  written  in  terms  of  the  scaling  number 
density,  rig,  all  the  scaling  quantities  are  related  to  each  other.  Choosing  a  value  for  iig  will  set 
the  values  for  the  other  scaling  quantities  accordingly.  At  this  point,  we  are  able  to  use  the 
dimensionless  forms  of  die  energy  terms  within  die  minimum  total  energy  program  and  use 
DVCPR  routine.  The  details  of  die  minimum  total  program  are  shown  in  Appendix  B. 

Model  Output 

By  picking  a  suitable  scaling  number  density,  we  can  generate  appropriate  values  for  the 
remaining  scaling  quantities.  Multiplying  these  scaling  quantities  by  dimensionless  numbers 
provides  the  input  for  the  DVCPR  program  to  begin  solving  die  dimensionless  differential 
equations  we  derived  in  Chapter  III.  The  program  then  divides  die  model’s  N-S  profile  into  a 
number  of  segments  and  comes  up  with  dimensionless  values  for  each  of  the  energy  variables  for 
each  segment.  These  values  can  then  be  plotted  on  a  graph  to  come  up  with  dimensionless  N-S 
profiles  of  die  energy  variables.  The  profiles  generated  by  this  program  only  depict  the  region 
between  die  bottom  surface  of  the  model,  or  neutral  sheet,  to  die  top  surface  or  die  plasma 
sheet  boundary  layer.  Although  die  southern  half  of  the  central  plasma  sheet  is  not  depicted  in 
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this  model,  the  profiles  should  by  symmetrical  from  north  to  south  for  the  number  densities,  ion 
velocities,  and  ion  temperatures.  Due  to  die  nature  of  die  magnetotail,  the  electric  and  magnetic 
fields  would  have  asymmetrical  profiles. 

By  setting  the  dimensionless  number  density  of  the  model  magnetotail  to  one  and  setting  the 
scaling  number  density  to  1.0  X  106  particles  per  cubic  meter  (typical  value  observed  within  the 
central  plasma  sheet  of  the  earth’s  magnetotail),  we  come  up  with  die  following  scaling 
quantities  for  the  other  model  energy  parameters. 


Scaling  Quantities 


n*  ->  (  1.0  X  106  [m‘3]  ) 
z,.  ->  ( 2.28  X  10s  [m]  ) 

Eg  ->  (4.12X103  [V/m]) 

Bs  ->  (  1.37  X  10'5  [t]) 

Cg  ->  (1.5  X  1014  [ (kg) (m4)/ (A) (s2)]  ) 


Next,  we  multiply  these  scaling  quantities  by  one  to  come  up  with  the  following  model  energy 
parameters,  with  units. 


Raw  Data 

= 

Dimensionless 

Value 

X 

Scaling  Quantity 

nRAW  > 

1.0  X106 

[m3] 

(1.0) 

X 

(  1.0  X  106  [m*3]  ) 

2  RAW  "> 

2.28  X  10s 

[m] 

(1.0) 

X 

(  2.28  X  10s  [m]  ) 

Eraw  "> 

4.12  X  103 

[V/m] 

(1.0) 

X 

(  4.12  X  103  [V/m]  ) 

®RAW  *> 

1.37  X  10  s 

[t] 

(1.0) 

X 

(  1.37  X  10  s  [t]  ) 

CRAW  •> 

1.5  X1014 

[kgm4/A3]  = 

(1.0) 

X 

(  1.5  X  1014  [kgm4/A3]  ) 

These  numbers,  when  entered  into  the  minimum  total  energy  program  (Appendix  B),  generate 


integrated  numerical  data  (Appendix  C)  and  a  dimensionless  N-S  profile  as  shown  in  Fig.  24. 
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Notice  that  the  ion  and  electron  number  densities  both  increase  as  they  get  closer  to  the  model 
magnetotail’s  neutral  sheet,  which  is  what  we  would  expect  to  see  in  the  earth's  magnetotail. 
Also  note  that  the  magnetic  field  strength  goes  to  zero  at  the  model  magnetotail’s  neutral  sheet- 
again  what  we  would  expect  to  see  in  die  earth’s  magnetotail. 

Figures  25  through  28  show  how  die  dimensionless  energy  parameters  vary  as  the  height  of 
the  model  increases  from  0.01,  a  very  thin  slab,  to  95  scaling  heights,  a  thick  slab  (numerical 
data  are  in  Appendix  C).  It  is  interesting  to  note  that  all  die  parameter  changes  occur  within 
three  (6.84  x  10**  m)  or  four  scale  heights  (9.12  x  10"*  m),  no  matter  how  high  we  make  the 
profile.  In  order  to  ensure  the  results  generated  by  the  minimum  energy  model  are  not 
inconsistent  with  the  laws  of  classical  physics,  we  compare  die  height  in  which  die  parameter 
changes  take  place  to  the  Debye  length  and  the  Larmor  radius  of  an  ion  used  within  the  model. 
The  equations  for  finding  an  ion’s  Debye  length  and  Larmor  radius  are 


Debye  Length  »  XD 

Larmor  Radius  =  r 


T 

69.0  (  —  )w  ; 

“h. 

(  “t,  X  Vfc,  ) 

(  q  X  B  ) 


Using  the  numerical  data  from  Appendix  C,  we  see  that  the  ion  Debye  length  is  approximately 
1542  m  while  the  Larmor  radius  ranges  in  value  from  2.29  x  10s  -  5.62  x  10^  m.  We  can 
ignore  the  Debye  length  since  it  is  smaller  than  the  Larmor  radius.  Looking  at  the  Larmor 
radius,  however,  we  can  see  that  its  values  have  the  same  order  of  magnitude  as  the  height  in 
which  the  parameters  change  in  the  model  profiles.  This  result  shows  that  the  minimum  energy 
model  is  not  inconsistent  with  die  laws  of  classical  physics. 


Magnetotail  Height  Magnetotal  Height 
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FIG.  24.  N-S  profiles  of  number  densities  (n,  and  iij),  x-component  of  the  magnetic  field  (Bx),  and 
the  z-component  of  the  electric  field  (E,,)  generated  from  the  minimum  total  energy  program. 


TABLE  II.  Scaling  Data  Corresponding  to  Fig.  (24). 


Scaling  Quantity 

Number  of 
Scaling  Quantities 

rij  (Number  Density) 

=  1.0  ; 

Zj  (Model  Magnetotail  Height) 

=  1.0; 

Eg  (Z-component  of  the  Electric  Field) 

-  1.0; 

Bs  (X-component  of  the  Magnetic  Field) 

=  1.0; 

Cs  (Magnetotail  Adiabatic  Constant) 

=  1.0  . 
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FIG.  25.  N-S  profiles  of  number  densities  (nj  and  i^),  x-component  of  die  magnetic  field  (Bx),  and 
the  z-component  of  the  electric  field  (iy  generated  from  die  minimum  total  energy  program. 


TABLE  III.  Scaling  Data  Corresponding  to  Fig.  (25). 


Scaling  Quantity 

Number  of 
Scaling  Quantities 

iij  (Number  Density) 

=  1.0; 

Zg  (Model  Magnetotail  Height) 

=  0.01  ; 

Eg  (Z-component  of  the  Electric  Field) 

=  1.0; 

Bs  (X-component  of  the  Magnetic  Field) 

=  1.0; 

Cs  (Magnetotail  Adiabatic  Constant) 

=  1.0  . 

Magnetotall  Haight  Magnetotan  Height 
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FIG.  26.  N-S  profiles  of  number  densities  (nj  and  r^),  x-component  of  the  magnetic  field  (Bx),  and 
the  z-component  of  the  electric  field  (Ej)  generated  from  the  minimum  total  energy  program. 


TABLE  IV.  Scaling  Data  Corresponding  to  Fig.  (26). 


Scaling  Quantity 

Number  of 

Scaling  Quantities 

rij  (Number  Density) 

=  1.0; 

z,  (Model  Magnetotail  Height) 

=  0.1  ; 

Eg  (Z-component  of  the  Electric  Field) 

=  1.0; 

Bs  (X-component  of  the  Magnetic  Field) 

=  1.0 ; 

Cs  (Magnetotail  Adiabatic  Constant) 

=  1.0. 

Magnatotail  Height  Magnetotail  Height 
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FIG.  27.  N-S  profiles  of  number  densities  (nj  and  nj,  x-component  of  the  magnetic  field  (Bx),  and 
file  z-component  of  the  electric  field  (E^)  generated  from  the  minimum  total  energy  program. 


TABLE  V.  Scaling  Data  Corresponding  to  Fig.  (27). 


Scaling  Quantity 

Number  of 

Scaling  Quantities 

rij  (Number  Density) 

=  1.0; 

Zj  (Model  Magnetotail  Height) 

=  10.0  ; 

Ej  (Z-component  of  the  Electric  Field) 

=  1.0; 

Bs  (X-component  of  the  Magnetic  Field) 

=  1.0  ; 

Cs  (Magnetotail  Adiabatic  Constant) 

=  1.0  . 

1.02 
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FIG.  28.  N-S  profiles  of  number  densities  (nj  and  n2),  x-component  of  the  magnetic  field  (Bx),  and 
the  z-component  of  the  electric  field  (Ez)  generated  from  the  minimum  total  energy  program. 


TABLE  VI.  Scaling  Data  Corresponding  to  Fig.  (28). 


Scaling  Quantity 

Number  of 
Scaling  Quantities 

(Number  Density) 

=  1.0; 

Zj  (Model  Magnetotail  Height) 

=  95.0 ; 

Ej  (Z-component  of  the  Electric  Field) 

-  1.0; 

Bs  (X-component  of  the  Magnetic  Field) 

=  1.0; 

Cs  (Magnetotail  Adiabatic  Constant) 

=  1.0. 

1.00 

1.02 
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In  order  to  come  up  with  real  values  that  represent  typical  values  of  the  earth’s  magnetotail, 
we  would  have  to  multiply  the  scaling  quantities  by  the  following  dimensionless  values. 


Raw  Data 

Dimensionless 
=  Value 

X 

Scaling  Quantity 

nRAW  *> 

1.0  X106 

[nf3] 

=  (  1.0  ) 

X 

(  1.0  X  106  [m'3]  ) 

ZRAW  m> 

2.96  X  107 

[m] 

=  (  1.30  X  102  ) 

X 

(  2.28  X  10s  [m]  ) 

ERAW  '> 

1.5  X10'5 

[V/m] 

=  (  3.64  X  10'9  ) 

X 

(  4.12  X  103  [V/m]  ) 

®RAW  *> 

2.0  XIC8 

[t] 

=  (  1.46  X  10  3  ) 

X 

(  1.37  X  10‘5  [t]) 

CRAW 

6.9  X1019 

[kgm4/A3] 

=  (  4.59  X  10'5  ) 

X 

(1.5  X  1014  [kgm4/A3]  ) 

Unfortunately,  the  minimum  total  energy  program  is  unable  to  handle  these  numbers  because 
the  system  height  (1.30  x  102  scale  heights)  is  much  larger  than  the  scaling  height  z,..  However, 
by  starting  with  a  small  height,  unity  for  example,  and  increasing  it  until  the  program  can  no 
longer  handle  the  input  values,  we  can  look  at  die  trends  in  the  numerical  data  (Appendix  Q 
and  profiles  (Figs.  29  -  34)  and  draw  conclusions  from  them. 

Remember,  the  top  and  bottom  values  of  each  of  these  profiles  are  constrained  by  the 
following  boundary  conditions: 


X,  (0) 

=  0  ; 

Mh) 

=  0; 

MO) 

=  0  ; 

X2(h) 

=  0; 

E,(0) 

=  0  ; 

Ef(h) 

-  0; 

Bx  (  0  ) 

=  0  ; 

B,(h) 

"  +B*. 

Magnetotail  Height  Magnetotail  Height 
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FIG.  29.  N-S  profiles  of  number  densities  (n}  and  i^),  x-component  of  the  magnetic  field  (Bx),  and 
the  z-component  of  the  electric  field  (Ez)  generated  from  the  minimum  total  energy  program. 


TABLE  VII.  Scaling  Data  Corresponding  to  Fig.  (29). 


Dimensioned  Parameter 

— 

Dimensionless 

Unit 

X 

Scaling 

Quantity 

nRAW 

1.0  X106 

[m3] 

(1.0) 

X 

(  1.0  X  106  [m*3]  ) 

ZRAW  "> 

2.28  X  10s 

[m] 

=: 

(1.0) 

X 

(  238  X  10s  [m]) 

E-raw  -> 

1.5  X10'5 

[V/m] 

s 

(  3.64  X  10'9  ) 

X 

(  4.12  X  103  [V/m]  ) 

®RAW  '> 

2.0  X10-* 

M  4  , 

= 

(  1.46  X  10'3  ) 

X 

(  1.37  X  10  s  [T]  ) 

CRAW  *> 

6.9  X1019 

[kgm4/A3] 

sss 

(  4.59  X  10'5  ) 

X 

(1.5  X  10'14  [kgm4/A3]  ) 
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FIG.  30.  N-S  profiles  of  number  densities  (n}  and  rij),  x-component  of  the  magnetic  field  (Bx),  and 
the  z-component  of  the  electric  field  (Ej)  generated  from  the  minimum  total  energy  program. 


TABLE  VIII.  Scaling  Data  Corresponding  to  Fig.  (30). 


Dimensioned  Parameter 

. 

Dimensionless 

Unit 

X 

Scaling 

Quantity 

nBAW  *> 

1.0  X106 

Inf3] 

(1.0) 

X 

(  1.0  X  106  [in3]  ) 

ZRAW  ■> 

4.56  X  10s 

[m] 

(2.0) 

X 

(  2.28  X  10s  [m]  ) 

eraw  *> 

1.5  X10'5 

(V/m) 

= 

(  3.64  X  10'9  ) 

X 

(  4.12  X  103  [V/m]  ) 

®RAW  '> 

2.0  XIO4* 

[t] 

= 

(  1.46  X  10'3  ) 

X 

(  1.37  X  10-®  [t]  ) 

^RAW  > 

6.9  X1019 

[kgm4/A3] 

= 

(  4.59  X  10  s  ) 

X 

(1.5  X  10'14  [kgm4/A3] ) 

Magnetotail  Height  Magnetotail  Height 
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FIG.  31.  N-S  profiles  of  number  densities  (n}  and  n^,  x-component  of  the  magnetic  field  (Bx),  and 
the  z-component  of  the  electric  field  (E,,)  generated  from  the  minimum  total  energy  program. 


TABLE  IX.  Scaling  Data  Corresponding  to  Fig.  (31). 


Dimensioned  Parameter  = 

Dimensionless 

Unit 

X 

Scaling 

Quantity 

nRAW  *> 

1.0  X  106 

[in'3] 

(1.0) 

X 

(  1.0  X  106 

[nr3]  ) 

ZRAW  '> 

6.84  X  10s 

[m] 

(3.0) 

X 

(  2.28  X  105 

[m]) 

ERAW  *> 

1.5  X10  5 

[V/m] 

(  3.64  X  10’9  ) 

X 

(  4.12  X  103 

[V/m]) 

ERAW  ‘> 

2.0  XIO* 

W  ,  ,  = 

(  1.46  X  10*3  ) 

X 

(  1.37  X  10  s 

It]  ) 

Craw  m> 

6.9  X  1019 

[kgm4/A3]  = 

(  4.59  X  10*  ) 

X 

(1.5  X1014 

[kgm4/A3] ) 
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FIG.  32.  N-S  profiles  of  number  densities  (nj  and  rij),  x-component  of  die  magnetic  field  (Bx),  and 
the  z-component  of  die  electric  field  (iy  generated  from  the  minimum  total  energy  program. 


TABLE  X.  Scaling  Data  Corresponding  to  Fig.  (32). 


Dimensioned  Parameter  = 

Dimensionless 

Unit 

X 

Scaling 

Quantity 

nRAW  "> 

1.0  X106 

[m-3] 

(1.0) 

X 

(  1.0  X  106 

[m*3]  ) 

ZRAW  "> 

9.12  X  10s 

[m] 

(4.0) 

X 

(  2.28  X  10s 

[m]) 

^RAW  '> 

1.5  X10‘s 

[V/m] 

(  3.64  X  10'9  ) 

X 

(  4.12  X  103 

[V/m]) 

®RAW  '> 

2.0  X10-* 

W  A  ,  = 

(  1.46  X  10‘3  ) 

X 

(  1.37  X  105 

It]  ) 

CRAW 

6.9  X1019 

[kgm4/A3]  = 

(  4.59  X  10  s  ) 

X 

(  1.5  X  1014 

[kgm4/A3]  ) 
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FIG.  33.  N-S  profiles  of  number  densities  (nj  and  n£,  x-component  of  the  magnetic  field  (Bx),  and 
the  z-component  of  die  electric  field  (Ez)  generated  from  the  minimum  total  energy  program. 


TABLE  XI.  Scaling  Data  Corresponding  to  Fig.  (33). 


Dimensionless  Scaling 

Dimensioned  Parameter  -  Unit  X  Quantity 


nRAW  ■> 

1.0  X106 

[nf3] 

(1.0) 

X 

(  1.0  X  106  [m‘3]  ) 

ZRAW  '> 

1.14  X  106 

[m] 

(5.0) 

X 

(  2.28  X  10s  [m]  ) 

ERAW  _> 

1.5  X10'5 

[V/m] 

(  3.64  X  10‘9  ) 

X 

(  4.12  X  103  [V/m]  ) 

®RAW  ‘> 

2.0  X10*8 

M  = 

(  1.46  X  10'3  ) 

X 

(  1.37  X  10  s  [t]  ) 

CRAW  '> 

6.9  X10'19 

[kgm4/A3]  = 

(  4.59  X  10'5  ) 

X 

(1.5  X  10'14  [kgm4/A3]  ) 
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FIG.  34.  N-S  profiles  of  number  densities  (nj  and  n£,  x-component  of  the  magnetic  field  (Bx),  and 
the  z-component  of  the  electric  field  (E^)  generated  from  the  minimum  total  energy  program. 


TABLE  XII.  Scaling  Data  Corresponding  to  Fig.  (34). 


Dimensionless  Scaling 

Dimensioned  Parameter  =  Unit  X  Quantity 


nRAW  ‘> 

1.0 

X  106 

[m'3] 

(1.0) 

X 

(  1.0  X  106  [m*3]  ) 

ZRAW  *> 

2.28  X  106 

[m] 

(  10.0  ) 

X 

(  2.28  X  10s  [m]  ) 

^RAW  _> 

1.5 

X  10'5 

(V/m] 

(  3.64  X  10‘9  ) 

X 

(  4.12  X  103  [V/m]  ) 

®RAW  "> 

2.0 

xio-8 

(t] 

(  1.46  X  10‘3  ) 

X 

(  1.37  X  10  s  [t]  ) 

CRAW  *> 

6.9 

X  1019 

[kgm4/A3]  = 

(  4.59  X  10  s  ) 

X 

(1.5  X  10'14  [kgm4/A3]  ) 
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Comparison  of  these  calculated  N-S  profiles  to  those  observed  in  the  earth’s  magnetotail,  as  well 
as  the  overall  assessment  of  the  energy  minimizing  technique  presented  in  this  thesis,  will  be 
discussed  in  the  next  chapter. 
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CHAPTER  V 
DATA  ANALYSIS 

Observed  Magnetotail  Parameters 

The  ultimate  goal  of  executing  die  minimum  total  energy  program  (see  Appendix  2)  is  to 
come  up  with  values  for  the  ion  and  electron  number  densities  (n^  and  n2)  and  for  die  electric 
and  magnetic  (Bx)  field  strengths  over  the  height  of  the  model  magnetotail  system. 

Plotting  these  on  a  graph  yields  N-S  profiles  of  these  model  parameters.  By  comparing  the 
model’s  calculated  profiles  with  those  obtained  from  typical  values  of  the  earth’s  magnetotail,  we 
should  be  able  to  determine  how  effective  the  minimum  total  energy  model  is  at  defining  a 
minimum  total  energy  state  of  the  earth’s  magnetotail. 

Figure  35  depicts  a  N-S  profile  of  some  key  parameters  typical  of  the  near  earth’s  magnetotail 
(20  -  30  Rg)  during  solar  quiet  times  and  with  a  southward  interplanetary  magnetic  field  (IMF). 
In  order  to  compare  the  profiles  in  Fig.  (35)  to  the  ones  generated  by  the  minimum  energy 
model,  we  must  first  linearize  the  logarithmic  scales  for  the  number  density,  flow  speed,  and  ion 
temperatures  (Figs.  36  -  38).  Since  this  N-S  profile  does  not  change  dramatically  for  a 
northward  pointing  IMF  (Lui,  1987),  we  will  assume  that  the  minimum  energy  model  depicts 
this  situation  also. 


Comparison  to  Model  Output 

At  first  glance,  the  model  generated  profiles  do  not  compare  very  well  to  the  profiles  obtained 
from  averaged  parameters  observed  in  the  earth’s  magnetotail,  but  there  are  some  similarities. 
Looking  at  the  model  generated  profiles  using  typical  magnetotail  values  (Figs.  29  •  34),  except 
for  height,  and  noting  that  we  are  only  looking  at  the  central  plasma  sheet  region,  we  can  make 
some  comparisons  between  the  two  profiles.  Notice  that  in  both  profiles,  model  generated  and 
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FIG.  35.  North-south  profiles  of  magnetic  field  (B),  number  density  (n),  plasma  bulk  flow  speed  (v), 
and  ion  temperature  (Tj)  in  Ae  midnight  region  of  the  magnetotail  [Lui,  1987]. 


FIG.  36.  Number  density  (n)  profile  from  Fig.  (35)  using  a  linear  number  density  scale. 
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observed,  the  ion  and  electron  number  densities  increase  rapidly  near  the  upper  edges  of  the 
boundary  (plasma  sheet  boundary  layer  for  the  observed  profile).  The  number  densities  then 
increase  less  rapidly  as  they  approach  the  center  or  bottom  boundary  (neutral  sheet)  of  each 
profile.  The  model  differs  from  the  obsei  ved  values  in  that  there  is  no  significant  increase  in  the 
number  densities  towards  the  center  of  the  system  (or  neutral  sheet). 

Another  similarity  occurs  in  the  magnetic  field  strength  profiles.  The  observed  magnetic  field 
strength  appears  to  lose  strength  more  quickly  in  the  top  half  of  the  central  plasma  sheet  (3  -  5 
Rg)  than  in  the  lower  half  (0  -  3  Rg).  The  same  situation  occurs  in  the  model  magnetotail,  but 
the  magnetic  field  strength  rate  changes  more  rapidly  in  the  upper  half  of  the  m  Jel  and  less  so 
(if  at  all)  in  the  lower  half  of  the  profile.  This  difference  becomes  more  dramatic  with  increased 
scale  height  (Figs  33  -  34). 

The  velocity  and  temperature  trends  for  both  the  magnetotail  and  model  profiles  also  show 
some  similarities.  For  instance,  both  systems  have  their  maximum  velocities  at  their  upper 
boundaries,  which  decrease  until  they  reach  their  minimum  at  the  lower  boundaries.  The 
temperature  trends,  on  the  other  hand,  show  the  opposite  trend  for  both  systems.  Both  systems 
have  their  lowest  temperatures  at  the  upper  boundaries,  which  increase  to  a  maximum  at  the 
lower  boundaries. 

Finally,  both  the  magnetotail  and  the  model  systems  both  display  a  degree  of  diamagnetism. 

In  other  words,  both  systems  have  a  tendency  to  try  to  get  rid  of  the  magnetic  field  energy  from 
their  centers.  Looking  at  the  profiles  for  both  systems  shows  that  the  magnetic  energy  is  at  a 
minimum  near  the  bottom  boundary  or  center  of  each  system.  In  the  magnetotail  system,  this 
magnetic  field  energy  increases  slowly  away  from  the  center,  then  about  halfway  (3  Rg)  between 
the  neutral  sheet  and  the  plasma  sheet  boundary  layer,  the  magnetic  field  strength  increases 
more  rapidly.  The  same  thing  happens  in  the  model  system,  but  the  rate  of  magnetic  field 
energy  strength  is  more  dramatic  with  increased  profile  height.  The  comparisons  between  the 
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magnetotail  and  model  systems  are  summarized  in  Table  XIII. 

As  mentioned  earlier,  all  the  energy  parameter  changes  appear  to  be  confined  to  the  top  of  all 
the  N-S  profiles  as  the  model  system’s  height  becomes  greater  than  ten  scaling  heights  (z  =  10  x 
Zj).  It  appears  that  any  further  increases  in  the  profile  heights  will  only  stretch  the  bottom 
regions  of  each  profile  whose  values  already  remain  relatively  constant  at  higher  heights.  This 
phenomenon  seems  to  suggest  that  the  model  is  depicting  a  situation  where  the  majority  of  the 
energy  changes  or  action  takes  place  within  die  outer  boundaries  or  sheaths  of  the  region  being 


Table  XIII.  Profile  Comparison  Summary  Between  Model  and  Magnetotail  Systems. 


Observed 

Model 

Number  density  increases  rapidly 
at  die  Plasma  Sheet  Boundary 

Layer,  less  so  near  the  Neutral 

Sheet. 

Number  density  increases  rapidly 
at  die  top  boundary,  less  rapidly 
at  die  bottom  boundary. 

Magnetic  field  strength  decreases 
quickly  near  the  Plasma  Sheet 

Boundary  Layer,  less  so  near  the 

Neutral  Sheet 

Magnetic  field  strength  decreases 
quickly  near  the  top  boundary, 
less  so  near  the  bottom  boundary. 

Plasma  bulk  flow  speeds  are  greatest 
near  the  Plasma  Sheet  Boundary 

Layer,  lowest  near  die  Neutral 

Sheet 

Ion  velocities  are  greatest  near 
the  top  boundary  of  the  model, 
lowest  near  the  bottom  boundary. 

Ion  temperatures  are  lowest  at  the 

Plasma  Sheet  Boundary  Layer,  highest 
at  die  Neutral  Sheet 

Ion  temperatures  are  lowest  at  die 
top  boundary  of  the  model,  highest 
at  the  bottom  boundary. 

Magnetic  field  displays  some  degree 
of  diamagnetism. 

Magnetic  field  displays  some 
degree  of  diamagnetism. 

Height  of  the  top  half  of  the  Ceqtral 

Plasma  Sheet  is  about  3.19  X 10' m. 

Maximum  height  of  die  model  is 
about  228  X  106  m. 
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modeled.  Closer  inspection  of  die  observed  energy  profiles  in  Fig.  35  shows  similar  trends 
within  or  near  the  Plasma  Sheet  Boundary  Layer  just  as  the  model  seems  to  indicate. 

There  appear  to  be  many  reasons  to  explain  the  differences  observed  between  the  two  profiles, 
and  one  of  the  more  obvious  may  be  the  geometry  of  the  model  system.  Although  the  Central 
Plasma  Sheet  of  the  earth’s  magnetotail  can  be  approximated  by  a  rectangular  slab,  its  height-to- 
width  proportions  are  not  infinite  as  depicted  by  the  model  system.  While  making  the  model’s 
length  and  width  dimensions  infinite  helped  keep  calculations  simple,  it  did  not  accurately 
portray  the  magnetotail  system. 

The  boundary  conditions  placed  on  the  model  system  may  also  have  restricted  it  too  much  to 
accurately  depict  the  magnetotail  system.  For  instance,  the  model  system  constrains  the  z 
component  of  the  electric  field  to  be  zero  at  both  the  bottom  and  top  boundaries.  Although  this 
constraint  keeps  the  model  simple,  it  may  not  reflect  what  is  actually  observed  within  the 
neutral  sheet  and  plasma  sheet  boundary  layer.  Some  of  the  boundary  conditions  may  have 
been  too  simple  due  to  some  of  the  assumptions  made.  For  example,  the  model  assumes  that 
tiie  z  component  of  the  magnetic  field  equals  zero  at  the  bottom  boundary  (or  neutral  sheet), 
which  may  not  be  true  in  the  magnetotail.  Similar  arguments  can  be  made  for  the  other 
boundary  conditions  used  in  the  model. 

Another  possible  reason  for  the  differences  observed  between  the  two  profiles  may  be  that  the 
model  may  have  been  oversimplified  and  does  not  depict  actual  magnetotail  conditions.  For 
example,  the  model  described  in  this  thesis  used  only  an  x  component  of  the  magnetic  field  (no 
y  component).  This  was  done  for  two  reasons.  The  first  was  to  simplify  the  calculations 
required  for  the  minimum  total  energy  program.  The  second  reason  was  that  for  a  first 
approximation,  the  region  that  the  model  was  trying  to  depict  (the  midnight  region  of  the 
central  plasma  sheet)  the  y  component  of  the  magnetic  field  appeared  to  be  negligible  compared 
to  the  x  component  We  cannot  say  the  same  tiling  about  the  z  component  of  the  magnetic 
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field.  Because  the  Central  Plasma  Sheet  is  the  likely  site  where  anti-parallel  magnetic  field  lines 
would  connect  between  the  northern  and  southern  tail  lobes,  the  z  component  of  the  magnetic 
field  does  exist  within  the  earth’s  magnetotail.  This  z  component  of  the  magnetic  field  would 
probably  have  an  appreciable  effect  on  the  outcome  of  die  model’s  output. 

Another  model  oversimplification  was  that  it  only  allowed  one  species  of  particles,  ions,  to 
move  about  within  die  model  system.  Again,  this  assumption  was  used  to  keep  the  calculations 
for  the  minimum  total  energy  simple.  This  assumption  may  have  affected  the  electron  number 
density  profile.  Given  the  same  initial  parameters  in  the  computer  program,  die  electron 
number  density  profiles  should  have  looked  like  die  ion  number  density  profiles. 

Finally,  the  model  system  was  assumed  to  be  in  a  steady-state  condition.  Observations  of  die 
earth’s  magnetotail  suggest  that  it  is  rarely,  if  ever,  in  a  steady-state  condition.  Therefore,  we 
have  a  situation  where  we  are  trying  to  compare  a  steady-state  model  system  to  a  non  steady- 
state  magnetotail.  If  it  were  possible  to  get  some  profiles  of  the  earth’s  magnetotail  parameters 
while  it  experienced  a  steady-state  condition,  however  briefly,  we  might  see  some  additional 
correlations  between  the  model  and  actual  magnetotail  systems. 


COPY  [LAB:  *  TC»  mr  TV)EP  flO*  T’ERMr*  ft  *4-  l.BG*~  ' 
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CHAPTER  V 
CONCLUSIONS 


Although  the  minimum  total  energy  model  presented  in  this  thesis  is  able  to  depict  many  of 
the  trends  observed  within  the  central  plasma  sheet  of  the  earth’s  magnetotail,  it  is  evident  that 
it  does  not  have  the  sufficient  detail  to  depict  it  accurately.  One  similarity  that  the  model  does 
have  with  the  earth’s  magnetotail,  however,  is  that  almost  all  of  the  parameter  changes  occur 
near  the  outer  or  top  boundary  of  the  model  system.  The  same  thing  is  true  in  the  earth’s 
central  plasma  sheet  in  that  the  greatest  parameter  changes  also  occur  near  the  outer  boundary 
or  plasma  sheet  boundary  layer.  Another  similarity  between  the  model  and  the  earth’s  central 


plasma  sheet  is  that  both  systems  show  diamagnetic  qualities.  In  other  words,  both  systems 
show  indications  of  trying  to  rid  their  centers  of  magnetic  field  energy.  \ 

Therefore,  the  model  seems  to  suggest  that  in  order  for  the  central  plasma  sheet  in  the  eartfaw 
magnetotail  to  reach  a  minimum  total  energy  state,  it  must  attempt  to  eliminate  as  much  enerfy 
from  the  center  of  its  system  as  possible.  The  model  demonstrates  this  fact  very  well  with  th. 
magnetic  field  energy.  It  appears  that  the  model  system  attempts  to  reach  a  minimum  total 
energy  state  by  keeping  the  magnetic  field  energy  to  a  minimum  near  the  center,  and  allowing 

t 

its  energy  to  increase  near  die  outer  boundary.  There  is  some  indication  that  the  earth’s 
magnetotail  system  also  shows  Some  indication  of  trying  to  keep  die  magnetic  field  strength  near 
the  neutral  sheet  to  a  minimum. 

These  similarities  seem  to  indicate  that  although  die  model  may  be  oversimplified  in  its 
present  form,  it  does  have  a  handle  on  the  basic  physics  governing  die  energies  located  wit^ 
the  central  plasma  sheet  However,  it  is  also  evident  that  further  modifications  are  still,  nee&tf 
to  make  this  minimum  total  energy  model  useful. 


wr  »  aoc  a*  ,,,  ^ 


Yunna  4 

/x  '  b 
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Areas  for  Future  Studies 

One  of  die  first  steps  to  improving  this  model  may  be  to  eliminate  or  modify  some  of  the 
assumptions  used,  specifically  those  discussed  in  Chapter  V.  Other  areas  of  improvement 
include: 

a)  changing  the  geometry  of  die  model  to  better  correspond  to  die  shape  of  the  central 
plasma  sheet  located  within  the  earth’s  magnetotail. 

b)  using  boundary  conditions  that  more  accurately  depict  conditions  observed  within 
the  plasma  sheet  boundary  layer. 

c)  adding  a  z  component  of  the  magnetic  field  into  die  model  system. 

d)  adding  a  capability  into  the  model  to  include  a  gradient  up  and  down  die  x-axis  of 
magnetotail  model.  In  this  way,  the  model  would  be  able  to  take  into  account  die 
processes  occurring  up  and  down  the  magnetotail  axis. 

e)  including  additional  processes  observed  within  the  earth's  magnetotail  into  the 

model. 

f)  comparing  model  output  with  profiles  generated  from  die  earth’s  magnetotail  during 
northward  IMF. 

Although  this  model,  in  its  present  form,  does  not  precisely  depict  the  earth’s  magnetotail 
system,  it  may  represent,  to  a  first  approximation,  other  systems  in  non-thermodynamic 
equilibriums  (i.e.  radiating  systems). 
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PURPOSE 

USAGE 

ARGUMENTS  - 


IMSL  ROUTINE  NAME  -  DVCPR 
(As  reproduced  from  IMSL  Library 
FORTRAN  Subroutines  for  Mathematics  and  Statistics) 


Solve  a  system  of  ordinary  differential  equations  with  boundary  conditions 
at  two  points,  using  a  variable  order,  variable  step  size  finite  difference 
method  with  deferred  corrections. 

Call  DVCPR  (N,  FCNI,  FCNJ,  FCNB,  XA,  XB,  NGMAX,  NGRID,  IP,  IR,  TOL,  X, 
Y,  IY,  ABT,  PAR,  WORK,  IWORK,  IER). 

N  -  Number  of  differential  equations  (INPUT). 

FCNI  -  Name  of  subroutine  for  evaluating  derivatives  (INPUT).  The 

subroutine  itself  must  also  be  provided  by  the  user  and  it  should 
be  of  the  following  form: 

SUBROUTINE  FCNI  (N,  X,  Y,  YPRIME) 

REAL  Y(N),  YPRIME (N)  ... 

FCNI  should  evaluate  YPRIME(l)  ...  YPRIME (N)  given  N,  X,  and 
Y(l)  ...  Y(N).  YPRIMEO)  is  the  derivative  of  Y(I)  with  respect  to 
X. 

FCNI  must  appear  in  an  external  statement  in  the  calling  program. 

FCNJ  -  Name  of  the  subroutine  for  evaluating  the  N  by  N  Jacobian  matrix 
of  partial  derivatives  (INPUT).  The  subroutine  itself  must  also  be 
provided  by  the  user  and  it  should  be  of  the  following  form: 

SUBROUTINE  FCNJ  (N,  X,  Y,  PD) 

REAL  Y(N),  PD(N,  N)  ... 

FCNJ  should  evaluate  PD(I,  J)  for  I,  J  =  1,  N,  given  N,  X,  and  Y(l) 
...  Y(N).  PDQ,  j)  is  the  partial  derivative  of  YPRIMEO)  with 
respect  to  Y(J). 

FCNJ  must  appear  in  an  external  statement  in  die  calling  program. 

FCNB  •  Name  of  the  subroutine  for  evaluating  the  boundary  conditions 
ONPUT).  The  subroutine  itself  must  also  be  provided  by  the  user 
and  it  should  be  of  the  following  form: 


SUBROUTINE  FCNB  (N,  YA,  YB,  F) 
REAL  YA(N),  YB(N),  F(N)  ... 
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XA.XB  - 

NGMAX  - 
NGRID  - 

IP 

IR 

TOL 

X 

Y 

IY 

ABT 

PAR 


FCNB  should  evaluate  F(l)  ...  F(N)  given  YA(1)  ...  YA(N),  YB(1)  ... 
YB(N).  YAG)  and  YB(I)  are  the  values  of  YG)  at  XA  and  XB, 
respectively,  and  the  boundary  conditions  are  defined  by 
FQ)  =  00.0, 1  =  1,  N.  The  initial  conditions  must  be  defined 
first,  then  the  coupled  conditions,  and  then  the  final  conditions. 

FCNB  must  appear  in  an  external  statement  in  the  calling 
program. 

Two  points  where  boundary  conditions  are  given  GNPUT).  XA 
must  be  less  than  XB. 

Maximum  number  of  grid  points  to  be  allowed  GNPUT). 

Number  of  points  in  the  input  grid  (counting  endpoints).  NGRID 
must  be  greater  than  3.  On  output,  NGRID  will  contain  the  final 
number  of  grid  points  GNPUT/OUTPUT). 

Number  of  initial  conditions  GNPUT).  IP  must  be  greater  than  or 
equal  to  0  and  less  than  N. 

Number  of  coupled  boundary  conditions  GNPUT).  IP  +  IR  must  be 
greater  than  0  but  less  than  or  equal  to  N. 

Relative  error  control  parameter  GNPUT).  The  computations  stop 
when  ABS  (ERROR  (J,  I))/AMAX1  (ABO  (Y(J,  I),  1.0)  is  less  than 
TOL  for  all  J  =  1  ...  N,  1  =  1  ...  NGRID,  where  ERROR  (J,  1)  is  the 
estimated  error  in  Y(J,  I). 

Vector  of  length  NGMAX  containing  the  final  grid.  If  PAR(4)  =  0, 
the  program  initialize  X  to  a  uniform  mesh  of  NGRID  points. 
Otherwise  the  user  must  supply  the  initial  grid  on  input 
GNPUT/OUTPUT) . 

Matrix  of  dimension  N  by  NGMAX  containing  the  computed 
solution  on  the  final  grid.  Y(J,  I)  will  return  an  approximation  to 
the  Jth  solution  component  at  XG).  If  PAR(4)  —  0,  the  program 
initialize  Y  to  zero.  Otherwise  the  user  must  supply  initial  values 
for  Y  GNPUT/OUTPUT). 

Row  dimension  of  matrix  Y,  exactly  as  specified  in  rite  dimension 
statement  GNPUT).  IY  must  be  greater  than  or  equal  to  N. 

Vector  of  length  N  containing,  in  its  Jth  component,  an  estimate  of 
the  maximum  absolute  error  over  rite  grid  points  for  rite  Jth 
solution  component  (OUTPUT). 

Options  vector  of  length  5  GNPUT).  If  PAR(l)  =  0  die  default 
options  are  used  and  the  remaining  components  are  ignored.  If 
PAR(l)  =  1,  all  remaining  components  of  PAR  must  be  given  a 
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value.  The  default  value  of  PAR (I)  in  each  case  is  zero. 

PAR(2)  greater  than  0  implies  that  continuation  is  to  be  done  for 
this  highly  nonlinear  problem.  It  is  assumed  that  the  user  has 
embedded  his  problem  in  a  one  parameter  family 
DY/DX  =  YPRIME  (X,  Y,  EPSNU) 

F(Y(A),  Y(B),  EPSNU)  =  0 

such  that  for  EPSNU  =  0,  the  problem  is  simple  (e.g.  linear),  and 
for  EPSNU  =  1,  the  original  problem  is  recovered.  The  program 
will  automatically  attempt  to  go  from  EPSNU  =  0  to  EPSNU  =  1. 
PAR(2)  is  the  starting  step  in  the  continuation.  The  step  may  be 
varied  by  DVCPR  but  a  lower  bound  the  stepsize  of  0.01  is 
imposed.  The  following  common  block  should  appear  in 
subroutines  FCNI,  FCNJ,  and  FCNB  - 

COMMON  /  Cl  /  EPSNU,  CONT 
REAL  EPSNU 
LOGICAL  CONT 

If  CONT  =  .TRUE,  vectors  YPRIME  in  subroutine  FCNI  and  F  in 
subroutine  FCNB  should  be  defined  by 

YPRIME  (I)  =  D  (YPRIME  (I))/D  (EPSNU) 

FO)  =  D  (F  fl))/D  (EPSNU) 

and  when  CONT  =  .FALSE.,  YPRIME  and  F  should  have  their 
normal  definitions. 

PAR(3)  =  1,  implies  that  intermediate  output  is  to  be  printed  (for 
debugging  purposes). 

PAR(4)  =  1,  implies  that  initial  values  for  X  and  Y  are  supplied  by 
the  user. 

PAR(5)  =  1,  implies  that  the  differential  equations  and  boundary 
conditions  are  linear,  and  the  algorithm  should  take  advantage  of 
this  fact 

WORK  -  Real  work  vector  of  length 

N*  (3  *N*NGMAX+4*N + 1) + NGMAX*  (7*N + 2) 

IWORK  -  Integer  work  vector  of  length 

2*N*NGMAX+N+NGMAX 

IER  -  Error  parameter  (OUTPUT). 

Terminal  error 
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IER  =  129  Illegal  values  for  N,  NGRID,  IP,  or  IP. 

IER  =  130  More  than  NGMAX  grid  points  are  needed  to  solve 
the  problem. 

IER  =  131  Newton’s  iteration  diverged. 

IER  =  132  Newton’s  iteration  reached  roundoff  error  level.  If 
requested  precision  is  not  attained,  this  means  that 
TOL  is  too  small. 
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APPENDIX  B:  Minimum  Total  Energy  Program 
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MINIMUM  TOTAL  ENERGY  PROGRAM 


REAL  Z(1001),  K(10,  1001),  WORK(384824),  ZA,  ZB,  PAR(5),  ABT(3),  TOL,  U(1001), 
ENERGY(lOOl),  VX(1001),  ZBR(lOOl),  NlR(lOOl),  N2R(1001),  EZR(lOOl),  BYR(lOOl), 
QR(1001) 

REAL  MRAW,  QRAW,  NRAW,  ZRAW,  ERAW,  BRAW,  CRAW,  VUGHT,  EPSILON,  MU,  NS,  ZS,  ES, 
BS,  QS,  CS,  VS,  MQN,  BOLT,  TRAW 

INTEGER  IWORK(21031),  N,  IK,  IER,  NGMAX,  NGRID,  IP,  IR 

INTEGER  I,  MRATIO 

COMMON/XXX/VAR,  C.  Cl,  C2,  GAM,  Gl,  G2,  ZA,  ZB,  DENSITY,  BO,  NO 
EXTERNAL  FCNI,  FCNJ,  FCNB 


************************************************************************************* 

**  OPEN  DATA  FILES  ** 

************************************************************************************* 


OPEN  (1,  FILE  -  TILENAME.DAT,  STATUS  - ’NEW’) 
OPEN  (2,  FILE  =  ’PLOT-NlJ>LT,  STATUS  -  ’NEW’) 
OPEN  (3,  FILE  -  ’PLOT-N2.PLT,  STATUS  -  ’NEW) 
OPEN  (4,  FILE  =  ’PLOT-EFID-PLT’,  STATUS  =  TIEW) 
OPEN  (5,  FILE  -  ’PLOT-BFLD.PLT,  STATUS  «  ’NEW) 
OPEN  (6,  FILE  =  ’PLOT-Q.PLT,  STATUS  «*  ’NEW’) 


[This  part  of  the  program  opens  files 
storing  numerical  and  plotting  data] 


ft************************************************************************************ 


INPUT  INITIAL  VALUES  (HAVING  DIMENSIONS) 


WRITE  (*,  *)  ’Enter  number  density  (in  particles  per  cubic  meter)’ 

READ  (*,*)  N^ 

[measured  n} 
ornj] 

WRITE  (*,  *)  ’Enter  height  of  magnetotaQ  (in  meters)’ 

READ  (*,*)  Zraw 

[measured  z] 

WRITE  (*,  *)  ’Enter  z-component  of  the  electric  field  (in  volts  per  meter)’ 
READ  (*,*)  Eraw 

[measured  E,] 

WRITE  (*,  *)  ’Enter  x-component  of  the  magnetic  field  (in  teslas)’ 

read  (*,*)  braw 

[measured  Bx] 

WRITE  (*,  *)  ’Enter  magnetotail  temperature  (in  degrees  Kelvin)’ 

READ  (*,  *)  Tmw 

[measured  T] 
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**  CONSTANTS  ** 

************************************************************************************* 


VUGHT 

GAM 

EPSILON 

MU 

BOLT 

Mraw 

Qraw 

CRAW 


3.0  X  10® 

5.0/ 3.0 
8.85  X  1012 
1.257  X  10"6 
1.381  X  10'23 
1.67  X  10'27 
1.602  X  1019 

(BOLT)  *  CTraw)*(Nraw  X  10°-°^) 


[Velocity  of  Light] 

[Gamma  =  Ratio  of  Specific  Heats] 
[Permittivity  of  Free  Space] 
[Permeability  of  Free  Space] 
[Boltzmann  Constant] 

[Mass  of  Ion  in  kgs] 

[Electron  Charge] 

[Calculated  Adiabatic  Constant] 


a************************************************************************************ 

**  SCALING  QUANTITIES  ** 

************************************************************************************* 


Nraw 


=  Inputted  from  above 


[Measured  Number  Density] 


NS 

Zs 

Es 

BS 

Qs 

Cs 


nraw  _ 

Sqrt  (Mraw  /  (MU  *  NRAW  *  Qraw  )) 
(Qraw  *  Ns  *  Zs)  /  (EPSILON) 

(Eg)  /  (VUGHT) 


(Qraw*Ns*Z?) 

(VUGHT2)  *  NSX1 


-GAM)  * 


Mraw 


[Scaled  num  den  =  measured  num  den] 
[Scaled  magnetotail  ht] 

[Scaled  E  fid  strength] 

[Scaled  B  fid  strength] 

[Scaled  Chg  conserve  term] 

[Scaled  adiabatic  constant] 


**  DIMENSIONLESS  VALUES  ** 

t************************************************************************************ 


Zraw 

eraw 

®raw 

Craw 


=  Inputted  from  above  (measured  magnetotail  height  in  meters) 

=  Inputted  from  above  (measured  electric  field  in  volts  per  meter) 
=  Inputted  from  above  (measured  magnetic  field  in  teslas) 

=  Calculated  from  above  (adiabatic  constant) 


N0 

Zo 

Eo 

B0 

Co 

Ci 

C2 

G, 

g2 


“  NRAW  /  Ns 
=  Zraw  /  Zg 

=  Eravv  /  Eg 

=  Braw/Bs 

=  CRAW/  Cs 
-  Co 

=  Co 

=  1  /  (Cj  *  GAM) 
=  1  /  (Cj  *  GAM) 


[Num  den  w/  units] 
[Magnetotail  ht  w/  units] 

[E  fid  w/  units] 

[B  fid  w/  units] 

[Adiabatic  constant  w/  units] 


[Adiabatic  constant  for  ions] 
[Adiabatic  constant  for  elections] 
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*******************************************  ****************************************** 

**  DVCPR  ARGUMENT  SETUP  ** 

************************************************************************************* 


N 

= 

10 

[Number  of  variables] 

Za 

= 

0 

[Height  of  bottom  surface  of  model] 

zB 

= 

zo 

[Height  of  top  surface  of  model] 

NGMAX 

= 

1001 

[Number  of  grid  points} 

NGRID 

= 

100 

[Number  of  horizontal  slices  in  model] 

IP 

= 

5 

[Number  of  initial  conditions] 

IR 

= 

0 

[Number  of  coupled  boundary  conds] 

IK 

= 

10 

[Number  of  output  columns] 

TOL 

= 

0.01 

[Relative  error  parameter] 

PAR(l) 

= 

1 

[No  defaults] 

PAR(2) 

= 

0 

PAR(3) 

= 

0 

[No  printed  output  for  debugging] 

PAR(4) 

= 

1 

[Values  for  X  and  Y  supplied  by  user] 

PAR(5) 

= 

0 

[Diff  Eqs.  and  B.C.S  are  not  linear] 

************************************************************************************* 

**  INITIAL  GUESS  ** 

************************************************************************************* 


RMM1 

= 

NGRID  -  1 

[This  part  of  the  program  divides  the 
magnetotail  into  evenly  spaced  ZQ) 

DO  101 
RIM1 

ZQ) 

= 

1,  NGRID 

I  -  1 

ZA  +  RIM1/RMM1  *  (ZB  -  ZA) 

horizontal  planes] 

K(l,  I) 

= 

NO 

[Estimated  ion  num  density  (N)  at  ht  I] 

K(2, 1) 

= 

NO 

[Estimated  electron  num  density  at  ht  I] 

K(3,I) 

= 

0 

[Estimated  electric  field  at  ht  I] 

K(4,0 

ZO^CBo/Zb) 

[Estimated  magnetic  field  at  ht  I] 

K(5,I) 

= 

Z0)  *  (N0  *  Zb)  /  Zb 

[Estimated  charge  conservation  at  ht  I] 

K(6, 1) 

= 

0 

[Lagrangian  equiv  of  n,  at  ht  I] 

K(7,0 

= 

0 

[Lagrangian  equiv  of  n2  at  ht  I] 

K(8, 1) 

= 

0 

[Lagrangian  equiv  of  E,  at  ht  I] 

K(9. 1) 

= 

0 

[Lagrangian  equiv  of  Bx  at  ht  I] 

K(10,  I) 

= 

0 

[Lagranfian  equiv  of  Qz  at  ht  I] 

10  CONTINUE 
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************************************************************************************* 

**  CALL  THE  IMSL  SUBROUTINE  DVCPR  ** 

************************************************************************************* 

CALL  DVCPR  (N,  FCNI,  FCNJ,  FCNB,  ZA,  ZB,  NGMAX,  NGRID,  IP,  IR,  TOL,  Z,  K,  IK,  ABT,  PAR, 
WORK,  IWORK,  IER) 

IF  (IER  .GT.  100)  GOTO  200 


************************************************************************************* 

**  TRANSFORM  DIMENSIONLESS  VARIABLE  TO  REAL  VALUES  ** 

************************************************************************************* 


DO  8001 

*  1,  NGRID 

[Real  Value  =  Dimensionless  value 

x  scaled  quantity] 

ZbrG) 

-  Z(I)  *Z<; 

[Yields  real  magnetotail  ht  at  ht  I] 

N1Ra) 

-  K(l,  I)  *  Nraw 

[Helds  real  ion  num  density  at  ht  I] 

N 2R0) 

=  K(2, 1)  *  NMVi 

[Yields  real  elect  num  density  at  ht  I] 

EzrCI) 

=  K(3,I)*ES 

[Yields  real  E  fid  strength  at  ht  I] 

ByrG) 

=  K(4, 1)  *  Bs 

[Yields  real  B  fld  strength  at  ht  I] 

QrG) 

=  K(5, 1)  *  Qs 

800  CONTINUE 


************************************************************************************* 

**  FORMAT  THE  DVCPR  OUTPUT  ** 

************************************************************************************* 

Once  the  system  of  ordinary  differential  equations  are  calculated  by  this  program  we  want  the 
out  put  to  be  formatted  as  follows: 


Interval 

Value  of 

Z  at  Ht  m 

Value  of 
nj  at  Ht  (I) 

Value  of 
n2_atHLQl 

Value  of 

Et  at  Ht  (I) 

Value  of 

B„  at  Ht  m 

Z(I) 

K  (1, 1) 

K  (2, 1) 

K  (3, 1) 

K(4,I) 

K  (5, 1) 

Value  of 

Value  of 

Value  of 

Value  of 

Value  of 

Interval 

0T  at  Ht  (I) 

*1  at  Ht  (!) 

A2«tHtg) 

*T  at  Ht  fl) 

A,  at  Ht  fl) 

Z(I) 

K  (5, 1) 

K  (6, 1) 

K(7,I) 

K(8,I) 

K  (9, 1) 

(This  data  can  be  transferred  into  a  plotting  routine  in  order  to  make  N-S  profiles  of  the 
desired  parameters). 
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************************************************************************************* 

**  CALCULATING  THE  MINIMUM  TOTAL  ENERGY  (U)  OF  THE  MAGNETOTAIL  ** 
************************************************************************************* 

DO  141  =  1,  NGRID 

Vy(I)  -  [-K(6, I)  /  Cj  *  GAM)]  *  [  K(l,  I)1_GAM]  *  Bx  +  K(9, 1) 

or  Iv*  "  + 

U(I)  =  0.5  *  K(l,  I)  *  [(V  O))2]  +  0.5  *  K(3, 1)2  +  0.5  *  K(4, 1)2  +  C,  *  K(l,  I)GAM  + 

C2  *  K(2, 1)GAM 


or 


I  U(I) 


E.2  B_2  , 


ENERGY0)  =  U(I) 

WRITE  (1,  145)  ’Energy  at  Ht  [z(I)l  is  [ENERGY0)]’ 

SUM  =  0 
DO  160  I  =1,  NGRID 
SUMNEXT  =  ENERGYQ)  +  SUM 
SUM  =  SUMNEXT 
UTOTAL  =  SUMNEXT 


or  [  U  (  TOTAL  ) 


Z»  ” 


WRITE  (1,  170)  The  minimum  value  of  U  with  a  magnetotail  ht  of  [Zq]  is  [UTOTAL]’ 
200  TYPE  *,  Hun  terminated  because  IER  =  [IER]’ 


END 
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a************************************************************************************ 

**  SUBROUTINE  FOR  CALCULATING  THE  DIFFERENTIAL  EQUATIONS  ** 

************************************************************************************* 


SUBROUTINE  FCNI  (N,  Z,  K,  KPRIME) 

REAL  X,  K(N),  KPRIME (N),  iij,  n2,  Ej,  Bz,  Qz,  Xj,  Ij,  Xg,  Xg,  Xq,  Vy 
INTEGER  N 

COMMON/XXX/VAR,  C,  Cj,  C2,  GAM,  G1(  G2>  ZA,  Z*  DENSITY,  Bq,  N0 


ni 

n2 

Ez 

Bx 

Q* 

*1 

x2 

^z 

*z 


K(l) 

K(2) 

K(3) 

K(4) 

K(5) 

K(6) 

K  (7) 

K(8) 

K(9) 

K(10) 

[  +Xj  *  Gj  *  (  Hj  )1'GAM  *  Bx  -  Xx  ] 


or  [  Vy 


,  *i  (  *4  )I"T  (  ) 

C,Y 


KPRIME(l)  =  I  G,  *  (  nj  l2"0"*  ] 

*  [  Ez  -  (  Vy  *  Bx  )  ] 

KPRIME  (2)  =  [  -G2  *  (  n2  )2_GAM  ]  *  [  Ez  ] 
KPRIME (3)  =  [  n,  -  n2  ] 

KPRIME(4)  =  [  nj  *  v  ] 


0,  j  .  ,  (  °i  )*  *(  E,  ~  v,B. )  1; 

dz  CjY 

or  [  *!  .  .  <^2HLSl>  1; 
dz  CjY 

°r  [  5  =  +  (  “1  “  “a  )  1; 

0F  1  ^  "  *  <  n«v,  >  1: 

or  1  ~  -  ♦(«!)]; 


KPRIME(5)  -  [  nj  ] 
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KPRIME(6) 


KPRIMEC7)  = 


KPRIME(8)  = 


KPRIME(9)  « 


KPRIME(IO) 


[  (  -.5  *  vy2  )  -  (  C,  *  GAM  *  (  nx  )GAM1  ) 

+  (  Xx  *  Gx  *  (  2  -  GAM  )  *  ( rij 

*  (  -  ^  +  (  vy  *  Bx  )  )  -  XE 

*  (  Xx  *  vy  )  -  XQ 

t  -Gj  *  GAM  *  (  n2  )GAM1  ] 

+  [  X2  *  G2  *  (  2  -  GAM  ) 

+  XE 


r  3X.  v  * 

°r  [-ST  "  “  “f  “Cl  Y  04)’-* 
oz  2 

clT 

x  (  -  B,  -  *px  )  -  x^ 

-  (  Vr  >  "  \  1; 

or  1  Ife 

.  ^a(2-Y) 

ClT 

x  (  n,  )  ■-*  (  E,  ) 

♦  \]: 


[  *E*  *  Xj  *  Gj  *  (  nj  )2jGAM  ] 

+  [X2*G2*(n2)2^AM] 


or  [  .  -  «  -  hlhTl 

to  ^  C,» 


♦ 


C,Y 


I  *BX  +  Xj  *Gj  *  (n,  )MAM*vy]  or  [f^5-  .  -  X|  vf  }]; 

dz  CjY 


[0] 


or  [  Q  .  (0)1; 

dc 


RETURN 

END 
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************************************************************************************* 

**  SUBROUTINE  FOR  CALCULATING  THE  JACOBIAN  MATRIX  ** 

************************************************************************************* 


SUBROUTINE  FCNJ  (N,  Z,  K,  PD) 

REAL  Z,  K(N),  PD(N,N),  n,.  n*  Er  Bx,  Qz,  X„  X2,  XE,  Xg,  Xq,  v  ^1,  *L, 

V  r  cb,  3X,  3BX 

INTEGER  N 

COMMON/XXX/VAR,  C,  C„  Gj,  GAM,  G,,  G^  ZA,  Zg,  DENSITY,  Bq,  N0 


K(l) 

K(2) 

K(3) 

KC4) 

K(5) 

K(6) 

K(7) 

K(8) 

K(9) 

K(10) 


[  +X,  *  G,  *  (  n,  *  Bx  -  Xx  ]  or  [  vy 


X,(b,)^(B,)  ]; 

C,Y 


DO  300  I  =  1,  10 
DO  310  J  =  1,10 
PD(I,J)  *  0.0 
310  CONTINUE 

300  CONTINUE 


=  [  +G,  *  X,  *  (  1  -  GAM  )  ] 
*  [  (  1  /  (  n,  )GAM  )  *  Bx  ] 


,  X,  (  1  -  y  )  (  B,  )  ]. 

C,y  (  n,  )  T 


[1  +G,*X1*(n1  J1-0™] 


[  i 


=  [+G,*(n,)1‘GAM*B  ] 


[^L  .  .  (V^LljU]; 


[  -1.0  ] 


or  [  -zr2-  -  -  (  10  )  l; 

d\ 
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PD(1,  1) 


PD(1,  3) 


PD(1,  4) 


PD(1,  6) 


PD(1,  9) 


PD(2,  2) 


PD(2,  3) 


=  [  G,  *  (  2  -  GAM  )  *  (  n,  )1-GAM  or  [  ^ 

*(  WB*>1 

*  [  G,  *  (  n,  *  Bx  )  ] 

eta, 


(  2  -  t  X  n,  )  1-T 

C»Y 


x  (  E,  -  v,Bx  ) 


(n»)2'T(^B«> 

- - ]; 

C,Y 


-  [  G,  *  (  n,  )2<5AM  ] 


C,Y 


=  t -G,  *  (  n,  )2<JAM 


*(BX* 


+  vy)  ] 


Oh)2^ 

C,t 

-  <  b,  x  ^  i; 


dv_ 


=  [  -G,  *  (  n,  )2<:AM  *  ( Bx  *  )  ]  or  [ 


ax, 


ax. 


(n,)2’T(B«)(  } 


]; 


0V_ 


-  ,  <  B,  X  ) 

[-G1*(n1)2-GAM*(B,*  or  [  , - 1; 

3X8,  C,Y 


-  [  -G2  *  (  2  -  GAM )  *  (  n2  )1<}AM  or  [  ^  -  -  (  2  ~  Y  X  n,  )*•*(  E,  ) 

aij  CjT 

*  (Ej)  ] 


=  t-G^Cn,)2™*]  or  .  -ih2!l]; 

3E,  CjY 


PD(3,  1)  =  [  +1.0  ] 
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or  [  ^  b  +  (  i.o  )  ]; 

an, 


PD(3,  2)  =  [  -1.0  ] 


or  [  ^  •=  -  (  1.0  )  3; 

At] 


PD(4,  1)  =  [(n,  *  ^1)  +vy] 

*4 


°r  ‘ 


PD(4,  4) 


-  I  n,  *  OZl  J 


[vD  w*_  i 

■=r  -  («,X^)]: 

od1  a»x 


PD(4,  6)  -  [  n,  *  ^JL  J 

a*. 


r  cB  dv  , 

or  I  — i  -  (  n,  X  -*  )  1; 


PD(4,  9)  =  [  n,  *  ] 


r  5B,  3v  . 

OT‘< 


PD(5,  1)  =  [+1.0  3 


or  [  *  ♦  (  1.0  )  1; 


110 


PD(6,  1) 


PD(6,  3) 


PD(6,  4) 


[(-v*  ) 


or  [ 


GAM-2 


-  (  1/Gj  )  *  (  GAM  -  1  )  *  ( tij  )' 

+  (  G,  *  A,  *  (  2  -  GAM  )  ) 

*  ((  1  -  GAM  )  *  (  fij 

*  (  -E,  +  vy  *  Bx )) 

+  LI  *  (  G.  *  (  n,  *  (  OUl  *  B_  ) 

*  (  2  -  GAM  )  -  (  A-  *  J 

“*  a, 


3A, 

8bj 


3v_ 

-(v,X^) 

-(C.YXY-1XB,)’'1 

+  Mini 

ClY 

x(l  ~yX®i)"t 

x  (  -E,  ♦  v,B,  ) 

♦  At(  n,  )1l^XBl) 

X  (  2  -  Y  )  -  A^(  |i  )]; 


-  [  -G,  *  Aj  *  (  2  -  GAM  ) 


or  .  -M^-TXn,)1^, 

*,  C,Y 


*(ni)WAMl 


-  ((-».*  ±1) 
r  ®. 

or  [ 

SB, 

dv 

+  (  G,  *  A,  *  (  2  -  GAM  ) 

.  M  *  ‘  T  > 

C.Y 

*  (  n,  )1<JAM  *  *  Bx  ) 

X  ( 3  ^  X  B, 

-CA  *^) 

"»  SB, 

3v. 

■**<*> 

+  (  Gj  *  A1  *  (  11,  ),-GAM  *  vy  )  ] 

l,(  2  -  r  x  «, )'"’( v. 

Ill 


PD(6,  6) 


PD(6,  8) 


PD(6,  9) 


[(-v  * 
y 


+  (  G,  *  X,  *  (  2  -  GAM) 


*  (  n,  )1_GAM  *  *  B.  ) 

dX, 


r  3X. 

or  [  _ L 

dX , 


dv 

xt  ( 2  -  r  ) 


C,Y 


3v_ 


*  (  *1  >  ,_T(  ^  X  B,  ) 

dX, 


-  (  Gj  *  (  2  -  GAM  )  ) 


*  (  n,  )1<im  *  (  -Ej  +  vy  *  Bx  ) 

_ 1 1 

B*  »i 


-  ( »  *  )  i 


_  (2  -  Y  ) 

ClY 

x  (  n,  )  »-»(  -E,  ♦  VyB,  ) 


dv  , 
-  x„  (  — *  )  ]; 
B-  ax, 


-  [-1.0] 


or  [ 


»L 

ax„ 


(  i.o  >  1; 


=  [  (  -V, 


* 

ax. 


+  (  G,  *  Xj  *  (  2  -  GAM) 


(n,  )1JGAM*  i*Bx) 


-( 


* 

ax. 


)-Vv 


r  ax. 
or  [  1 


ax^ 


dv 

-(v'x<> 

A,  (  2  -  Y  ) 


C,Y 


dv. 


(n‘)!"T(^XB«) 


3v  . 


r  5X. 

or  [  _ L 

dXn 


PD(6,  10)  =  [  -1.0  ] 


-  ( i  .o )  1; 


PD(7.  2)  =  [  (  -1  /  G2  )  *  (  GAM  -  1) 


*  (  n2  )GAM’2  +  (  X2  *  G2  ) 


*  (  2  -  GAM  )  *  (  1  -  GAM  ) 


*  (  n2  )<ifM  *  Ej  ] 


PD(7,  3)  =  [  G2  *  X2  *  (  2  -  GAM  ) 


*  (  n2  )1-GAM  ] 


PD(7,7)  =  [  G2  *  (  2  -  GAM  ) 


*  (  *  Ez  ] 


P D(7,  8)  =  [+1.0] 


or  [5l 


or  [  121 


-(CjyXy-D 


*<» 2)Y'2  +  FL 


x(2  -  yXI  ~  Y) 
x  (  n*  )  T(  ^  )  ]; 


ax,  X,  (  2  -  Y  X  %  )  ^ 


(2-TXn,)H(Et)]. 

C*Y 
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PD(8,  1) 

=  [  -G,  *  Aj  *  (  2  -  GAM  ) 

or 

[  - 

1,  (  2  -  y  X  n,  ?-* 

Aij 

C,Y 

PD(8,  2) 

=  [G2*  A2*  (2 -GAM) 

or 

[  f!^5»  = 

Aj  ( 2  -  y  x  nj  ?-* 

ALj 

^Y 

*  (  n2  )1<JAM  ] 

PD(8,  3) 

-  [-1.0] 

or 

-  ( i.o )  1; 

PD(8,  6) 

=  [  -Gj  *  (  llj  )2_GAM  ] 

or 

[  9Xb>  « 

C.T 

PD(8,  7) 

II 

r-n 

♦ 

w 

j 

or 

[  _ 

" 

PlY 

dv_ 


PD(9  1)  -  [  Gj  *  A,  *  (  n,  *  (  Z0l  )  or  [  _5»  ,  + _ gj 

So,  Ax,  Cjy 


]; 


+  G,  *  A,  *  (  2  -  GAM  ) 


A,(  2  -  Y  X  »H  )»•*(  v,  )  y 
C,Y 


*  C  n,  )1_GAM  *  vv  ] 
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PD(9,  4)  -  [  -1  +  G,  *  X,  *  (  n,  ) 


*(  ^Z)] 


2-GAM 


or  [ 


ay 


3B. 


CiY 


3; 


PD(9,  6)  -  [  Gj  *  (  nj  ) 


2-GAM 


or  [ 


dk^ 

ax, 


*(X,*  +vy)] 

ax,  y 


*, ( ”*  ^  I? > 

♦ _ “l_ 

C»Y 

+  (Mfl(vL)]; 

C,Y 


dv_ 


PD(9,  9)  -  [  G,  *  X,  *  (  n,  *  J?L  ]  or  [  «  + _ _ _ ]; 


d\ 


ax„ 


C.Y 


RETURN 

END 
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************************************************************************************* 

“  SUBROUTINE  FOR  ENTERING  BOUNDARY  CONDITIONS  ** 

************************************************************************************* 


SUBROUTINE  FCNB  (N,  KA,  KB,  F) 

REAL  KA(N),  KB(N),  F(N) 

INTEGER  N 

COMMON/XXVVAR,  C,  Clf  C^,  GAM,  Gj,  G2>  ZA,  Zg,  DENSITY,  Bq,  N0 


F(l) 

= 

KA(6) 

or 

[  B.C 

at  sfc  A  is 

- 

o  ]; 

F(2) 

= 

KA(7) 

or 

[  B.C 

at  sfc  A  is 

at 

o  ]; 

F(3) 

= 

KA(3) 

or 

[  B.C. 

at  sfc  A  is 

E. 

= 

o  1; 

F(4) 

= 

KA(4) 

or 

[  B.C 

at  sfc  A  is 

B, 

- 

o  1; 

F(5) 

= 

KA(5) 

or 

I  B.C 

at  sfc  A  is 

Q, 

- 

o  ]; 

F(6) 

s 

KB(6) 

or 

[  B.C. 

at  sfc  B  is 

s 

o  ]; 

F  (7) 

= 

KB(7) 

or 

[  B.C 

at  sfc  B  is 

as 

o  ]; 

F(8) 

* 

KB(3) 

or 

t  B.C 

at  sfc  B  is 

K 

X 

o  1; 

F(9) 

= 

KB(4)  -  B0 

or 

t  B.C 

at  afc  B  is 

B, 

- 

+  b.  1; 

F(10) 

= 

KB(5)-(N0*Zb) 

or 

I  B.C 

at  sfc  B  is 

Q, 

m 

+  (  Vb)]: 

(Note:  Surface  A  is  die  bottom  surface  of  die  model  magnetotail  which  in  this  particular  case 
represents  the  center  or  neutral  sheet  of  die  earth's  magnetotail.  Surface  B  is  the  top  of  the 
model  magnetotail  and  represents  the  upper  boundary  of  the  earth’s  magnetotail.) 


RETURN 


END 


APPENDIX  C:  Model  Data 
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Multiplying  Scaling  Quantities  by  One. 


These  tables  show  the  calculated  values  of  each  energy  variable  in  the  minimum  total  energy 
equation.  The  value  of  each  variable  is  added  to  the  values  calculated  from  the  previous  scaling 


height. 


°o  -  *0  EZq  —  1.0;  Bxq  —  1.0;  Cq  —  1.0. 

Integrating  Energy  Parameters  from  0.0  to  1.0. 

I 

zfD 

MOL 

_n,GL 

_SafQ_ 

JUD. 

B  M. 

Vyffi. 

Temofll 

1 

0.00 

2.56 

1.09 

1.01 

0.00 

0.00 

0.85 

1.06 

5 

0.04 

12.81 

5.44 

5.03 

0.01 

0.09 

4.27 

5.30 

10 

0.09 

25.63 

10.88 

10.06 

0.04 

0.42 

8.55 

10.59 

20 

0.19 

51.33 

21.71 

20.12 

0.16 

1.79 

17.17 

21.16 

30 

0.29 

77.21 

32.46 

30.16 

0.35 

4.10 

25.95 

31.66 

40 

0.39 

103.34 

43.06 

40.19 

0.61 

7.38 

34.99 

42.08 

50 

0.50 

129.82 

53.47 

50.21 

0.92 

11.64 

44.37 

52.36 

60 

0.60 

156.72 

63.62 

60.20 

1.25 

16.90 

54.19 

62.48 

70 

0.70 

184.12 

73.43 

70.18 

1.59 

23.17 

64.60 

72.37 

80 

0.80 

212.07 

82.83 

80.13 

1.89 

30.49 

75.73 

81.98 

90 

0.90 

240.57 

91.70 

90.07 

2.11 

38.89 

87.77 

91.23 

100 

1.00 

269.56 

99.94 

99.99 

2.19 

48.37 

100.97 

100.03 

°o  = 

1.0;  Zq  = 

0.01;  EZq  = 

s  1.0;  Bxq 

=  1.0;  Cq 

=  1.0. 

Integrating  Energy  Parameters  from  0.0  to  0.01. 


I 

zfT> 

JJQL 

JliHL 

J120L 

_E*GL 

VyQL 

TemDfll 

1 

0.00 

5491 

l.ll 

1.00 

0.00 

0.00 

99.58 

1.07 

5 

0.0004 

27452 

5.53 

5.00 

0.0001 

0.11 

497.91 

5.36 

10 

0.0009 

54866 

11.06 

10.00 

0.0005 

0.50 

995.82 

10.71 

20 

0.0019 

109406 

22.05 

20.00 

0.0020 

2.11 

1991.44 

21.38 

30 

0.0029 

163278 

32.92 

30.00 

0.0045 

4.83 

2987.46 

31.96 

40 

0.0039 

216150 

43.58 

40.00 

0.0077 

8.63 

3983.29 

42.41 

50 

0.0050 

267703 

53.97 

50.00 

0.0115 

13.50 

4979.12 

52.69 

60 

0.0060 

317633 

64.03 

60.00 

0.0156 

19.39 

5974.95 

62.75 

70 

0.0070 

365657 

73.71 

70.00 

0.0195 

26.28 

6970.80 

72.55 

80 

0.0080 

411516 

82.96 

80.00 

0.0228 

34.13 

7966.65 

81.05 

90 

0.0090 

454976 

91.72 

90.00 

0.0252 

42.88 

8962.50 

91.22 

100 

0.0100 

495832 

99.95 

100.00 

0.0260 

52.49 

9958.38 

100.02 
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3.  For  =  1.0;  Zq  =  0.1;  EZq  *  1.0;  Bxq  =  1.0;  Cq  =  1.0. 


Integrating  Energy  Parameters  from  0.0  to  1.0. 


I 

z(11 

JZQL 

JliHL 

JlaQL 

_E*0L 

VyQL 

Tenu>fT) 

1 

0.00 

56.90 

l.n 

1.00 

0.00 

0.00 

9.94 

1.07 

5 

0.004 

284.44 

5.53 

5.00 

0.001 

0.11 

49.71 

5.36 

10 

0.009 

568.51 

11.06 

10.00 

0.005 

0.50 

99.43 

10.71 

20 

0.019 

1133.92 

22.05 

20.00 

0.020 

2.11 

19886 

21.37 

30 

0.029 

1692.96 

32.91 

30.00 

0.045 

4.82 

298.32 

31.96 

40 

0.039 

2242.46 

43.57 

40.00 

0.077 

8.62 

397.80 

42.41 

50 

0.050 

2779.35 

53.96 

50.00 

0.115 

13.48 

497.32 

52.69 

60 

0.060 

3300.67 

64.03 

60.00 

0.155 

19.36 

596.89 

62.74 

70 

0.070 

3803.63 

73.71 

70.00 

0.194 

26.24 

696.53 

72.55 

80 

0.080 

4285.63 

82.96 

80.00 

0.228 

34.09 

796.25 

82.05 

90 

0.090 

4744.31 

91.72 

90.00 

0.251 

42.83 

896.06 

91.22 

100 

0.100 

5177.52 

99.95 

100.00 

0.260 

52.44 

996.00 

100.02 

4.  For  ^  =  1.0;  Zq  -  10.0;  Ezq  =  1.0;  Bxq  *=  1.0;  Cq  =  1.0. 


Integrating  Energy  Parameters  from  0.0  to  10.0. 


I 

_z£LL 

J2QL 

JSiJlL 

-gaQI 

USL 

BJLL 

VyOL 

TemDfl) 

1 

0.00 

2.03 

1.01 

1.01 

0.00 

0.00 

9.6x10'* 

1.01 

5 

0.405 

10.15 

5.04 

5.04 

9.0x10"* 

9.9x10"* 

4.9x10"* 

5.03 

10 

0.911 

20.30 

10.08 

10.08 

4.4x10"* 

4.7x10"* 

1.1x10"* 

10.07 

20 

1.923 

40.61 

20.16 

20.16 

2.4x10"* 

2.5x10"* 

3.4x10"* 

20.14 

30 

2.935 

60.91 

30.25 

30.24 

8.8x10"* 

85x10"* 

9.3x10"* 

30.21 

40 

3.947 

81.21 

40.33 

40.32 

2.8x10* 

2.5X10"2 

2.6xl0"2 

40^8 

50 

4.960 

101.52 

50.41 

50.40 

8.6X10"3 

7.0xl0"2 

7.0x1 0"2 

50.35 

60 

5.972 

121.82 

60.50 

60.48 

2.6xl0"2 

1.9X10"1 

1.9X10"1 

60.43 

70 

6.984 

142.13 

70.61 

70.53 

7.8xl0'2 

5.3X10"1 

5.2X10"1 

70.51 

80 

7.996 

162.50 

80.75 

80.53 

Z2X10"1 

1.46 

1.43 

80.62 

90 

9.009 

183.28 

90.85 

90.38 

5.7X101 

4.00 

3.93 

90.70 

100 

10.021 

206.51 

99.88 

99.98 

9.5X10"1 

10.00 

11.14 

100.05 
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5.  For  no  -  1.0;  z*,  =  95.0;  Ez*,  *  1.0;  Bxo  -  1.0;  Cq  =  1.0. 


Integrating  Energy  Parameters  from  0.0  to  9S.0. 


zm 

um. 

JliflL 

jlaOL 

_i*GL 

VyQL 

TemDtn 

l 

0.00 

2.01 

1.01 

1.01 

0.00 

0.00 

6.4xl018 

1.00 

5 

3.85 

10.03 

5.00 

5.00 

4.4x1  O'11 

4.1xl0*12 

4.0xl012 

5.01 

10 

8.67 

20.06 

10.01 

10.01 

5.8xl012 

2-lxlO11 

3.5xl012 

10.02 

20 

18.30 

40.12 

20.02 

20.02 

8.2x1  O'11 

1.4xl010 

8.7X1012 

20.04 

30 

27.94 

60.18 

30.03 

30.03 

-7.7xlOn 

4.7xl010 

8.3xlO*n 

30.06 

40 

37.57 

80.24 

40.04 

40.04 

9.6x1  O*11 

l^xlO*9 

2.5xl010 

40.08 

50 

47.21 

100.30 

50.05 

50.05 

lJbclO'11 

2.4xl0‘9 

5.9xl010 

50.12 

60 

56.84 

120.36 

60.05 

60.05 

-4.9X1011 

4.4xl09 

1.8xl0"9 

60.13 

70 

66.47 

140.42 

70.06 

70.06 

-7.7xl010 

8.4xl0'9 

5.5xlO"9 

70.15 

80 

76.11 

160.48 

80.07 

80.07 

-2.0x1  O'9 

l.lxlO"* 

3.8X10"9 

80.17 

90 

85.74 

180.54 

90.08 

90.08 

l.OxlO-5 

l.lxlO"4 

l.lxlO-4 

90.19 

100 

95.37 

201.24 

99.81 

99.97 

9.0x1 0‘2 

1.66 

1.84 

100.01 
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Typical  Magnetotail  Values. 


6.  a)  For  =  1.0;  Zq  =  1.0;  Ezq  -  3.64xl0'9;  Bxq  =  1.46xl0"3;  Cq  =  4.59x10  s. 


Integrating  Energy  Terms  from  0.0  to  1.0. 


I 

zm 

_UflL 

JQiHL 

jlaOL 

-BxflL 

VyQL 

1 

0.00 

95x10  s 

1.00 

1.00 

0.00 

0.00 

15xl0"3 

5 

0.04 

4.6x10** 

5.01 

5.01 

7.8X10*8 

1.3x10** 

65xl0*3 

10 

0.09 

9.3x10** 

10.02 

10.02 

35xl0*7 

5.6x10** 

1.2xl0‘2 

20 

0.19 

1.9x1 0'3 

20.04 

20.04 

1.5X10-6 

2.4x1 0"3 

2.5x1  O'2 

30 

0519 

2.8xl0"3 

30.06 

30.05 

3.5x10"® 

5.5xl0"3 

3.8xl0*2 

40 

0.39 

3.7xl0*3 

40.07 

40.07 

6.4x10*® 

9.9xl0'3 

5.1xl0'2 

50 

0.50 

4.6xl0*3 

50.08 

50.08 

1.0x10  s 

1.6xlC2 

6.4xl0*2 

60 

0.60 

5.6xl0"3 

60.08 

60.08 

1.5x10  s 

2.3x1 02 

7.9xlC2 

70 

0.70 

6.5X10"3 

70.08 

70.07 

2.2x10  s 

3.2xl02 

9.4xl02 

80 

0.80 

7.5xl03 

80.07 

80.07 

3.0X10*5 

4.2xl02 

l.lxlO1 

90 

0.90 

8.4xl03 

90.04 

90.04 

4.1x10-® 

5. 4x1  O'2 

1.3X10*1 

100 

1.00 

9.3xl0'3 

99.99 

99.99 

55x10-® 

6.7x1 0"2 

1.5X10’1 

b)  For 

U  (  I  )  -  E  +  \  *  %  +  <V  ♦  <V  • 

o  2  2  2 


Energy  Variables*  Percentage  of  Ufll. 


_JL 

zm 

JJflL 

JliY^2_ 

(%) 

(%) 

B  2/2 
(%) 

Q iil 
(%) 

1 

0.00 

959x10  s 

0.08 

0.00 

0.00 

0.50 

5 

0.04 

4.65x10** 

0.08 

1.94xl0"12 

5.08x10“® 

0.50 

10 

0.09 

9.30x10** 

0.08 

95 9xl012 

541X10*5 

0.50 

20 

0.19 

1.86xl0’3 

0.08 

4.12x1c11 

1.05x10"* 

0.50 

30 

059 

2.79X103 

0.08 

9.90x1  O'11 

545x10** 

0.50 

40 

0.39 

3.72xl0"3 

0.09 

1.88xl010 

4.48x10** 

0.50 

50 

0.50 

4.65x1c3 

0.09 

3.19X10*10 

7.18x10** 

0.50 

60 

0.60 

5.59xl0"3 

0.09 

5.04X1010 

1.06xl0"3 

0.49 

70 

0.70 

652xlC3 

0.10 

7.62X1010 

1.49xl0"3 

0.49 

80 

0.80 

7.45xlC3 

0.10 

1.12X10-9 

501xl0'3 

0.49 

90 

0.90 

8.39xl0"3 

0.11 

1.63xl0‘9 

563xlC3 

0.49 

100 

1.00 

9.33xlC3 

0.12 

251x10*° 

3.38x1  O'3 

0.49 

The  value  for  die  minimum  total  energy  of  this  model  system  is  9.33xl0'3. 


Tempfll 

4.6xl0‘5 

2.3x10** 

4.6x10** 

95x10** 

1.4xl0'3 

1.8xl0"3 

2.3xl0"3 

58xl0*3 

35xl0*3 

3.7xl0'3 

4.1xl0"3 

4.6xl0'3 


_Cn?T 

(%) 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

050 

0.49 

0.49 

0.49 

0.49 

0.49 
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Typical  Magnetotail  Values. 


7.  a)  For  -  1.0;  Zq  =  2.0;  Ez«,  -  3.64xl0*9;  Bxq  -  1.46xl0*3;  Cq  -  4.59x10  s. 


Integrating  Energy  Terms  from  0.0  to  2.0. 


1 

_zflL 

urn. 

jiiHi 

jjaOL 

_E*0L 

-5*01- 

VyOL 

TernDfl") 

1 

0.00 

9.2x10  s 

1.00 

1.00 

0.00 

0.00 

4.0x10“* 

4.6x10  s 

5 

0.08 

4.6x10“* 

5.01 

5.01 

1.6x10*® 

8.1X10*5 

2.0xl0*3 

2.3x10*® 

10 

0.18 

9.2x1c4 

10.02 

10.02 

7.4x10*® 

3.7X1C4 

4.QxlC3 

4.6x10* 

20 

0.38 

1.8xlC3 

20.03 

20.03 

3.2xl0*7 

1.6xl0‘3 

8.2xl0*3 

9.2x10“* 

30 

0.59 

2.8xl0*3 

30.04 

30.04 

7.9xl0*7 

3.6xl0*3 

1.3x1 0*2 

1.4xl0"3 

40 

0.79 

3.7xl0*3 

40.06 

40.06 

1.6X10*6 

6.7x10*® 

1.8xl0*2 

1.8xl0*3 

50 

0.99 

4.6xl0‘3 

50.07 

50.07 

2.7x10*® 

l.lxlO*2 

2.3xl0*2 

2.3xl0*3 

60 

1.19 

5.5xl0*3 

60.07 

60.07 

4.6x10*® 

1.6xl0*2 

3.0xl0*2 

2.8xl0*3 

70 

1.40 

6.5xl0'3 

70.07 

70.07 

7.3x10® 

2.3xl0*2 

3.8xl0"2 

3-2xl0*3 

80 

1.60 

7.4xl0"3 

80.06 

80.06 

1.1x10  s 

3.2xl0*2 

4.8xl0*2 

3.7xl0*3 

90 

1.80 

8.3xlC3 

90.04 

90.04 

1.8x10  s 

4.2xl0*2 

5.9xl0"2 

4.1xl0*3 

100 

2.00 

9.2xl0'3 

99.99 

99.99 

2.6x10  s 

5.6xl0*2 

7.3xlC2 

4.6xl0*3 

b)  For 

“(I)  -  £  ^*x*¥*cV  ,cV- 


Energy  Variables*  Percentage  of  Udl. 


_L 

zfl) 

Ufl) 

(%) 

j Zn. 

(%) 

Bx2/2 

Gill 

(%) 

Cn?T 

(%) 

(%) 

1 

0.00 

9.29x1  O'5 

8.7xlC* 

0.00 

0.00 

0.50 

0.50 

5 

0.08 

4.61x10“* 

8.7xlC* 

8.63xl0*14 

2.14xlC® 

0.50 

0.50 

10 

0.18 

9.22x10** 

8.8x10“* 

4.19xl0*13 

1.02X10*5 

0.50 

0.50 

20 

0.38 

1.84xl0*3 

9.1xlC* 

1.99X10*12 

4.55X10*5 

0.50 

0.50 

30 

0.59 

2.77X10*3 

9.8xlC* 

5.41xl0*12 

1.09x10** 

0.50 

0.50 

40 

0.79 

3.69xlC3 

l.lxlO*3 

1.21x1c1 1 

2.08x10-* 

0.50 

0.50 

50 

0.99 

4.61X10*3 

1.2xl0*3 

2.53xl0*n 

3.52x10-* 

0.50 

0.50 

60 

1.19 

5.53xl0*3 

1.4xlO*3 

S.llxlO*11 

5.56x10** 

0.50 

0.50 

70 

1.40 

6.46xlC3 

1.7xl0*3 

1.02xl0“10 

8.42x10“* 

0.50 

0.50 

80 

1.60 

7.39xl0*3 

2.1xl0*3 

2.05xl0“10 

1.24X10*3 

0.50 

0.50 

90 

1.80 

8.31X10*3 

2.7xl0*3 

4.14xlO“10 

1.80X10*3 

0.50 

0.50 

100 

2.00 

9.24xl0*3 

3.4xl0*3 

8.05X10*10 

2.57X10*3 

0.50 

0.50 

The  value  for  die  minimum  total  energy  of  this  model  system  is  9.24x1  O'2. 
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Typical  Magnetotail  Values. 


8.  a)  For  n*,  =  1.0;  Zq  =  3.0;  Ezq  =  3.64xl0"9;  Bxq  «  1.46xl0‘3;  Cq  *  4.59x10  s. 


Integrating  Energy  Terms  from  0.0  to  3.0. 


I 

zm 

urn 

jiiHi 

JlaOL 

.£*01 

-B*0L 

VyOI 

Temofll 

1 

0.00 

9.2x10  s 

1.00 

1.00 

0.00 

0.00 

1.4x10“* 

4.6x10  s 

5 

0.12 

4.6x10"* 

5.01 

5.01 

3.2xl0"9 

4.4x10  s 

7.3x10“* 

2.3x10“* 

10 

0.27 

9.2x10"* 

10.01 

10.01 

1.5x10"® 

2.0x10“* 

1.5xl0*3 

4.6x10“* 

20 

0.58 

1.8xl0"3 

20.02 

20.02 

6.8x10"® 

8.6x10“* 

3.1xl0‘3 

9.2x10“* 

30 

0.88 

2.8xl0"3 

30.03 

30.03 

1.8xl0"7 

2.0xl0"3 

4.9xl0"3 

1.4X10"3 

40 

1.18 

3.7xl0"3 

40.04 

40.04 

3.9x1 0"7 

3.9xl0‘3 

7.3xl0"3 

1.8xl0"3 

50 

1.49 

4.6X10*3 

50.05 

50.05 

7.9xl0"7 

6.5xl0"3 

l.OxlO"2 

2.3xl0"3 

60 

1.79 

5.5xl0"3 

60.06 

60.06 

1.5x10"® 

l.OxlO*2 

1.4xl0"2 

2.8xl0"3 

70 

2.10 

6.4xl0"3 

70.06 

70.06 

2.8x10"® 

1.5xl0"2 

1.9xl0"2 

3.2xlO"3 

80 

2.40 

7.4xl0"3 

80.06 

80.06 

5.3x10"® 

2.2xl0"2 

2.6xl0"2 

3.7xl0"3 

90 

2.70 

8.3xl0"3 

90.04 

90.04 

9.8x10"® 

3.1xl0"2 

3.6xl0"2 

4.1xl0"3 

100 

3.00 

9.2xl0‘3 

99.99 

99.99 

1.7X10"5 

4.4xl0"2 

4.9xl0"2 

4.6xl0*3 

b)  For 

U(I)  -  E 

0  2  2  2 


Energy  Variables'  Percentage  of  Ufll. 


I 

zfll 

um 

jaiY^2_ 

(%) 

E^2_ 

(%) 

Bx^2_ 

(%) 

.ShL 

(%) 

Cn2T 

(%) 

1 

0.00 

9.21x10  s 

1.1x10“* 

0.00 

0.00 

0.50 

0.50 

5 

0.12 

4.60x10“* 

1.1x10"* 

3.34xl0"ls 

6.32xl0*7 

0.50 

0.50 

10 

0.27 

9.21x10“* 

1.2x10“* 

1.67X10"1* 

3.04x10"® 

0.50 

0.50 

20 

0.58 

1.84xl0"3 

1.3xl0“* 

8.96X10"1* 

1.39x10"® 

0.50 

0.50 

30 

0.88 

2.76xl0"3 

1.5x10“* 

2.94xl0"13 

3.51xl0"® 

0.50 

0.50 

40 

1.18 

3.68xl0"3 

1.9x10“* 

8.54xl0"13 

7.17X10"5 

0.50 

0.50 

50 

1.49 

4.60X10"3 

2^x10"* 

2.43X10"12 

1.33x10“* 

0.50 

0.50 

60 

1.79 

5.53xl0"3 

3.5x10“* 

6.94X10"12 

2.33x10"* 

0.50 

0.50 

70 

2.09 

6.45xl0"3 

5.1x10“* 

2.02x1 0"11 

3.99x10“* 

0.50 

0.50 

80 

2.40 

7.37xl0"3 

7.9x10“* 

5.96X10"11 

6.74x10“* 

0.50 

0.50 

90 

2.70 

8.30xl0"3 

1.3xl0"3 

1.79xlO"10 

1.13xl0"3 

0.50 

0.50 

100 

3.00 

9.22X10"3 

2.0xl0"3 

5.28X10"10 

1.91xl0"3 

0.50 

0.50 

The  value  for  the  minimum  total  energy  of  this  model  system  is  9.22xl0'3. 
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Typical  Magnetotail  Values. 


9.  a)  For  n<j  =  1.0;  Zq  =  4.0;  Ezq  =  3.64xl0"9;  Bxq  =  1.46xl0"3;  C0  =  4.59x10  s. 


Integrating  Energy  Terms  from  0.0  to  4.0. 


I 

UHL 

jiiOl  n2nv 

J*QL 

VyOL 

Temofll 

1 

0.00 

9.2x10  s 

1.00 

1.00 

0.00 

0.00 

5.3xl0"s 

4.6xl0'5 

5 

0.16 

4.6x10"* 

5.00 

5.00 

5.8xlO"10 

2.2x10  s 

2.7X10"4 

2.3x10"* 

10 

0.36 

9.2x10"* 

10.01 

10.01 

2.7xl0"9 

9.8x10  s 

5.4x10"* 

4.6x10"* 

20 

0.77 

l.SxlO'3 

20.02 

20.02 

1.3x10"® 

4.3X10"4 

1.2xl0"3 

9.2x10"* 

30 

1.17 

2.8xl0"3 

30.03 

30.03 

3.9x10"® 

l.lxlO"3 

2.0xl0"3 

1.4xl0"3 

40 

1.58 

3.7xl0'3 

40.03 

40.03 

9.8x10"® 

2.1xl0"3 

3.1xl0"3 

l.SxlO"3 

SO 

1.98 

4.6x1  O'3 

50.04 

50.04 

2.3xl0"7 

3.6xl0"3 

4.8xl0"3 

2.3xl0"3 

60 

2.39 

5.5xl0"3 

60.05 

60.05 

5.3xl0"7 

6.0xl0"3 

7.2xl0"3 

2.8xl0"3 

70 

2.79 

6.4xl0"3 

70.05 

70.05 

1.2x10"® 

9.6xl0"3 

l.lxlO'2 

3.2xl0"3 

80 

3.20 

7.4xl0"3 

80.05 

80.05 

2.7x10"® 

1.5xl0'2 

1.6xl0"2 

3.7xl0"3 

90 

3.60 

8.3xl0"3 

90.04 

90.04 

6.1X10"6 

2.3xl0"2 

2.5xl0'2 

4.1xl0'3 

100 

4.00 

9.2xl0'3 

99.99 

99.99 

1.3x10  s 

3.5xl0*2 

3.7xl0"2 

4.6xl0"3 

b)  For 


U(I)  - 

1  it1 

E 

0  2 

2  2 

♦  CStt,T  +  CnjT 

• 

Energy  Variables’  Percentage  of  Ufll. 

I 

zfll 

JJflL 

_n,Yy!Z2. 

3te. 

Cn,T 

-QI2L 

(%) 

c%) 

c%) 

(%) 

(%) 

1 

0.00 

9.20X10"5 

1.5x10  s 

0.00 

0.00 

0.50 

0.50 

5 

0.16 

4.60x10"* 

l.SxlO"5 

1.09X10"16 

1.51X10"7 

0.50 

050 

10 

0.36 

9-20x10"* 

1.6xlO"5 

5.65xlO"1® 

7.35xl0"7 

0.50 

0.50 

20 

0.77 

1.84xl0"3 

1.9xl0"5 

3.56xlO"ls 

3.50x10"® 

0.50 

0.50 

30 

1.17 

2.76X10"3 

2.5xlO's 

1.51X10"14 

9.47x10"® 

0.50 

0.50 

40 

1.58 

3.68xl0"3 

3.7xlO'5 

6.07X10"1* 

2.13x10  s 

0.50 

0.50 

50 

1.98 

4.60X10"3 

6.0xl0"5 

2.49X10"13 

4.43xl0"5 

0.50 

0.50 

60 

2.39 

5.52xl0"3 

1.0x10“* 

1.05xl0"12 

8.96x10  s 

0.50 

0.50 

70 

2.79 

6.44X10"3 

2.0x10"* 

4.56xl0"12 

1.80x10"* 

0.50 

0.50 

80 

3.20 

7.36xl0"3 

3.8x10"* 

2.01x10" 11 

3.60x10"* 

0.50 

0.50 

90 

3.60 

8.29xl0"3 

7.5x10“* 

9.04x10"” 

7.27x10"* 

0.50 

0.50 

100 

4.00 

9.22xl0"3 

1.5xlO"3 

4.15xlO"10 

1.48xl0"3 

0.50 

0.50 

The  value  for  the  minimum  total  energy  of  this  model  system  is  9.22xl0*3. 
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Typical  Magnetotail  Values. 


10.  a)  For  n<j  =  1.0;  Zq  =  5.0;  EZq  =  3.64xl09;  Bxq  =  1.46xl0*3;  q,  =  4.59x10  s. 


Integrating  Energy  Terms  from  0.0  to  5.0. 


I 

JJflL 

_n,ILL 

JhOL 

EzOL 

VyOL 

1 

0.00 

9.2x10  s 

1.00 

1.00 

0.00 

0.00 

1.9x10  s 

5 

0.20 

4.6x10"* 

5.00 

5.00 

9.7x1 0"n 

9.9x10"* 

9.8x10  s 

10 

0.46 

9.2x10"* 

10.01 

10.01 

4.7xl010 

4.5x10  s 

2.0x10"* 

20 

0.96 

1.8xl0*3 

20.02 

20.02 

2.5xl0*9 

2.0x10"* 

4.5x10"* 

30 

1.47 

2.8xl0"3 

30.03 

30.03 

8.4xl0"9 

5.1x10"* 

8.2x10"* 

40 

1.97 

3.7xl0"3 

40.03 

40.03 

2.5x10'* 

l.lxlO"3 

1.4xl0"3 

50 

2.48 

4.6xl0'3 

50.04 

50.04 

6.9x10"* 

2.0xl03 

2.3xl0'3 

60 

2.99 

5.5xl0'3 

60.05 

60.05 

1.9xl0'7 

3.5xl0"3 

3.9xl03 

70 

3.49 

6.4xl0*3 

70.05 

70.05 

5.3xl0'7 

6.1xl0'3 

6.5xl0'3 

80 

4.00 

7.4x1 03 

80.05 

80.05 

1.5x10-* 

l.OxlO"2 

l.lxlO"2 

90 

4.50 

8.3xl0"3 

90.04 

90.04 

4.0x10-* 

1.7xl0'2 

1.8xl0"2 

100 

5.01 

9.2xl0"3 

99.99 

99.99 

1.0x10  s 

2.9xl0"2 

3.0xl0"2 

b)  For 

u<n  -  E  ^ 

0  2  2  2 


Energy  Variables'  Percentage  of  UCD. 


I 

z(l) 

UfLL 

JllY^2. 

(%) 

E%2_ 

(%) 

3x212 

(%) 

CnjT 

(%) 

1 

0.00 

9.20x10  s 

2.1x10-* 

0.00 

0.00 

0.50 

5 

0.20 

4.60x10"* 

2.1x10-* 

3.21X101® 

3.19x10-® 

0.50 

10 

0.46 

9.20x10"* 

2.2X10"* 

1.72xl0"17 

1.57xl0*7 

0.50 

20 

0.96 

1.84xl0'3 

2.8x10-* 

1.33xl0‘16 

7.89xl0"7 

0.50 

30 

1.47 

2.76xl0'3 

4.4x10-* 

7.64X101* 

2.33x10-* 

0.50 

40 

1.97 

3.68xlO'3 

8.0x10-* 

4.45xl0ls 

5.93x10-* 

0.50 

50 

2.48 

4.60xl0"3 

1.6X10-5 

2.69xl014 

1.42X10*5 

0.50 

60 

2.99 

5.52xl0'3 

3.6xl0'5 

1.70xl0'13 

3.39x10  s 

0.50 

70 

3.49 

6.44x1  O'3 

8.4x10  s 

l.lOxlO'12 

8.13x10  s 

0.50 

80 

3.99 

7.36xlO'3 

2.0x10"* 

7.30xl012 

1.97x10"* 

0.50 

90 

4.50 

8.28X10'3 

4.9x10"* 

4.90X1011 

4.84x10"* 

0.50 

100 

5.00 

9.21xl0'3 

1.2xl0"3 

3.54xl010 

1.20xl0‘3 

0.50 

The  value  for  the  minimum  total  energy  of  this  model  system  is  9-20xl0"3. 


Temp(T) 

4.6x10  s 
2.3x10"* 
4.6x10"* 
9.2x10"* 
1.4xl0'3 
1.8xl0'3 
2.3xl0"3 
2.8xl0"3 
3.2xl0'3 
3.7xl0‘3 
4.1xl0'3 
4.6xl0'3 


Cn2T 

(%) 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 
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Typical  Magnetotail  Values. 


11.  a)  For  =  1.0;  Zq  =  10.0;  Ezq  =  3.64xl0'9;  Bxq  =  1.46xl0'3;  C0  =  4.59x10  s 


Integrating  Energy  Terms  from  0.0  to  10.0. 


I 

zfl) 

_UQL 

JhflL  n2m 

_i*QL 

-BxHL 

vyffi_ 

1 

0.00 

9.2x10  s 

1.00 

1.00 

0.00 

0.00 

1.3x1  O'7 

5 

0.40 

4.6x10"* 

5.00 

5.00 

3.8xl013 

1.3xl0'7 

6.7xl0'7 

10 

0.91 

9.2x10"* 

10.01 

10.01 

4.1xl013 

6.3x1  O'7 

1.5X10-6 

20 

1.92 

1.8xl0"3 

20.02 

20.02 

1.7xl012 

3.4x10"* 

4.6x10* 

30 

2.94 

2.8xl0'3 

30.03 

30.03 

4.7xl012 

l.lxlO"5 

1.3x10  s 

40 

3.95 

3.7xl0'3 

40.03 

40.03 

2.0x10'" 

3.3xl0"5 

3.5xl0's 

50 

4.96 

4.6xl0'3 

50.04 

50.04 

2.1xl010 

9.4xl0'5 

9.6x10  s 

60 

5.97 

5.5xl0'3 

60.05 

60.05 

1.7xl0*9 

2.6x10"* 

2.6x10"* 

70 

6.98 

6.4xl0'3 

70.05 

70.05 

1.3x10-* 

7.2x10"* 

7.2x10"* 

80 

8.00 

7.4xl0'3 

80.05 

80.05 

9.6x1  O'* 

2.0xl03 

2.Qxl03 

90 

9.00 

8.3xl0'3 

90.04 

90.04 

7.1xl0'7 

5.5x1  O'3 

5.5x1  O'3 

100 

10.02 

9.2xl0'3 

99.99 

99.99 

5.2x10* 

1.5xl0'2 

1.5xlC2 

or 

U(I) 

I 

-  E 
0 

-V. 

2 

E,2  +  B,2 
2  2 

+  Cn,T  +  Cn2T  . 

Energy  Variables’  Percentage  of  U(T). 


I 

JJQL 

JliYy^2_ 

(%) 

3%2_ 

(%) 

(%) 

CnjT 

(%) 

1 

0.00 

9.20x10  s 

9.1x10'" 

0.00 

0.00 

0.50 

5 

0.40 

4.60x10"* 

9.6x10'" 

7.88X10'23 

5.80xl0'12 

0.50 

10 

0.91 

9.20x10"* 

1.2xlO"10 

1.91X10'21 

3.18x10'" 

0.50 

20 

1.92 

1.84xl0'3 

3.4xlO'10 

4.78xl0'21 

2.45xlO'10 

0.50 

30 

2.94 

2.76X10'3 

l.SxlO'9 

5.72X10'21 

1.42X10'9 

0.50 

40 

3.95 

3.68xl0'3 

8.4xl0"9 

1.54X10'20 

8.26xlC9 

0.50 

50 

4.96 

4.60xl0'3 

5.1x10-* 

5.44xl0'19 

5.03x10-* 

0.50 

60 

5.97 

5.52xl0'3 

3.2xl0"7 

2.69X10'17 

3.18xl0'7 

0.50 

70 

6.98 

6.44X10'3 

2.1x10-* 

1.33xl0'15 

2.07x10-* 

0.50 

80 

8.00 

7.35xl0'3 

1.4xl0'5 

6.84xl0'14 

1.37x10  s 

0.50 

90 

9.01 

8.28xl0'3 

9.3x10  s 

3.91X1012 

9.24x10  s 

0.50 

100 

10.02 

9.20xl0'3 

6.3x10"* 

2.68xlO'10 

6.30x10"* 

0.50 

The  value  for  the  minimum  total  energy  of  this  model  system  is  9.20xl0"3. 


Tempfll 

4.6x10  s 
2.3x10"* 
4.6x10"* 
9.2x10"* 
1.4xl0"3 
1.8xl0'3 
2.3xl0'3 
2.8xl0"3 
3.2xl0"3 
3.7xl0'3 
4.1xl0"3 
4.6xl0'3 


Cn?Y 

(%) 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 


APPENDIX  D:  Bibliography 


127 


BIBLIOGRAPHY 


Adkins,  C.  J.,  Equilibrium  Thermodynamics,  Second  Edition,  McGraw-Hill  Book  Company  (UK) 
Limited,  Maidenhead,  England,  1968. 

Akasofu,  S.  I.,  Polar  and  Magnetospheric  Substorms,  Springer-Verlag,  New  York,  1968. 

Akasofu,  S.  I.  and  L  J.  Lanzerotd,  "The  Earth’s  Magnetosphere,"  Physics  Today,  28-35,  1975. 

Akasofu,  S.  I.,  "Concluding  Remarks,"  in  Magnetospheric  Physics,  edited  by  Lui,  A.  T.  Y.,  John 
Hopkins  University  Press,  Baltimore,  1987. 

Axford,  W.  I.,  "Tail  of  die  Magnetosphere,"/.  Geophys.  Res.,  70,  1231,  1965. 

Bim,  J.,  "The  Boundary  Value  Problem  of  Magnetotail  Equilibrium,"/.  Geophys.  Res.,  97(A11), 
19,441-19,450,  1991. 

Burkhart,  G.  R.,  J.  F.  Drake,  P.  B.  Dusenbery,  and  T.  W.  Speiser,  "A  Particle  Model  for 
Magnetotail  Neutral  Sheet  Equilibria,"/.  Geophys.  Res.,  97(A9),  13,799-13,815,  1992. 

CarovQlano,  R.  L.,  J.  F.  McClay,  and  H.  R.  Radoski,  "Physics  of  the  Magnetosphere,"  Springer- 
Verlag  New  York  Inc.,  New  York,  1967. 

Cattell,  C.  A.,  C.  W.  Carlson,  W.  Baumjohann,  and  H.  Luhr,  "The  MHD  Structure  of  the 
Plasmasheet  Boundary:  (1)  Tangential  Momentum  Balance  and  Consistency  With  Slow  Mode 
Shocks,"/.  Geophys.  Res.  Ltrs.,  19(20),  2083-2086,  1992. 

Clauer,  C.  R.  and  R.  L.  McPherron,  "Mapping  the  Local  Time-Universal  Time  Development  of 
Magnetospheric  Substonns  Using  Mid-latitude  Magnetic  Observations,"  /.  Geophys.  Res.,  79(19), 
2811-2820,  1974. 

Cheng,  C.  Z.,  "Magnetospheric  Equilibrium  With  Anisotropic  Pressure,"  /.  Geophys.  Res.,  97(A2), 
1497-1510, 1992. 

Cowley,  S.  W.  H.  and  D.  J.  Southwood,  "Some  Properties  of  a  Steady-State  Geomagnetic  Tail," 

/.  Geophys.  Res.  Ltrs.,  7(10),  833-836,  1980. 

Eastman,  T.  E.,  L.  A.  Frank,  and  C  Y.  Huang,  "The  Boundary  Layer  as  the  Primary  Transport 
Regions  of  die  Earth’s  Magnetotail,’  /.  Geophys.  Res.,  90,  9541, 1985. 

Erickson,  G.  M.,  "A  Quasi-Static  Magnetospheric  Convection  Model  in  Two  Dimensions,”  /. 
Geophys.  Res.,  97(A5),  6505-6522,  1992. 

Erickson,  G.  M.  and  R.  A.  Wolf,  "Is  Steady  Convection  Possible  in  the  Earth’s  Magnetotail?"  /. 
Geophys.  Res.  Ltrs.,  7(11),  897-900, 1980. 


128 


Fairfield,  D.  H.,*Magnetotail  Energy  Storage  and  the  Variability  of  the  Magnetotail  Current 
Sheet,*  in  Magnetic  Reconnection  in  Space  and  Laboratory  Plasmas,  Geophysical  Monograph  30, 
edited  by  Hones,  E.  W.  Jr.,  168-176,  American  Geophysical  Union,  Washington  D.C.,  1984. 


Fairfield,  D.  H.,  "Structure  of  the  Geomagnetic  Tail,*  in  Magnetospheric  Physics,  edited  by  Lui,  A. 
T.  Y.,  Chapter  2,  23-33,  John  Hopkins  University  Press,  Baltimore,  1987. 

Goldstein,  H.,  Classical  Mechanics,  Second  Edition,  Addison-Wesley  Publishing  Co.,  Reading, 
Massachusetts,  1980. 

Gosling,  J.  T.,  D.  N.  Baker,  and  E.  W.  Hones,  Jr.,  'Journeys  of  a  Spacecraft,*  Los  Alamos  Science, 
33-53,  1984. 

Hau,  L.  N.  and  G.  H.  Voigt,  "Loss  of  MHD  Equilibrium  Caused  by  the  Enhancement  of  the 
Magnetic  By  Component  in  Earth’s  Magnetotail,*  J.  Geophys.  Res.,  97(A6),  8707-8711,  1992. 

Heikkila,  W.  J.,  Tmpusive  Plasma  Transport  through  the  Magnetotpause,"  Geophys.  Res.  Lett.,  9, 
877,  1982 

Hill,  T.  W.  and  A.  J.  Dessler,  "Plasma  Motions  in  Planetary  Magnetospheres,"  Science,  252,  410- 
415,  1991. 

Hones,  E.  W.  Jr.,  The  Earth’s  Magnetotail,"  Scientific  American,  40-47,  1986. 

Kerr,  R.  A.,  "Probing  the  Long  Tail  of  the  Magnetosphere,"  Science,  226,  1298-1299,  1984. 

King,  J.  W.  and  Newman,  W.  S.,  Solar-Terrestrial  Physics,  Academic  Press,  Inc.,  London,  1967. 

Knecht,  D.  J.  and  B.  M.  Shuman,  "The  Geomagnetic  Field,"  in  Handbook  of  Geophysics  and  the 
Space  Environment,  edited  by  A.  S.  Jursa,  Chapter  4,  Air  Force  Geophysics  Laboratory,  United 
States  Air  Force,  1985. 

Lanzerotti,  L  J.,  "Earth’s  Magnetic  Environment,"  Sky  &  Telescope,  360-362,  1988. 

Lanzerotti,  L  J.  and  Krimigis,  S.  M.,  "Comparative  Magnetospheres,"  Physics  Today,  24-34, 1985. 

Lee,  D.  Y.  and  R.  A.  Wolf,  "Is  the  Earth’s  Magnetotail  Balloon  Unstable?"  J.  Geophys.  Res., 
97(A12),  19,251-19,257,  1992. 

Lui,  A.  T.  Y.,  "Observations  on  the  Fluid  Aspects  of  Magnetotail  Dynamics,"  in  Magnetospheric 
Physics,  edited  by  Lui,  A.  T.  Y.,  Chapter  3,  101-118,  John  Hopkins  University  Press,  Baltimore, 
1987. 

Lui,  A.  T.  Y.,  "Road  Map  to  Magnetotail  Domains,"  in  Magnetospheric  Physics,  edited  by  Lui,  A.  T. 
Y.,  Chapter  1,  3-9,  John  Hopkins  University  Press,  Baltimore,  1987. 

Lundin,  R.,  I.  Sandahl,  and  J.  Woch,  "The  Contribution  of  the  Boundary  Layer  EMF  to 
Magnetospheric  Substorms,’  in  Magnetospheric  Substorms,  Geophysical  Monograph  64,  edited  by 
Kan,  J.  R.,  Potemra,  T.  A.,  Kokubun  S.,  and  Iijima  T.,  American  Geophysical  Union,  Washington 
D.C.,  1991. 


129 


McPherron,  R.  L.,  "Physical  Process  Producing  Magnetospheric  Substorms  and  Magnetic  Storms," 
in  Geomagnetism,  Vol  4,  edited  by  J.  A.  Jacobs,  Chapter  7,  Academic  Press,  London’,  1991. 

National  Research  Council,  Solar-Terrestrial  Research  for  the  1980’s,  National  Academy  Press, 
1981. 

Ness,  N.  F.,  "Magnetotail  Research:  The  Early  Years,"  in  Magnetospheric  Physics,  edited  by  Lui,  A. 
T.  Y.,  Chapter  1,  11-20,  John  Hopkins  University  Press,  Baltimore,  1987. 

Nishida,  A.,  Geomagnetic  Diagnosis  of  the  Magnetosphere,  Springer-Verlag,  New  York,  1978. 

Nishida,  A.,  "Reconnection  in  Earth’s  Magnetotail:  An  Overview,"  in  Magnetic  Reconnection  in 
Space  and  Laboratory  Plasmas,  Geophysical  Monograph  30,  edited  by  Hones,  E.  W.  Jr.,  159-167, 
American  Geophysical  Union,  Washington  D.C.,  1984. 

Peterson,  I.,  "Startling  Tales  From  the  Magnetotail,*  Science  News,  139,  356,  1991. 

Pilipp  W.  and  G.  Morfill,  "The  Plasma  Mantle  as  the  Origin  of  the  Plasma  Sheet,"  in 
Magnetospheric  Particles  and  Fields,  edited  by  B.  M.  McCormac,  55,  D.  Reidel  Publ.  Co., 

Dordrecht,  Holland,  1974. 

Pippard,  A.  B.,  The  Elements  of  Classical  Thermodynamics,  Cambridge,  London,  1957. 

Schindler,  K.  and  J.  Bim.,  "Magnetotail  Theory,"  Space  Science  Reviews,  44,  307-355,  1986. 

Scholer,  M.,  "General  Discussion  on  Magnetotail  Kinetics,"  in  Magnetospheric  Physics,  edited  by 
Lui,  A.  T.  Y.,  Chapter  4,  257-259,  John  Hopkins  University  Press,  Baltimore,  1987. 

Sibeck,  D.  G.,  "Evidence  for  Flux  Ropes  in  the  Earth’s  Magnetotail,’  in  Physics  of  Magnetic  Flux 
Ropes,  Geophysical  Monograph  58,  edited  by  Russell,  C.  T.,  Priest,  E.  R.,  Lee,  L.  C.,  Chapter  _, 
American  Geophysical  Union,  Washington  D.C.,  1990. 

Speiser,  T.  W.,  "Kinetic  Aspect  of  Tail  Dynamics:  Theory  and  Simulation,"  in  Magnetospheric 
Physics,  edited  by  Lui,  A.  T.  Y.,  Chapter  4,  277-285,  John  Hopkins  University  Press,  Baltimore, 
1987. 

Tascione,  Thomas  F.,  "Magnetosphere,"  in  Introduction  to  the  Space  Environment,  Chapter  5,  Orbit 
Book  Company,  Malabar,  Florida,  1988. 

Tsurutani,  B.  T.  and  von  Rosenvinge,  T.  T.,  "ISEE-3  Distant  Geotail  Results,"  J.  Geophys.  Res. 

Ltrs.,  11(10),  1027-1029,  1984. 

Voigt,  G.  H.  and  R.  V.  Hilmer,  "The  Influence  of  the  IMF  By  Component  on  the  Earth’s  Magneto- 
Hydrostatic  Magnetotail,*  in  Magnetospheric  Physics,  edited  oy  Lui,  A.  T.  Y.,  Chapter  2,  91-97, 
John  Hopkins  University  Press,  Baltimore,  1987. 

Walker,  R.  J.,  T.  Ogino,  and  M.  Ashour-Abdalla,  "Simulating  the  Magnetosphere:  The  Structure 
of  the  Magnetotail,"  in  Solar  System  Plasma  Physics,  Geophysical  Monograph  54,  edited  by  Waite, 
J.  H.  Jr.,  J.  L.  Burch,  and  R.  L.  Moore,  61-68,  American  Geophysical  Union,  Washington  D.C., 
1989. 


130 


Williams,  D.  J.,  "Ring  Current  and  Radiation  Belts/  Rev.  of  Geophys.  25(3),  570-578, 1987. 

Young,  R.  and  E.  Hameiri,  “Approximate  Magnetotail  Equilibria  With  Parallel  Flow/  /.  Geophys. 
Res.,  97(A11),  16,789-16,802,  1992. 

Zwickl,  R.  D.,  L.  F.  Bargatze,  D.  N.  Baker,  C.  R.  Clauer,  and  R.  L.  McPherron,  “An  Evaluation  of 
die  Total  Magnetospheric  Energy  Output  Parameter,  UT ,  in  Magnetospheric  Physics,  edited  by 
Lui,  A.  T.  Y.,  Chapter  3,  155-161,  John  Hopkins  University  Press,  Baltimore,  1987. 


